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Preface 


Biomechanics of the urinary bladder as a subject for investigation is relatively 
young yet medical problems associated with it are as old as mankind. To unravel 
intricate mechanisms and thus develop effective treatments to bladder diseases 
requires a dynamic analysis of multifunctional integrative processes at all struc- 
tural levels. 

The predominant paradigm of the basic medical science is fundamentally 
reductionist. This way of thinking is described as “divide and conquer” and is 
rooted in the assumption that complex problems can be solved by splitting them 
into smaller, simpler, more tractable parts. Although this approach has been 
responsible for tremendous successes in the continuous revolution and evolution of 
molecular technologies and has created great expectations for application to 
improving diagnosing, treating, and preventing diseases, it is quite obvious that the 
merely reductionist methodology is a necessary, but not a sufficient, condition for 
gaining a comprehensive understanding of physiological mechanisms of interest. 
More specifically, it is critical to establish functional links between constructive 
elements and to study the biological system from the holistic perspective. 

Considerable progress has been made during the last decade in our ability to 
model different parts of the human body. With improvements in computing 
technology it has become possible to simulate and study the functions of the heart, 
the stomach, the small and large intestines, the uterus, etc., with astonishing levels 
of anatomical and physiological realism, and to assess performance under normal 
and pathological conditions. However, the systematic analysis and modeling of the 
human urinary bladder is beginning to attract the attention of scientists. 

The aim of this book is to provide an update on recent achievements in the field 
of biomechanics of the urinary bladder and to present a systematic development of 
a mathematical model of the organ as a thin soft biological shell. This book is not a 
textbook or workbook! We rather encourage a reader to treat it as a guidebook that 
opens the way to many unsolved problems (s)he can explore and addresses 
questions to think about. Therefore, this book could be useful to postgraduate 
students and researchers interested in the applications of computational 
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mathematics and solid mechanics to modern problems of biomedical engineering 
and medicine. 

A brief overview of the anatomy, physiology, and mechanics of the human 
urinary bladder with an emphasis on experimental facts required to justify 
assumptions and to formulate constructive hypotheses of models is presented in 
Chap. 1. The reader should not expect, though, to find a complete biomedical 
survey on the subject and is advised to consult special literature on the subject for 
further information. 

The current trend in mathematical modeling of biomechanics of the bladder is 
reviewed in Chap. 2. The existing models treat the organ as a thin elastic mem- 
brane under the general assumptions of geometrical and physical nonlinearity. 
However, they are limited to analysis of the effects of different shapes of the 
bladder on stress-strain distribution during the filling stage and do not offer the 
desired insight into the physiological processes which explain its function. 

The basic concepts of the theory of surfaces and thin soft shells essential to a 
subsequent understanding of the mathematical model of the bladder are discussed 
in Chaps. 3 and 4. The material is focused on correctness and depth of conceptual 
arguments without resorting to advanced mathematics so that the reader can easily 
comprehend the material. 

The following Chap. 5 is dedicated to constitutive model for the biological 
tissue. It is considered as a chemically reactive mechanical continuum. The fun- 
damental theoretical concepts and their ramifications are presented in ways that 
provide both their significance and biological validity. 

A one-dimensional morphostructurally relevant model of the detrusor fasciculus 
is developed and studied numerically in Chap. 6. Attention is given to analyzing 
the generation and propagation of the electromechanical wave along the myofiber 
under normal physiological conditions and after application of different pharma- 
cological compounds. 

In Chap. 7 a theoretical framework for the analysis of integrated physiological 
phenomena in the urinary bladder is developed. The emphasis is given to modeling 
of intramural regulatory pathways, i.e., neuronal arrangements and their inter- 
connections with the upper regulatory centers. The Hodgkin-Huxley formalism is 
adopted to describe the dynamics of electrical processes in neurons and synapses. 

In the following Chaps. 7 and 8 processes of electrochemical coupling at 
cholinergic and adrenergic synapses and modeling the spectrum of possible 
responses produced by acetylcholine and adrenaline on the fasciculus are 
analyzed. 

In Chap. 9 a mathematical formulation of the problem of competitive antago- 
nists, allosteric interactions, and allosteric modulation of competitive agonists/ 
antagonists, and numerically studied pharmacokinetics of drugs used in clinical 
practice to manage bladder dysfunction are presented. 

A mathematical model of the human urinary bladder—an integrated mechanical 
and self-controlled neuroregulatory system—as a thin soft biological shell; the 
results of stress-strain distribution, electromechanical activity, and pharmacolog- 
ical modulation of the organ are presented in Chap. 10. It is important to note that 
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the model contains numerous parameters and constants that have not yet been 
evaluated experimentally, e.g., mechanical properties of the human bladder under 
biaxial loading, constants of chemical reactions, and electrical properties of the 
detrusor. These parameters and constants were adjusted during numerical simu- 
lations to resemble the physiological and diseased states of the organ. Therefore, 
the results are not aimed to achieve an accurate quantitative representation but 
rather offer a descriptive qualitative evaluation of biomechanical processes in the 
bladder during filling and voiding. 

The last chapter, Chap. 11 focuses on the latest achievements in the field of 
mathematical modeling in urology and outlines current challenges. 

By presenting the underlying biological, chemical, and physical processes of 
bladder function within the unified context of mathematical models, we hope to have 
provided the reader with the knowledge and insights needed to solve medical 
problems related to organ dysfunction. We have done our best to make things 
straight. 

This book has depended on the assistance, advice, and encouragement of many 
people. To all helpers we are deeply grateful. Our special thanks go to my colleagues 
Dr. D. Malcolm and Ms. C. Squires for their invaluable comments and suggestions 
on the improvement of the manuscript. Our gratitude extends to our families for their 
continuous support and finally, to the staff at Springer, Heidelberg, Germany which 
has published this book. 


R. N. Miftahof 
H. G. Nam 
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Chapter 1 
The Bladder as a Dynamic System 


In order to understand the nature of things, we must begin by 
asking, not whether a thing is good or bad, [...] but of what 
kind it is? 

J. Maxwell 


1.1 Anatomy of the Human Urinary Bladder 


The human urinary bladder is a musculomembranous hollow organ located deeply in 
the pelvic cavity. Anatomically, the organ is divided into three major parts: the apex, 
the body, and the fundus (base), consisting of the trigone and the neck (Fig. 1.1). The 
fundus is imbedded in the prostate in males and in the musculofibrous tissue in 
females and is intimately attached to the internal urinary sphincter through the neck. 
The fibrous fascia endopelvina provides an additional connection between the base 
and the pelvic wall and the rectum. The entire body of the organ is enclosed in loose 
fatty tissue of the paravesicular fossa. The apex is covered by a thin stretchable 
peritoneum. The latter forms a series of folds—the false ligaments—which do not 
bear any biomechanical significance. The parts of the bladder are interconnected by 
the anterior, posterior, superior, right, and left lateral walls to form a smooth surface. 

Histomorphologically, the wall of the human urinary bladder consists of four 
layers: the mucous (urothelium), submucous, muscular, and serous layer. Details 
on the morphology can be found in many textbooks and research monographs 
(Elbadawi 1991; DeLancey et al. 2002; Campbell-Walsh 2007). Thus, only a few 
aspects which are relevant to the biomechanics of the organ are discussed here. 
The innermost urothelium (tunica mucosa) is made out of polyhedral in shape cells 
of stratified transitional epithelium, including basal cells, intermediate cells, and 
umbrella cells. The outer umbrella cell layer interfaces with urine and forms the 
primary barrier that includes a mucin/glycosaminoglycan layer, which may pre- 
vent bacterial attachment and diffusion of urine components across the epithelium, 
and an apical plasma membrane with low permeability to urea and water. In 
addition, umbrella cell tight junctions form a tight seal between adjacent cells, and 
are comprised of multiple claudin species which regulate paracellular transport. 
The uroepithelium maintains the barrier even as the bladder undergoes cycles of 
filling and voiding. This accommodation likely reflects the ability of the highly 
wrinkled mucosal surface of the bladder to unfold, and the increases in mucosal 
surface area that result from fusion of a population of subapical discoidal/fusiform 
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Fig. 1.1 Anatomical divisions of the human urinary bladder 


vesicles with the apical plasma membrane of the umbrella cell layer. Upon 
voiding, the mucosa refolds, and the membrane added to the apical surface of the 
umbrella cells is thought to be recovered by endocytosis (Apodaca et al. 2007). 

The submucous layer (tunica submucosa) contains a large number of collagen 
and elastin fibers, myofibroblasts, and areolar tissue. Scanning electron micros- 
copy studies showed that fibrils are densely interwoven and make a loose network 
of the three-dimensional stroma for the muscular and mucous layers (Murakumo 
et al. 1995). Because of the fine arrangement of the collagen and elastin fibers, the 
tunica mucosa possesses a great level of extensibility, i.e., it folds into rugae, when 
the bladder is empty, and stretches flat when the organ is filled with urine. 

The muscular layer (tunica muscularis) is the most prominent layer and is 
comprised of three layers: internal, middle, and external. Smooth muscle fibers of 
the inner and external layers run longitudinally—from the fundus to the apex— 
while the muscle elements of the middle layer have predominantly a circumfer- 
ential orientation. The three layers together form the detrusor muscle. The muscle 
fascicles and cells are firmly covered with collagen sheaths. They provide a 
supporting connective tissue frame to the organ and guarantee its deformability. 
Elastic fibers are, on the other hand, sparse throughout the bladder wall, except for 
denser networks around the blood vessels and muscle fascicles. 

The serosa (tunica adventia) is partial—it only covers the superior and lateral 
walls—and is derived from the peritoneum. The latter contains wavy collagen 
bundles piled up in a sheet, and intercalated by clusters of adipose cells. Ultr- 
asonographic measurements of the wall thickness, h, of the human urinary bladder 
have revealed: h = 3.3 + 1.1 mm (Hakenberg et al. 2000). It remains relatively 
constant throughout different regions of the organ. 

Blood supply to the organ is by the superior, middle, and inferior vesical arteries, 
derived from the anterior trunk of the hypogastric artery. Additionally, in the male— 
the obturator and inferior gluteal—and in the female—the uterine and vaginal 
arteries, respectively—supply small visceral branches to the bladder. The venous 
drainage is through a complicated plexus that empties in the hypogastric veins. 
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The bladder receives autonomic innervation—sympathetic and parasympa- 
thetic—and somatic, mediated by pudendal nerves. The sympathetic nerves arise 
in the thoracolumbar segments of the spinal cord, whereas the parasympathetic and 
somatic nerves originate in the sacrum. Along with their efferent function, each of 
the above nerves conveys afferent signals about bladder distension. Light and 
electron microscopic studies have revealed the existence of intramural ganglia in 
the human urinary bladder (Gilpin et al. 1983; Dixon et al. 1983, 1997; Birder 
et al. 2010). These intrinsic neurons are believed to be an extension of the pelvic 
plexus and may be involved in an integrative function of bladder activity. 

The human urinary bladder, as described above, represents a multicomponent 
dynamic system. Optimal spatiotemporal arrangements among its anatomical and 
cellular/subcellular components guarantee normal function of the organ. It must be 
emphasized, though, that a variety of patterns of its behavior are determined by 
constant interactions among structural elements and any instability in their inte- 
grated function results in development of symptoms, e.g. urgency, frequency, 
incontinence, hesitancy, postmicturition dribble, which, if not corrected, may lead 
to serious pathological conditions. The following discussion summarizes the key 
features of the main components which give rise to physiological performance of 
the bladder. 


1.2 Detrusor Muscle 
1.2.1 Morphology of Smooth Muscle 


The morphostructural unit of the detrusor is the smooth muscle cell—the myocyte. 
It has a characteristic spindle-like shape measuring ~ 100—300 um in length and 
~ 5-6 um in diameter. Its cytoplasm contains a centrally located nucleus, intra- 
cellular thin «- and f-actin (~6 nm), intermediate, mainly desmin (~ 10 nm) 
(Malmqvist et al. 1991a, 1991b), and thick (~ 20-25 nm) filaments, mitochondria 
and fairly sparse elements of the sarcoplasmic reticulum. Thin «- and fi-actin 
filaments are arranged into a lattice that is attached to the cell membrane at the 
sites of dense bands (plaques). They guarantee the integrity, strength, and high 
degree of deformability of the bladder wall and provide binding sites for myosin 
thick filaments of the SMIB and SM2B types (Martin et al. 2007). 

Regularly spaced dense bands are comprised of multifunctional proteins: 
integrins, desmin, vincullin, tensin, calponin, nonmuscle f- and y-actins, and fil- 
amin (Mabuchi et al. 1997; Small et al. 1986; Small and Gimona 1997). They 
establish direct structural and functional contacts between the intracellular cyto- 
skeleton and the extracellular matrix (ECM). The anchoring plaques play the 
essential role in transmitting forces of contraction-relaxation in the tissue, and act 
as mechanosensors in gene expression signaling pathways, cell migration, growth, 
and adaptation (Geiger and Ginsberg 1991; Yamada and Geiger 1997; Zamir and 
Geiger 2001). 
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Electron microscopy and freeze fracture studies convincingly demonstrated that 
individual myocytes are interconnected by small and irregular gap junctions. 
Confocal immunofluorescence, Western blot techniques, transcriptase-PCR reac- 
tion, and in situ hybridization methods showed that they are formed mainly by the 
subunit proteins connexin-43 and 45. They provide the structural basis for cyto- 
plasmic continuity, mediate the movement of ions and small molecules, and 
support synchronization and long-range integration in the detrusor (Fry et al. 1999; 
Wang et al. 2006; John et al. 2003; Neuhaus et al. 2002; Hashitani et al. 2001, 
2004). Myocytes are arranged into smooth muscle fasciculi, ~300 + 100 um, and 
further are assembled into bundles, —1—2 mm in length. 

Immunohistochemical evidence demonstrated the presence mainly of collagen 
types I, III, and IV, elastin fibers, laminin, osteopontin, fibronectin, and integrins 
(241-3, avf3, «5p1) in the lamina propria of the normal bladder (Wilson et al. 
1996). The three-dimensional hierarchy of folding and coiling of the fine fibrillar 
matrix in concert with adhesive proteins ensures the property of the detrusor as a 
myogenic syncytium (Fry et al. 2004; Rubinstein et al. 2007). It offers crucial 
mechanical characteristics such as high compliance, even stress-strain distribu- 
tion, and coordinated phasic contractility during filling and emptying (Nagatomi 
et al. 2007; Wognum et al. 2009). In addition, continuous remodeling of the 
stromal network allows the organ to respond acutely and efficiently to prolonged 
periods of strain by adjusting its function and structure through dynamic myocyte— 
ECM interactions and altering signaling pathways (Aitken and Bágli 2009). 


1.2.2 Electromechanical Activity of the Detrusor 


The contractile apparatus of detrusor myocytes consists of thin-actin and thick- 
myosin filaments, a family of special proteins and kinases, e.g., light chain myosin, 
tropomyosin, calmodulin, A-caldesmon, calponin, myosin light chain kinase, and 
myosin phosphatase. Actin filaments are single helical coils of actin associated 
with tropomyosin and caldesmon. Myosin filaments are made out of two coil rod- 
like structure heavy chains with a globular head domain. A principal determinant 
of the dynamics of contractions is free cytosolic calcium (Ca^) that triggers the 
cyclic actin-myosin complex formation. Two types of contractions—tonic and 
phasic—are produced by the detrusor. Thus, during the late stage of bladder filling 
the muscle generates tonic contractions, and it undergoes phasic contractions 
during bladder emptying. 

Contractility of the detrusor is controlled by spontaneous and/or induced 
electrical processes. Their repertoire depends on the balanced function of 
plasmalemmal ion channels: L- and T-type Ca?^, Ca?^-activated K*, voltage- 
dependent K*, and CIT channels. The presence of L- and T-type Ca?* channels in 
the human bladder has been confirmed by electrophysiological and pharmaco- 
logical studies (Wegener et al. 2004; Tomoda et al. 2005; Kajioka et al. 2002; 
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Uckert et al. 2000; Sui et al. 2003; Elliott et al. 1996; Badawi et al. 2006; Hashitani 
and Brading 2003a, b; Hollywood et al. 2003). They are formed of five distinct 
subunits: o4, %2, f, 6, and y. The «ı-subunit contains the channel pore, voltage 
sensor, and drug binding sites, while o», fj, 6, and y-subunits modulate the chan- 
nel's permeability. L-type channels possess characteristics of long-lasting, high- 
voltage-dependent channels and ensure the main influx of extracellular calcium 
ions, Cat, during depolarization. 

Three subfamilies of T-type Ca channels, which differ in their « subunits, 
have been identified in the detrusor. They are activated at low voltage and remain 
open for a short period of time (Badawi et al. 2006; Sui et al. 2001, 2009; Perez- 
Reyez 2003). Experimental data suggest that channels are responsible for the 
generation of spikes, pacemaker activity, and play a key role in regulation of the 
frequency of phasic contractions (Meng and Cha 2009; Brading 2006; Sui et al. 
2006). 

Potassium channels constitute a superfamily of four channels: the large 
Ca?*-activated K* (BKc,), small conductance (SKc,), voltage-gated (K,), and 
ATP-sensitive (Karp) potassium channels. The BKc, channel is made of six 
transmembrane proteins. The channel's sensitivity to calcium and activity is 
regulated by phosphorylation of the pore-forming « subunit. It offers a mechanism 
whereby cyclic nucleotides and protein kinase C modulate channel function (Tian 
et al. 2008). Two types of K, channels—delayed rectifying, and rapidly inacti- 
vating—are identified. They are formed by a single unit of six transmembrane 
proteins and the pore-hairpin loop. The channels remain uncoupled at low [Ca?*] 
and switch to a calcium sensor mode with the rise in intracellular calcium. 
Together with SKc, channels, they determine the resting membrane potential, 
action potential repolarization, excitability, and muscle contractility (Herrera et al. 
2000, 2001; Hashitani and Brading 2003a, b; Layne et al. 2010; Hristov et al. 
2011). 

Although intracellular ATP-gated K* channels (KArp) have been implicated in 
regulation of the resting membrane potential and spontaneous mechanical con- 
traction in cells, their overall contribution to electrical activity is considered to be 
relatively low. The modulation and function of these channels in normal physio- 
logical and diseased conditions have not been fully studied yet. 

No Na* channels have been identified in the human detrusor muscle. 

The distinct role of Cl” channels is unclear due to the uncertainty of their 
molecular identity. Calcium-activated chloride currents have been recorded on the 
isolated detrusor smooth muscle cells. They are evoked by elevation of Ca?* and 
have distinctive biophysical properties (Hartzell et al. 2005; Chen 2005; Dutzler 
2007). However, until the structure of these channels is resolved, we can only 
speculate on their specific functional role. 

The resting membrane potential, V", of human bladder smooth muscle cells ranges 
between —55 and —38 mV. Estimated and direct measurements of the input mem- 
brane resistance and capacitance of myocytes have shown Rm œ 125 + 49 MOQ-cm? 
and Cm œ 1.0 uF/cm?, respectively (Hashitani and Brading 2003a, b; Sui et al. 2006; 
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Meng and Cha 2009). There is a controversy with regard to the existence of spon- 
taneous slow wave activity. Some authors (Brading 2006) claim that intracellular 
recording from isolated and intact strips of detrusor do not show low amplitude 
resting potential oscillations consistent with slow waves. In contrast, the traces of 
simultaneous recordings of mechanical and intracellular electrical activity in human 
detrusor smooth muscle obtained by Visser and van Mastright (1999, 2000, 2001) 
convincingly demonstrate spontaneous fluctuations of the resting membrane 
potential of amplitude ~ 8-10 mV at a wide range of frequencies: v = 0.33—25 (Hz). 
The detrusor muscle produces spontaneous action potentials (APs) or spikes of 
magnitudes ~34—46.5 (mV) and v ~ 0.07—0.28 (Hz) (Visser and van Mastrigt 
1999; Hashitani and Brading 2003a, b). They occur as single, clusters or bursts of 
3-20 APs. Each spike has a relatively constant duration, — 1.3 s, a characteristic 
slow rising phase of depolarization, — 0.6 s, and followed by a fast after hyperpo- 
larization phase, —0.7 s. Spontaneous APs are resistant to tetrodotoxin (TTX), 
caffeine, ryanodine, thapsigargin, and cyclopiazonic acid, suggesting that extrinsic 
innervation and intracellular calcium stores do not contribute to their generation. 
However, spikes are abolished by L-type Ca?* channels blockers, e.g., nifedipine, 
verapamil, or in calcium-free solutions, indicating that they are of the intrinsic 
(intramural) origin. 

The conduction velocity, v; of APs in mammals has been evaluated using the 
electromyographic mapping technique. The results showed that the maximum v, 
in the rabbit detrusor is 3 cm/s (Kinder et al. 1998), depending on the site and 
physiological status of the organ. Electrical coupling and passive cable properties 
of detrusor muscle cells from a pig bladder were studied with the two-electrode 
method (Hashitani and Brading 2003a, b). Although the results are inconclusive, it 
is possible to assume that APs have the preferred direction of propagation along 
the axis of the muscle cell over a short distance. The spread of excitation in the 
transverse direction is poor. The general harmony of observations with the ana- 
tomical structure and distribution of gap junctions suggests that the detrusor 
syncytium possesses properties of electrical anisotropy. However, no direct 
attempts to measure preferential conductivity in the human bladder have been 
carried out. 


1.2.3 Pacemaker Activity 


There is increasing experimental evidence demonstrating that myofibroblasts— 
interstitial cells (IC)—modulate spontaneous electrical activity of the bladder. 
Using methods of transmission electron microscopy, immunostaining, and c-kit 
receptor labeling, ICs have been found to be abundantly distributed immediately 
below the urothelium and between detrusor cells and smooth muscle bundles 
(Kubota et al. 2011; van der Aa et al. 2004; Klemm et al. 1999; McCloskey and 
Gurney 2002; Hashitani et al. 2001; Hashitani 2006; McHale et al. 2006). 
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According to their location, ICs are divided into there subpopulations: (1) 
boundary IC—adjacent to the boundary of the bladder, (2) intramuscular IC— 
scattered among smooth muscle cells within muscle bundles, and (3) interbundle 
IC—distributed in connective tissues. They form close connections with intra- 
mural nerves and respond positively to various chemical mediators (McCloskey 
2010). Based mainly on morphological similarities with IC of Cajal found in the 
gastrointestinal tract, where their role as pacemakers is "established", it has been 
hypothesized that myofibroblasts in the bladder act as pacemakers. This view is 
supported by evidence that application of imatinib mesylate—a selective c-kit 
antagonist—disrupted spontaneous electrical activity in the organ (Kubota et al. 
2004). On the other hand, experiments on single and groups of smooth muscle 
cells reveal that they are able to produce spontaneous discharges even without ICs. 
Moreover, it has been shown that Ca?* transients in ICs occur independently of 
those of smooth muscles even when synchronous calcium waves sweep across 
muscle bundles (Hashitani et al. 2004). Therefore, there is reason to believe that 
ICs play a role in mediating the propagation of actions potentials and not in 
providing the focus for their generation (Hashitani et al. 2004). 

A comparative analysis of behaviors of isolated cells and muscle strips from 
different regions of the bladder also suggests that the trigone myocytes may serve 
as the precursor for spontaneous electromechanical activity (Roosen et al. 2009). 
However, the concept is based on speculative assumptions about morphofunctional 
relationships and has not been fully tested experimentally. 

The generation of strong regular electrical discharges is essential for the 
development of coordinated forceful contractions in the bladder during micturi- 
tion. It is most likely achieved through the dense intramural parasympathetic 
innervation of the wall and the network of intramural ganglia rather than through 
the syncytial cable properties guaranteed by existing adherens and gap junctions. 
Thus, sildenafil, a phosphodiesterase type 5 inhibitor, suppresses spontaneous 
contractions of intact detrusor but does not affect spontaneous activity in single 
muscle bundles. Three types of spontaneous activity have been recorded: small 
amplitude membrane potentials, APs, and slow oscillatory conductance changes. 
However, no affirmative data on the dynamics of the propagation of excitation 
within the human detrusor to support or refute this view were available up to the 
time of writing the book. 


1.3 Neurohormonal Regulatory System 
1.3.1 Anatomical Considerations 


The storage and periodic elimination of urine by the bladder depend on the 
activity of smooth and striated muscles in the bladder and the urethral outlet. 
The coordination between these organs is mediated by a complex neural control 
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system and involves the brain, the spinal cord, and the major pelvic and intra- 
mural ganglia. 

A meta-analysis of positron-emission tomography and fMRI studies of the 
human and animal brains have demonstrated the existence of non-specific—‘level- 
setting’—and highly specific centers that regulate filling and voiding mechanisms. 
The former are comprised of neuron populations located in the medullary raphe 
nuclei, the locus coeruleus, and the brain stem; the latter includes the pontine 
micturition centre (PMC), the periaqueductal grey (PAG), the preoptic and caudal 
hypothalamus, thalamus, the insula, and the cortical regions, namely, the pre- 
frontal cortex, the anterior cingulated gyrus, and the supplementary motor areas. 

Neuronal tracts in the spinal cord include sympathetic, parasympathetic, and 
somatic elements. Parasympathetic and sympathetic preganglionic neurons 
(PGNs) are found in the intermediate gray matter (laminae V-VII) sacral and 
lumbar segments, respectively, of the spinal cord. They project dendrites into the 
dorsal commissure, the lateral funiculus, and the lateral dorsal horn and exhibit an 
extensive axon collateral branching that is distributed bilaterally. The somatic 
motor neurons that innervate the external urethral sphincter are located in the 
ventral horn (lamina IX) with a diverse arrangement of transverse dendrites. 

The parasympathetic and sympathetic effector nerves to the bladder exit the 
sacral spinal cord in the anterior roots S2-S4 and T10-L2, respectively, and course 
through the pelvic and hypogastric nerves to the major pelvic ganglion. Pelvic 
ganglia are unique in that both preganglionic and postganglionic parasympathetic 
and sympathetic neurons are colocalized within the same ganglion capsule. Their 
postganglionic processes extend into the bladder wall where they contact intra- 
mural ganglia. 

The ganglia are dispersed at random throughout the detrusor muscle bundles 
and the adventitial connective tissue. The ganglia in detrusor contain a variable 
number of cells ranging from one to six neurons although groups of 2-50 neurons 
have been described. The neurons possess typical multipolar characteristics and 
are interconnected via frequent axosomatic and less common axodendritic syn- 
apses to form an extensive neuronal network (Gilpin et al. 1983; Pittam et al. 1987; 
Dixon et al. 1997; Apodaca et al. 2007; Birder et al. 2010). The cells found in the 
adventitia have up to 20 nerve cells and are associated with large nerve trunks. 
These are undoubtedly representing preganglionic inputs, either sympathetic or 
parasympathetic. Although there is little evidence for sympathetic (inhibitory) 
innervation of human detrusor muscle, it has been proposed that an inhibitory 
effect is achieved through prejunctional inhibition, via activation of possibly f2- 
adrenoceptors (ARs), of nerve pulse conduction at the levels of intramural ganglia 
and/or the nerve terminal (Mattiasson et al. 1987; Andersson and Arner 2004; Fry 
et al. 2006). It has also been suggested that excitatory presynaptic modulation of 
parasympathetic neurons is mediated by serotonin via 5-HT receptors (Candura 
et al. 1996). 

The sensory signals from the urinary bladder are conveyed to the spinal cord by 
the same pelvic, hypogastric, and pudendal nerves. They comprise thin myelinaed 
AÓ fibers that are present mainly in the detrusor, and unmyelinated C-fibers that 
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originate largely from the uroepithelium and lamina propria (Morrison et al. 2006; 
Yoshimura et al. 2008). Both types of sensory fibers are diffusely and uniformly 
distributed throughout the organ. The Ao fiber endings are considered to be the 
primary mediators of the physiological sensation of organ fullness and referred to 
as “tension receptors”. In contrast, the C fibers become active at the high level of 
bladder capacity. It has been proposed that their main roles are the sensation of 
urgency and nociception (Michel and Chapple 2009). The cell bodies of the 
afferent nerve fibers are in the dorsal root ganglia (DRG) at the lower thoraco- 
lumbar and sacral spinal segments. Some of the afferent fibers synapse on intra- 
mural ganglia and it has been suggested that they complete integrative circuits 
capable of coordinating local autonomous activity of the organ (Coolsaet et al. 
1993; Drake et al. 2001, 2003). 

Recent studies indicate that the urothelium functions as a transducer of some 
chemical and physical stimuli in the bladder. Multiple receptors for acetylcholine 
(ACh), noradrenaline (NA), adenosine S5'-triphosphate (ATP), neurotrophins, 
bradykinin, purines, etc., ion mechanosensitive Na* channels, myofibroblasts, and 
primary afferent and efferent neurons have been identified near to and within the 
sub- and urothelial layers. Through the channels and receptors, the urothelium is 
able to respond to various sensory inputs with the release of neurotransmitters— 
ACh, NA, ATP, prostaglandins, prostacyclin, nitric oxide (NO), and cytokines 
among others (Birder et al. 2010). 


1.3.2 Neurotransmission 


Based on the main neurotransmitter involved in the dynamics of signal 
transduction mechanisms, neurons are classified as cholinergic (the main neuro- 
transmitter is ACh) and adrenergic NA, respectively. However, the proposed 
classification does not reflect reality because pre- and post-ganglionic fibers of 
nerve cells co-localize and co-release multiple non-cholinergic and non-adrenergic 
transmitters, e.g., 5-hydroxy-tryptamine (5-HT), NO, substance P (SP), neurokinin 
A, B (NKA, B), adenosine - 5’ triphosphate, etc. Somatic motor nerves, however, 
are predominantly cholinergic. 

ACh is the major excitatory neurotransmitter in the human urinary bladder. It is 
synthesized in neurons by the enzyme choline acetyltransferase from choline and 
acetyl-CoA. ACh molecules are stored in vesicular form in the presynaptic nerve 
terminals. Recently, it has been established that there are two main sources of ACh 
release in the bladder: a neuronal source, which originates from intrinsic neurons, 
axons, and extrinsic neurons, and a non-neuronal source, which emanates from 
urothelial cells (Yoshida et al. 2004, 2008; Zagorodnyuk et al. 2010). The main 
trigger for its release is cytosolic calcium, Ca, which is increased during 
depolarization of the nerve terminal. Free Cat binds to reactive calcium centers 
on the vesicles, causing their fusion with the presynaptic cell membrane, and 
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vesicular ACh exocytosis into the synaptic cleft. The process has a quantum 
character. The excess free fraction of ACh reacts with acetylcholinesterase enzyme 
which is abundant in the synaptic cleft. Two types of enzyme are identified—true 
(specific), and pseudo (non-specific) cholinesterases. The former exhibits the 
highest activity against the free ACh than against other choline esters, and the 
latter is responsible for controlling the effect of neurotransmitter acting as a local 
hormone. Two binding sites at the active center of the enzyme have been found: 
the anionic site which reacts with the cationic head of the ACh molecule, and the 
esteric site which combines with its ester group. Both types of enzyme catalyze the 
hydrolysis of ACh, but only true cholinesterase acts on acetyl-fj-methylcholine. 
The free fraction of ACh diffuses to the postsynaptic membrane where it binds to 
muscarinic receptors. 

Although the bladder expresses the full complement of muscarinic receptors 
(u1-u5), RT-PCR and Western blot analyses demonstrated the dominance of type 2 
(u2) and type 3 (u3), compared to 44, H4, and us, receptors (Mansfield et al. 2005; 
Sigala et al. 2002; Yamaguchi et al. 1996; Alexander et al. 2001; Kotlikoff et al. 
1999). Quantitatively, 2 type receptors outnumber the 445 type at a 3:1 ratio but the 
functional affinity is shown to be greater for u3 compared to u2 receptors 
(Chess-Williams et al. 2001). Details about their distribution within the organ are 
still emerging and the existing experimental evidence suggests two main sites of 
location: the detrusor smooth muscle and IC of suburothelium. This evidence, 
along with clinical observations, supports the current hypothesis that u3 receptors 
are mainly responsible for the normal rhythmic electromechanical activity of the 
organ (Fry et al. 2010). Their regulatory effects are achieved through direct 
stimulation of smooth muscle and ICC cells (Andersson et al. 1991, 1995; 
Andersson 2002). Non-neuronal cells of the urothelium, traditionally viewed as 
passive barriers, have been proposed to play an essential part in bladder sensory 
mechanisms. They specifically respond to chemical and mechanical stretch stimuli 
and thus offer a reciprocal chemical communication with intramural nerves and the 
detrusor (de Groat 2004; Apodaca 2004; Hanna-Mitchell and Brider 2008). 

ACh is co-released along with ATP from parasympathetic nerve varicosities. 
One molecule of ATP contains three phosphate groups, and is produced by ATP- 
synthase from inorganic phosphate and adenosine diphosphate or adenosine 
monophosphate. In the human bladder, ATP acts via purinergic P2X, type 
receptors that function as non-specific cation channels (Elneil et al. 2001; O'Reilly 
et al. 2001; Hardy et al. 2000). Hydrolysis of ATP ultimately yields adenosine, 
which exerts its own effects on the detrusor through P1 type receptors. Four 
subtypes have been described so far: Aj, A24, Aog, and A, that are coupled to 
G-proteins either positively (A2) or negatively (Aj/3). All types of receptors 
unequivocally regulate adenylyl cyclase (AC) activity. Despite its co-localization 
and co-release with ACh, the functional role of ATP in normal bladder mechanics 
is not significant. 

NA, also known as norepinephrine, is a hormone and a neurotransmitter. It is 
synthesized in the cytosol of adrenergic neurons and adrenal gland medullary cells. 
The main precursor—the amino acid tyrosine is converted in an enzymatic 
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pathway into a series of intermediates, namely L-dopa, dopamine, NA, and ulti- 
mately adrenaline. Newly synthesized molecules are stored in chromaffin granules. 
Intracellular rise of calcium triggers its release by exocytosis. NA acts via two 
main groups of G-protein coupled ARs, of « and f) types. Low concentrations of 
&14 and x;p-mRNAs ARs are present in the fundus, trigone, and base of the organ 
(Malloy et al. 1998). They are linked to the Ga); protein system and PKA 
intracellular signaling pathways. 

Using cloning, in situ hybridization, and pharmacological methods, it was 
demonstrated that the detrusor expresses three subtypes of ARs-f,, f», and fs. 
They are coupled to Ga,, and Gj, proteins and are associated with activation of the 
cAMP-dependent intracellular pathway. Selective f;5-AR agonists cause effective 
inhibition of bladder contractions. This effect cannot be achieved with the use of 
specific f, and f^ agonists indicating that mainly fjj-ARs mediate relaxation 
(Igawa et al. 1999, 2001; Takeda et al. 1999; Yamaguchi 2002). Results of recent 
in vitro investigations have suggested though that cAMP plays a minor role in f- 
adrenergic receptor-medicated relaxation (Frazier et al. 2005). Specifically, it has 
been shown that a f-AR agonist can activate potassium channels, and thus exerts 
its inhibitory effect (Kobayashi et al. 2000). 

NA exerts negative feedback to downregulate its own synthesis at the pre- 
synaptic %2 ARs. Excess neurotransmitter is removed by deamination and degra- 
dation by monoamine oxidase and catechol-O-methyltransferase enzymes, 
respectively, to inactive metabolic products. 

5'-hydroxytriptamine (serotonin) also has been suggested to play a role in 
detrusor contractility. Its is formed from the amino acid r-tryptophan which 
undergoes metabolic conversions controlled by tryptophan hydroxylase 1 and 
amino acid decarboxylase enzymes. Neurons do not synthesize 5-HT, but uptake 
exogenous neurotransmitter through high affinity uptake sites located on their 
axons and terminals. 

Serotonin has multiple sites of action, including the central nervous system, the 
dorsal and ventral horns in the lumbosacral spinal cord, and presynaptical para- 
sympathetic nerve terminals (Tonini et al. 1994; D'Agostino et al. 2006). In the 
human bladder it acts at 5-HT554 type receptors. The stimulation of 5-HT; 
receptors enhances contractility of the detrusor directly, while the activation of 
5-HT3 4 receptors triggers contraction through the neuronal mechanism. Con- 
versely, the inhibition of cholinergic neurotransmission at nerve terminals and 
relaxation of the bladder wall are controlled by 5-HT,;4 receptors (D'Agostino 
et al. 2006). Serotonin is metabolized in the liver by two-step enzymatic oxidation 
to 5-hydroxyindoleacetic acid (5-HIAA). Monoamine oxidase and aldehyde 
dehydrogenase are essential enzymes that control the process. 

Tachykinins, including substance P (SP), and neurokinins A and B (NKA, 
NKB) are present in the human bladder. Radioligand binding, autographic and 
pharmacological studies have helped identify three types of receptors—NK1, 
NK2, and NK3—in the primary afferent nerves of intramural ganglia and the 
detrusor muscle (Burcher et al. 2000; Giuliani et al. 1993; Palea et al. 1996; Smet 
et al. 1997). Tachykinins cause contractions of the organ via the excitation of NK2 
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receptors. The effect is linked to activation of L-type Ca?* channels (Warner et al. 
2003; Maggi et al. 1989; Zeng et al. 1995). As an alternative supplementary 
mechanism, the involvement of the phospholipase C (PLC) pathway is proposed. 

In addition, many neuropeptides, e.g., endothelin 1, vasoactive intestinal 
polypeptide, somatostatin, angiotensin I, calcitonin gene-related peptide, NO, 
prostaglandins (PGI, Es, F24), thromboxane A>, and neuropeptide Y among 
others, have been demonstrated to be synthesized, stored, and released in the 
organ. Their functional role has to be established but it has been hypothesized that 
their role as neuromodulators is small. Therefore, they will not be considered in 
this book. 


1.3.3 Electrophysiological Characteristics of Neurons 


Studies of the electrical properties of pelvic ganglia neurons are limited to animals, 
specimens only, mainly rats and cats. In vitro preparations, which involve the 
enzymatic dissociation of cells, and in vivo intracellular recordings from intact 
ganglia, have provided a relatively concise description of electrophysiological 
characteristics and patterns of activity of neurons (Tabatabai et al. 1986; Pittam 
et al. 1987; Yoshimura and de Groat 1996; Felix et al. 1998; Lee et al. 2002; 
Kanjhan et al. 2003; Tan et al. 2007). The resting membrane potential for most of 
the neurons changes between —47.8 + —53.4 (mV) and the input resistance: 
—35 + 68 (MQ). The cells exhibit phasic and tonic types of discharges of 
amplitude 60-90 (mV) followed by a short period of afterhyperpolarization 
(AHP). Further, phasic neurons demonstrate rapidly (<4 APs) and slow adapting 
(24 APs) responses to an injected depolarizing current. Tonic neurons fire with 
increasing frequency, 10 < v < 80 Hz, which is proportional to the stimulus 
strength, and the firing rate remains relatively stable for pulses applied at 
v < 0.1 Hz. However, discharge patterns change in both types of neurons as a 
result of increase in the frequency of excitation. Thus, the tonic cells generate APs 
of a “tonic bursting” characteristic, while the slowly adapting phasic cells produce 
spikes that resemble the “excitable” reaction in response to high frequency 
repetitive stimulation. In addition, the majority of phasic and tonic neurons in the 
major pelvic ganglion show rebound APs at the end of hyperpolarization. 

Some tonic type neurons also exhibit spontaneous activity with production of 
low ~5 mV, and high amplitude 60-90 (mV) APs. They can be blocked by 
application of TTX, suggesting this activity is not intrinsic in nature. 

Data on electrical characteristics and patterns of behavior of the urinary bladder 
intramural ganglia are even sparser (Pittam et al. 1987; Hanani and Maudlej 1995). 
There are many similarities between the pelvic and intramural ganglia neurons. 
Intracellular recordings from neurons have shown that their resting membrane 
potential ranges between —40 and —60 (mV), and the input resistance is ~58 MQ. 
The injection of depolarizing currents into cells at frequencies >0.1 Hz causes the 
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generation of either tonic or phasic spiking, similar to those observed in the major 
pelvic ganglion. APs are usually followed by the period of fast AHP. 

Intramural neurons are also capable of generating spontaneous low (1-3 mV) 
and high amplitude APs (50-60 mV) at a frequency of 60 Hz. Compared to the 
spontaneous activity recorded in pelvic ganglia, the produced spikes are not 
affected by hexamethonium, suggesting that they are intrinsic in origin. 


1.4 Morphofunctional States in the Bladder 


1.4.1 Bladder Filling 


The bladder functions—storage and elimination of urine—are based on a 
coordinated interplay of mechanical reactions of the detrusor and the outflow 
region—the internal and external sphincters. Dynamic connections among intrinsic 
electrical, chemical (neurohormonal), and mechanical processes in the urinary 
bladder are required for the organ to function as a physiological entity. Existing 
ascending and descending pathways carry information between the brain, the 
spinal cord, and major pelvic and intramural ganglia and guarantee the effective 
regulation of bladder function. These forward—feedback interactions modulate the 
expression and activity of surface membrane and nuclear receptors, extra- and 
intra-cellular signaling pathways, ion channels, and contractile proteins, and define 
a variety of patterns of behavior. Electro-chemo-mechanical coupling is a 
sequence of events that results in mechanical activity of the detrusor. It is preceded 
by a wave of depolarization of the cell membrane and/or ligand-receptor complex 
formation that triggers a sequence of intracellular processes leading to relaxation 
or contraction of the muscle. 

Storage reflexes, known collectively as the ‘guarding reflex’, are organized 
primarily in the spinal cord and activated during bladder filling. Throughout the 
phase the detrusor remains relaxed and the internal/external sphincters contracted, 
thus preventing involuntary organ emptying. Such coordination is achieved by 
activation of the intramural stretch receptors whose distension produces low- 
amplitude afferent signals that are conveyed through the pelvic nerves to the spinal 
cord. Here, they are organized by interneuronal circuitry to induce concurrent 
excitation to the sympathetic outflow in the hypogastric nerve and intramural 
ganglia of the organ. Some input from the lateral pons have been suggested to have 
a role in involuntary sphincter control. Additionally, the presumed connection 
between central and lateral PAG probably enables higher centers to control 
excitation to the PMC. During bladder filling PMC neurons and the parasympa- 
thetic pathway are turned off, but at a critical level of distension the afferent 
signaling switches them to maximal activity. 

The main neurotransmitter involved at this stage is noradrenalin. It binds 
positively to f3-Ga,, and Gj.-coupled ARs located on the bodies of 
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parasympathetic neurons and smooth muscle cells and activates the (AC) enzyme. 
The latter catalyzes the conversion of ATP into cyclic adenosine monophosphate 
which subsequently activates the protein kinase A (PKA) pathway. Binding of 
cAMP to the PKA enzyme causes the transfer of ATP terminal phosphates to 
myosin light chain kinase. As a result, there is a decrease in the affinity of myosin 
light chain kinase (MLCK) for the calcium-calmodulin complex with eventual 
detrusor relaxation. It has been hypothesized that PKA may augment the effect by 
inhibiting PLC, intracellular Ca** turnover, and by affecting directly the perme- 
ability for BKc, channels. 

PKA activity is controlled entirely by cAMP. Under low concentrations of the 
enzyme, it remains catalytically inactive. The level of cAMP is regulated both by 
the activity of AC and by phosphodiesterases that degrade it to 5'-AMP. There is 
growing experimental evidence that phosphodiesterases 1, 4, and 5 (PDEI, 4, and 
5), play a role in relaxation of the detrusor induced by parasympathetic stimulation 
(Truss et al. 1996a, b, 2001; Laurenza et al. 1989; Uckert et al. 2002; Wheeler 
et al. 2005; Oger et al. 2007; Filippi et al. 2007). Another mechanism that 
decreases the production of cAMP is the activation of Gag, proteins which may 
directly inhibit AC through the MAPK signaling. However, no experiments have 
been conducted either to confirm or reject the mechanisms. 

Another concomitant sequence of events related to the neurochemical processes 
described above is electrical hyperpolarization of the detrusor. The generation of 
the inhibitory postsynaptic potential causes an increase in the BKc,, Ky, and Cl” 
channels’ open state probability and a reduction in permeability for voltage-gated 
L- and T-type Ca?* channels with resulting hyperpolarization of the cell. Deac- 
tivation of the L- and T-type Ca? channels affects the influx of extracellular 
calcium ions inside the cell. Non-selective cation channels, Ca?^ activated C17 
channels, and Na* /ca?* exchanger do not play a major role in Ca^* influx 
(Frasier et al. 2008). However, a basal level of free cytosolic calcium, 50-100 nM, 
and the release of Ca?t from the SR through Ca?'-induced calcium release, 
ryanodine, and inositol 1,4,5-trisphosphate (IP3) receptor activation are insufficient 
to initiate or to sustain contractions. The detrusor muscle remains relaxed. 


1.4.2 Bladder Voiding 


The following micturition phase is mediated by reflexes that start in the brain and 
depends entirely on neuromyogenic activity of detrusor. Functional imaging studies 
of voiding have demonstrated that the process is associated with activation in the 
prefrontal cortex, the insula, the hypothalamus, the PAG, and PMC. The PAG 
receives inputs from sacral afferents and pass them to the PMC. The PMC and the 
paraventricular nucleus of the hypothalamus, in turn, connect non-specifically to the 
]umbosacral parasympathetic nuclei and autonomic preganglionic motor neurons in 
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the spinal cord. From there, the signal is transmitted to the parasympathetic neurons 
of intramural ganglia. 

The excitation of the parasympathetic nerves, the release of ACh (predomi- 
nantly), and the depolarization of the smooth muscle membrane result in phasic 
mechanical contractions of detrusor. Although the urinary bladder expresses dif- 
ferent types of muscarinic receptors, contractions of the normal detrusor appear to 
occur largely via stimulation of u3 muscarinic receptors. The receptors are coupled 
to a range of intracellular signaling pathways that involve Gio and Gg/;; proteins 
(Fetscher et al. 2002; Hegde and Eglen 1999). Activation of u3 receptors specif- 
ically triggers the cascade of G;, protein reactions that leads to stimulation of 
PLC with the subsequent hydrolysis of phosphoinositide, formation of inositol- 
1,4,5-triphosphate (IP3) and 1,2-diacylglycerol (DAG), mobilization of intracel- 
lular Ca?* ions from the sarcoplasmic reticulum, and, finally, contraction. Various 
other signaling mechanisms have been demonstrated to be linked to u3 receptor 
activation—Rho, tyrosine, and mitogen-activated protein kinases, and phospho- 
lipases D and A5, among others. However, results of numerous in vitro and in vivo 
studies suggest that they play a minor role in muscarinic receptor agonist-induced 
contraction of the organ (Frasier et al. 2008). Excitation of u receptors and 
stimulation the Gj) protein system causes the downstream inhibition of the cyclic 
adenosine monophosphate (CAMP) production, which is the main relay in the 
adrenergic pathway, as well as the direct inactivation of membrane voltage- 
dependent Ca?* -K* channels. The final result of the above cholinergic translations 
is tight control of detrusor contractility. 

Contractions of the bladder are comprised of localized “micro-contractions” 
that occur in single or multiple discrete regions, electrical waves, and micro- 
stretches. Tension, T = 1.0 — 2.0 (g), is recorded in the organ at different loca- 
tions—the detrusor and the bladder base (Brading 2006; Coolsaet 1985). Their 
intensity is small, compared to contractions that are evoked by stimulation of the 
nervous system or triggered by mechanical stretching, and thus cannot effect 
significant changes in the shape and intravesicular pressure, p, of the bladder. The 
evaluation of p as a function of volume, V, showed that in a man with normal 
bladder compliance p(0) — 0 — 0.8 kPa and p(300) — 1.46 — 3 kPa (numbers 
in brackets correspond to V — 0 and 300 ml, respectively) (Malbrain and Deeren 
2006; Schmidt et al. 2002), the phenomenon has been attributed partly to spon- 
taneous contractility that allows the individual muscle bundles to adjust their 
length in response to filling and partly to viscoelastic properties of the tissue. 

The average amount of urine that the bladder accommodates at a time is ~0.3- 
0.5 1. It varies from person to person and may be affected by underlying patho- 
logical conditions. Thus, the capacity in severe cases of fibrosis reduces to 0.2 1 
and in case of obstruction it may reach 2.5-3 1. The largest documented volume, 
9.35 1, was recorded in a patient with severe obstruction due to cancer of the neck 
of the bladder (Leslie 2009). 

Despite its complex relationships with the rectum, uterus, pelvic bones, blood 
vessels, and ducts, the bladder expands during urinary filling without significant 
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resistance from surrounding structures. The actual size, position, and relations of 
the bladder to other organs vary according to the amount of fluid it contains. Thus, 
the empty bladder is placed entirely within the pelvis and has the form of a 
flattened tetrahedron with its vertex tilted forward. The configuration of the 
bladder during distension has been studied extensively using computed tomogra- 
phy (CT) and magnetic resonance imaging (MRI) (Roeske et al. 1995; Miralbell 
et al. 1998; Villeirs et al. 2004; Frank et al. 2006; Pinkawa et al. 2007). However, 
because each imaging process takes several minutes, the traditional methods of 
data acquisition could only offer results which are qualitative in nature. Only 
recently, with high-resolution MRI technology combined with computer-based 
image rendering has it become possible to get a comprehensive insight into the 
dynamics of shape changes of the organ (Krywonosa et al. 2010). It has been 
demonstrated clearly that during the filling stage, the bladder can assume diverse 
forms ranging from the most commonly observed ellipsoid to various irregular 
shapes. Radii of a fully distended bladder measured in the three principal direc- 
tions, i.e., the anterior—posterior, R;, superior-inferior, R., and left-right, R,, are 
R, = 6 cm, R: c 4.3 cm, and R, ~ 4.3 cm, respectively. 

Voiding begins when a critical level of bladder distension is reached and strong 
afferent signals arising from stretch receptors switch the nervous system to max- 
imal activity. In humans the process is voluntary. It involves multiple functional 
areas in the brain and the rostral brainstem, and implies an interruption of the tonic 
suppression of the PAG cells input to the PMC. Excitation passes from the PMC 
along intraspinal descending pathways to the sacral segment where it activates the 
parasympathetic outflow to the bladder. Stimulation of cholinergic pathways 
results in contraction of the detrusor. During micturition, force generation and 
muscle shortening are initiated relatively quickly and highly synchronized, and 
occur over a large length range. This efficiency is mainly achieved through the 
dense parasympathetic intramural network and intracellular mechanisms—IP3 
production and Ca?^ sensitization (Fry et al. 2004; Takahashi et al. 2004). A 
transient rise in [Ca2*] to 1 uM and the subsequent formation of a Ca?* —cal- 
modulin—active MLCKcomplex induce the myosin phosphorylation and con- 
traction. Many excellent reviews on the dynamics of smooth muscle mechanics are 
available in the literature and the interested reader is advised to consult them 
directly for further details. Strong contractions with patent sphincters cause an 
increase in intravesicular pressure, p — 7.8 kPa and a change in configuration of 
the bladder: it becomes more rounded. With relaxation of the sphincter and urine 
outflow, the pressure falls to the resting level and the organ attains the undeformed 
configuration. 
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Chapter 2 


Investigations into Biomechanics 
of the Bladder 


At each stage entirely new laws, concepts, and generalizations 
are necessary, requiring inspiration and creativity to just as 
great a degree as in the previous one. 

P. Anderson 


2.1 Biomechanics of the Detrusor 


A knowledge of the mechanical properties of the tissue of the bladder wall is crucial 
for the integration of motor functions into a biologically plausible biomechanical 
model. The combined study of urine flow in the bladder, the urinary sphincter, and the 
urethra is called urodynamics. The cystometric technique was first reported by Rose 
(1927) and has been extensively used since then for both clinical and research 
purposes. Today, urodynamic investigations remain “the gold standard” and the 
most definitive tests for the evaluation of organ function objectively through a 
physiological micturition cycle. They comprise a series of tests which imply the real- 
time monitoring of changes in bladder volume, intravesical Pyes, and abdominal Paba, 
pressures, uroflowmetry, surface electromyography, and video-urodynamics during 
subsequent phases of artificial bladder filling and emptying (Schmidt et al. 2002; 
Nitti 2005). All recordings are made via a urethral catheter and a rectal balloon and 
acquired data are used to calculate the cumulative bladder capacity V, detrusor 
pressure paei(Daet = Pves — Paba), and its contraction strength T, urine outflow rate Q, 
and urethral opening pressure and resistance. Although the method provides valuable 
quantitative information about the overall behavior of the bladder one has to realize 
that these studies are restricted to the assessment of micturition parameters only 
which are dependent essentially on theoretical concepts—physiological, mathe- 
matical, computational—employed in the evaluation of functionality of the organ. 
As a part of urodynamic studies detrusor electromyography (EMG) focuses 
primarily on recording electrical smooth muscle activity in living animals during 
active bladder emptying using surface electrodes, attached to the bladder wall, or 
needle electrodes, directly inserted into the detrusor muscle. This method is being 
utilized for fundamental studies of bladder muscle physiology, testing the integrity 
of neural circuits, and intrinsic/extrinsic control mechanisms under normal and 
pathological conditions (Kinder et al. 1997; 1999, 2001; Ballaro et al. 2001). EMG 
recordings offer higher spatial resolution and better dynamic estimates of bladder 
function than the assessment of pressure changes only. Thus, spontaneous slow 
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wave and repetitive spiking activities originating from the detrusor itself have been 
demonstrated convincingly in in situ experiments on rabbit bladders (Kinder et al. 
1997, 2001). The multiple electrode technique positioning along the wall has 
allowed estimation of the maximum conduction velocity of an arbitrary spike, 
which is —3 cm/s. The time domain and power spectrum analyses of recorded 
data clearly indicate the correlation of detrusor electromyogenic activity with 
intravesical pressure rise (Scheepe et al. 1999). Despite significant advancements 
in animal EMG studies, no electromyographic recordings have been obtained and 
verified from human detrusor. 

Although a large number of urodynamic studies on isolated and intact animal 
and human bladders have been carried out to measure their electromechanical 
activity and an enormous number of publications are available in the open liter- 
ature, these results do not allow reconstruction of active and passive uniaxial and 
biaxial force-stretch ratio relationships, evaluation of structural changes in the 
detrusor and surrounding tissues, and assessment of the spatiotemporal dynamics 
of variations in mechanical properties during the tension bearing process. There- 
fore, it is not surprising that credible models of the bladder tissue as a mechanical 
biologically active continuum have not been constructed yet. 

Most experiments on the bladder tissue under simple and complex loading pro- 
tocols have been conducted on animals and only a few studies have been dedicated to 
the investigation of the human organ per se (Alexander 1971, 1976; Kondo et al. 
1972; Coolsaet et al. 1975a, b, 1976; van Mastright et al. 1978; van Mastrigt and 
Nagtegaal 1981; Andersson et al. 1989; Venegas et al. 1991; Wagg and Fry 1999; 
Finkbeiner 1999; Sacks 2000; Gloeckner et al. 2002; Gloeckner 2003; Korossis et al. 
2009; Nagatomi et al. 2008; Parekh et al. 2010; Wognum 2010). Linear strips for 
uniaxial stretching were usually collected from different regions of the organ. Since 
the experiments were performed on the segments removed from the host, it was 
assumed that the muscle fibers were fully relaxed and the mechanical contribution 
was attributed to mechanochemically inert components of smooth muscle cells 
along with elastin and collagen fibers. In vitro quasi-static and dynamic tension tests 
were performed along two structurally defined orthogonal directions of anisot- 
ropy—the longitudinal (4j) and circumferential (4.). Their orientation coincided 
with the long and circumferential axes of the bladder, respectively. Assuming the 
homogeneity of the stress and strain fields and the incompressibility of the tissue, the 
passive force and stretch ratios p 17 Ac,1) were calculated. The interpolation of 


data in the preferred axes of structural anisotropy yield 
Trey = erlexper(aey 1) — 15 Aen >], (2.1) 


where c;, c2 are mechanical constants. 

Experimental results have demonstrated that the tissue has nonlinear, pseudo- 
elastic properties and is similar to other biological materials. Analysis of the 
T (ĉc) curves has shown a characteristic “triphasic” response with a nonlinear 
transition between the low and high elastic states (Fig. 2.1). The bladder tissue is 
compliant at low levels of stretching (phase 1) followed by a highly nonlinear 
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Fig. 2.1 A representative 
force-stretch ratio curve for a 
soft biological tissue 


1.0 A 


transitory phase 2. Specimens demonstrate pure linear elastic behavior for high 
levels of stretching (phase 3). The representative data for the pig, female rat, and 
dog urinary bladders are given in Table 2.1. 

Overall analysis across different species has revealed that the bladder wall has a 
considerable inherent inhomogeneity in its material properties, and it does not 
stretch equally in all directions—it is more compliant circumferentially than 
longitudinally. It is noteworthy that while insignificant differences between the 
loading and unloading curves are present due to “biological hysteresis”, the force— 
stretch ratio responses are independent of the stretching rate. 

Histoarchitectural correlations with the dynamics of stress-strain development 
in the bladder wall have revealed that uncoiling of ECM collagen fibers and small 
randomly oriented crack growth already begins at early stages of filling of the 
organ. They steadily increase in size as the distension of the bladder proceeds 
(Fig. 2.2a, b). There is a disrupture in dense packaging of the fibrillary—collagen 
and elastin—matrix with the expansion and confluence of multiple small fractures. 
The distribution and orientation of elastin fibers in the bladder wall is region and 
direction dependent. Most of the elastin is present in the ventral and lateral regions 
and appear to be oriented predominantly circumferentially. The detrusor muscle 
and collagen fibers though are most compact within the lower body and trigone 
regions and are the least affected by distension. 

Viscoelastic properties of the bladder wall tissue were studied extensively on 
uniaxially loaded strips in vitro and whole organ in vivo. The ramp and quasi-static 
loading protocols were employed in experimental settings. The quasi-linear vis- 
coelastic model was used to describe the strain history dependence and hysteresis 
(Fung 1993). It assumed that the relaxation function K(4, t) is the product of the 
pseudoelastic response T(A) and a reduced relaxation function G(r) 


XU ems I T[A(t — 1)] Cts) Or (2.2) 


where 


| do ea[X(t/12) — X(t/21)] 
G(t) — kealia] for tı <t<1, (2.3) 
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Table 2.1 The representative data for the pig, female rat, and dog urinary bladders 


Species Phase 1 Phase 2 Phase 3 
Pig 10 < 4,4 < 1.25 12354 < 1.8 Imax ~ 3 0.2 

T^, ~ 0 + 0.8 kPa T?, ~ 0.8 + 32 kPa Th. = 4 + 0.5 kPa 
Female rat 1.0 < À., < 1.15 1.15 <A.) < 1.28 Amax ~ 1.3 0.1 

T^, ~ 0 + 10 kPa T^, ~ 10 + 65 kPa T? x = 100 + 25 kPa 
Dog 1.0 < 21 < 1.05 1.05 < day < 1.15 Imax ~ 1.6 + 0.2 

T^, ~ 0 + 38 kPa T^, ~ 38 + 82 kPa TP, = 140 + 20 kPa 
and 


X(t/t) = Je "'ft)dt where (t/t) X m. (2.4) 
0 


In the above, cq is the decay parameter, and T1, T2 are the fast and slow time 
constants, respectively. 

Results of stress relaxation studies have revealed indifference in biomaterial 
responses to quasi-static, ramp-and-hold and oscillatory modes of loading along 
the structural axes of anisotropy. There is a shift of the stiffness and damping 
curves toward smaller frequencies of applied load and a decrease in the slope with 
higher stress levels indicating that larger stresses result in less relaxation and the 
damping is more effective at smaller frequencies. 

Biaxial tests to investigate in vitro pseudoelastic characteristics of the bladder 
wall tissue of different animals under quasi-static and dynamic loadings were 
conducted on square-shaped specimens. These studies allowed the deduction of 
full in-plane mechanical properties of the tissue. The edges of the specimens were 
aligned parallel and perpendicular to the orientation of the longitudinal and cir- 
cular smooth muscle fibers. The experimental protocol to obtain force-stretch ratio 
curves T? (4., Àj) used constant stretch ratios of 2): A. 


The in-plane passive T? , forces under biaxial loading are calculated as 


OpW 
P — 
ie ae ay UY (2.5) 


The most general form of the pseudo-strain energy density function W is 


pw =} lesa — 1? + 2c4(À4 T I) (4e = 1) + ca(Ac = i 
(2.6) 


+c6 exp(es( — 1)? -Ees( = 1)? + 2co(21 — = 1))], 


where p is the density of the undeformed tissue. 
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Bladders of pigs, rats, and dogs under biaxial loading exhibit a complex response 
including nonlinear pseudoelasticity, transverse anisotropy, and finite deformability, 
and no dependence on the stretch rate. The curves T? 1 (Acs 41) show that as the stretch 


ratio in one direction increases gradually, the extensibility along the other decreases 
(Fig. 2.3). There is a concomitant increase in the stiffness of the biomaterial. The 
maximum force the tissue can bear during the biaxial tests depends on the ratio 4;:/,. 
Experiments have shown that the shear force applied to the tissue is significantly less 


10? T? axe 1 Compared with the stretch force. Unfortunately, no experimental results 


have been obtained from the human urinary bladder for the analysis. 

Investigations into uniaxial and biaxial mechanical properties of actively con- 
tracting tissue remain the challenge in biomechanics. To the moment this book was 
written there were no experimental data available on in-plane active behavior of 
the wall of the urinary bladder. The main problem is to keep specimens physio- 
logically viable and stable, i.e., for in vitro samples to reproduce myoelectrical 
patterns that are consistent with those observed in vivo. Thus, it is practically 
impossible to sustain and control simultaneously spiking and contractile activity of 
smooth muscle syncytia. 

Over the years, surprisingly little attention has been given to the problem of 
constructing constitutive models of the bladder wall. The most common type of 
mathematical models used for soft tissues are phenomenological, and are usually 
of a polynomial or exponential form (Palmas and Rigato 1967; Kondo and Susset 
1973; Regnier et al. 1983). Thus, earlier proposed descriptions of viscoelastic 
properties of the tissue employed a combination of Maxwell and Hooke elements. 
Elastic, time and viscous parameters, and constants of the models were evaluated 
from experimentally recorded cystometry curves. Despite their robustness and 
practicality, they do not capture the underlying mechanisms of biomaterial 
behavior and, therefore, have failed to integrate information about tissue compo- 
sition and structure with its mechanical properties. 

The first detailed research on a multiphase structural constitutive model 
development for soft tissue with applications to the urinary bladder mechanics has 
been undertaken by Wognum (2010). Assuming that: 


Fig. 2.3 The strain energy 
function contour and fit of 
experimental biaxial data 
(black lines) for the rat 
bladder. Reproduced with 
permission from Wognum 
(Wognum 2010) 
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1. The soft tissue is an idealized network of muscle and undulated (in undeformed 
state) collagen fibers embedded into a compliant ground matrix; 

2. The mechanical net response is the sum of responses of individual fibers; and 

3. The tissue is incompressible; and 

4. The strain energy density function W satisfies the decomposition. 


W = dbecmWecm(E) + smWsm(E), (2.8) 


where gcm, dsm are the extracellular matrix and smooth muscle volume frac- 
tions, respectively, and E is the Green—Lagrange strain tensor, the Piola—Kirchhoff 
stress S is obtained as 


S(E) = ——ImC™! = $gcuSgcM (Esc) + @sySsm(Esm) ^ lnC |. — (2.9) 


Here lm is the Lagrange multiplier and C is the right Cauchy-Green strain 
tensor. The subsequent recruitment of fibers during loading in an ensemble of 
weight bearing elements suggests that a complete stress in the ECM and SM are 


n/2 Eos 
S(Ej) = nds J RÔ) J POEA 7 @7d0i = ECM, SM, (2.10) 
—n/2 0 


where the parameters and functions are referred to the collagen/smooth muscle fiber: 


n is the modulus, $, is the volumetric fraction, R( 0) is the distribution function, 0 is 
the angle orientation in the undeformed configuration, Eens is the fiber ensemble 
slack strain, D(x) is the recruitment function, r is the orientation vector. The exact 


forms of R(0) and D(x) are assumed to be known a priori. Thus, for collagen fibers 
the common representation for collagen fibers is bimodal beta distribution, and for 
smooth muscle—the two-parameter exponential model is widely adopted. 

The general nature of the approach has a potential to explain underlying 
remodeling mechanisms of individual constituents, i.e., the process of uncoiling, 
straightening, and reorientation of fibers along the direction of the applied force, 
under normal physiological conditions and to estimate their role in various 
pathologies. 


2.2 Modeling of the Bladder 


Publications on mathematical modeling of the biomechanics of the urinary bladder 
are scarce. The primary focus has been on modeling of the bladder as a reservoir 
(Regnier et al. 1983; van Beek 1997; Vlastelica et al. 2007; Spirka and Damaser 
2007; Korkmaz and Rogg 2007; Fry et al. 2011), its neuronal control mechanisms 
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(Hosein and Griffiths 1990; Kinder et al. 1999; Bastiaanssen et al. 1996a, b; van 
Duin et al. 2000; Fernandez et al. 2004; Pérez et al. 2008), and the simulation of 
urine-bladder-urethra interactions during filling and micturition (Damaser and 
Lehman 1993; Fletcher et al. 1997; Damaser 1999; Kren et al. 2001; Celik et al. 
2007; Zang et al. 2010; Kim et al. 2011). All existing models, owing to the number 
of uncertainties, suffer from considerable biological naiveté that makes them 
irrelevant in clinical applications. As a result, no comprehensible explanation of 
normal and pathological lower urinary tract behavior has been proposed so far 
based on the in silico approach. 

The urinary bladder has been typically treated as a thin shell subjected to 
external and internal quasi-static loads. Initially, it was suggested that the organ 
could be approximated by a geometrically simple shape—a sphere. Assuming 
incompressibility, homogeneity, mechanical isotropy, and physical linearity of the 
wall, Regnier et al. (1983) established the dependence of the inflation pressure on 
the extension ratios 4; (i = 1, 2, 3) 


A 
= —W(,4,4 ^jd4 2.11 
y= [4 uro ài, (2.11 
1 


where the following form of strain energy function W was 


W(A,4,4-7) = c(po, p) + poP(p, 4)e 9. (2.12) 


Here po, p are the initial intravesicular pressure and the rate of pressure change, 
respectively, P(p, 2) is a third degree polynomial, c(po, p) is a constant. The cal- 
culated Lagrangian stresses demonstrated a good correlation between the predicted 
and experimental stress-strain relations recorded from uniaxial loading tests of 
excised dog bladder strips, as well as slow and rapid cystometry studies performed 
on the whole organ. 

Damaser and Lehman (1993) investigated the role of more complex bladder 
shapes—the prolate and oblate spheroids—on the intravesicular pressure-volume 
P(V) response during filling. Results showed that the p(V) curves obtained for 
prolate and oblate spheroids for a wide range of eccentricity, 0.1 « £ « 0.9, did not 
deviate significantly from the sphere. Later, Damaser (1999) incorporated low 
muscle tone in the model and reconfirmed their earlier findings. 

Křen et al. (2001) studied the interaction of urine with the elastic urethra and 
bladder during filling and micturition. The bladder-urethra system was modeled as 
an ellipsoid connected to a cylindrical tube. The wall of the bladder and the 
urethra, were considered to be isotropic and linear elastic. The urine was a 
Newtonian fluid and the flow was nonstationary and turbulent. To reproduce 
dynamic (external body and volume related) forces, "artificial" loading terms 
were included in the governing system of equations. Although authors provide 
graphical outputs for the urine velocity outflow, pressure distribution and changes 
in the configuration of the urinary tract during voiding, the results of simulations 
are inconclusive, i.e., no detailed quantitative information could be derived from 
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the data presented about the dynamics of micturition. A similar quasi-empirical 
analysis of variation in the intravesicular pressure, bladder volume, internal 
sphincter area, and urine outflow was conducted by Celik et al. (2007). In the 
model parasympathetic signaling to the detrusor muscle and internal sphincter was 
achieved through an a priori given signal. It was chosen to be a unit square 
function for Vy «V < Veit, where Verit is the critical volume when micturition 
starts, and a sinusoidal oscillatory function to mimic intermittent voiding. Using 
experimental curves, e.g., cystograms, and adjusting a number of computational 
parameters during simulations, the authors achieved satisfactory resemblance 
quantitatively and qualitatively to clinical observations. 

An inverse problem to determine the stress-strain relationships of the bladder 
wall was considered by Korkmaz and Rogg (2007). Based on the assumption that 
the organ is a thin spherical shell subjected to complex loading, the wall tension, c, 
was calculated from 


r 


9 = 55 (P — Paa), (2.13) 


where Papq is the intra-abdominal pressure, r is radius, and 6 is the thickness of the 
sphere. Using actual cystometry and uroflowmetry readings as input data it was 
possible to reproduce the nonlinear profile of the c(4) curve (4—stretch ratio) 
either in the longitudinal or circumferential directions. 

van Beek (1997) modeled the bladder as a thick-walled sphere. A special 
emphasis was given to the evaluation of the effect of muscle fibre orientation and 
physical nonlinearity of the tissue on the dynamics of deformation of the organ. 
Three different geometries of fiber winding were studied: the circumferential, 
longitudinal, and oblique. The total effective Cauchy stress, o^, was decomposed 


o° = Gb t Oa, (2.14) 


where o5 is the passive stress is a result of stretch of the connective tissue ele- 
ments, c4 is the active stress generated by muscle cells. The passive principal 
stress was obtained as 

OW(E) 


= ae (2.15) 


where the strain energy function, W, was chosen in the form 
W(E) = bile + boll + clexp(ailg + aze) — 1]. 


Here Jg, Ig are the first and second invariants of E, and a), a2,b,,b2,c are 
empirical constants. 

The classical Hill model composed of a parallel passive elastic, c}, and a series 
of passive elastic and contractile elements was used to simulate the active stress. 
The activation, A,(t) deactivation, Aq(t), and the length dependence, Aj(/,), 
functions were introduced to mimic the effect of regulatory mechanisms. They 
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were derived from experimental curves on single smooth muscle cells of a pig 
bladder 


1 1 : 4? t < te 
i= tae =, A= 1 + [(te — #)/ta] , Q6) 
O,t>te 


Aill) = (le — (Lmax — bw)) (be (Imax + dw) 


Here ¢,, fa are the time constants of contraction and relaxation, respectively, J, 
is the actual length of a contractile element, and Jax is the length of a contractile 
element at maximal active force, lw is the rate of curve width. 

The cavity pressure during filling was calculated from 


45 
| f0WjoA 
|J 4-1 


A 


p dà, (2.17) 


where the limits of integration are referred to the stretch ratios of outer and inner 
radii of the sphere, and during micturition it changed according to urine flow Q, 
p= i? + Ro, (2.18) 
Ot 
here i and R are the inertia and urethral resistance. 

Results of numerical simulations in the case of circular winding of muscle 
fibers demonstrated a gradual change of the spherical configuration of the shell 
(bladder) to a prolate spheroid during the collection phase and its reversal to the 
initial state at the end of void. Analysis of stress distribution showed the highest o° 
values at the inner wall in the equator region and the lowest at the poles. The 
bladder with longitudinal fiber geometry attained the shape of an oblate spheroid at 
the end of filling with the pattern of stress-strain distribution inverse to the case of 
circumferential geometry. 

While this was a worthy effort by the author, the investigation suffers from the 
lack of an accurate mathematical formulation of the problem of the dynamics of 
thick shell. Moreover, the designed numerical algorithm proved to have poor con- 
vergence and stability. Therefore, the results per se are of a limited scientific value. 

A conceptually analogous approach was employed by Vlastelica et al. (2007) in 
their numerical modeling of the urinary bladder as an axisymmetric thin biolog- 
ically active shell. The initial configuration of the organ closely resembled the one 
observed clinically. The wall tissue was assumed to be homogeneous and to 
possess nonlinear elastic properties. However, the supposition is valid for “pas- 
sive" biomaterials only and becomes invalid when any muscle activity is present 
or suspected. Results of simulations of the passive relaxation and detrusor acti- 
vation phases showed that maximum circumferential stress occurred at the lower 
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part of the bladder and at the fundus. The body of the organ experienced a uniform 
axial and circumferential stress distribution throughout. 

Although the mathematical models described above have provided useful 
insights into the investigation of clinical disorders related to bladder hypertony and 
urethral obstruction, they lack valid anatomical and physiological inputs, i.e., the 
neural regulatory elements, internal/external bladder and urethral sphincters, urine 
flow dynamics and neurotransmission mechanisms, among others. The first plau- 
sible attempt to address these issues was undertaken by Bastiaanssen et al. (1996a, b) 
in their myocybernetic model of the lower urinary tract. Authors exercised a system 
biology approach to simulate the bladder as a multihierarchical system that inte- 
grated neural network, biomechanical, and urodynamic components in one quan- 
titative model. The dynamics of neural excitatory and inhibitory stimuli that trigger 
detrusor and urethral sphincter contraction and relaxation, respectively, are 
described in terms of normalized activity functions fip, fas 


z dfap 
dt 


= oP —fp(l—oP) fa € [0,1] (2.19) 


2 = @s—fas fas € [0,1], (2.20) 


where Tt is a time constant, w?, œP, c, are the excitatory (subscript e), inhibitory 
(i), and excitatory somatic (s) inputs ((@?, œP, @,) € [0, 1]). 

Assuming quasi-static isometric contractions the detrusor was treated as a 
homogeneous incompressible viscoelastic continuum. The total actual tensile 
stresses in the smooth muscle and sphincter was decomposed into a sum of the 
active and passive stresses 


Op = kara (fad Omax 0v (Vo) 01 (Ip) s (Ip) = Ove(VD)); (2.21) 


Os = Kick (fasOmaxFv(Vs)oi(ls) + Ge(Is)). (2.22) 


Here, all functions are normalized on [0,1], max, ei(/p), oi(/s) are the maximal 
and actual isometric stresses of the detrusor and urethral sphincter, respectively, 
Ovy(vp), ev(vs) are velocity-dependent stresses obtained at optimal length of 
muscle fibers that satisfy the Hill equation, o.(/p), c«(Is), @ve(vp) are the elastic 
and viscoelastic stresses, Karea, Kthick are the ratio of the cross-sectional area of the 
detrusor taken at optimal muscle fiber length to the actual area and thickness of the 
sphincter, respectively. 

The bladder was considered to be a thin spherical isotropic shell and the 
sphincter was a cylinder. The intravesicular and active sphincter pressures were 
calculated as: pps) = p(s) In(rpis /rp(s): where Ts): lbs) are the outer and 
internal radii of the shell and the sphincter, respectively. With regard to the passive 
elastic properties of the urethral sphincter, an exponential relationship between the 
radius of the urethra and urine pressure was selected. 
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Steady flow assumption was made to describe flow in different parts of the 
bladder-urethra system. The flow Q was obtained from 


Pu ; 
ES Z 2.23 
T (o tet (2.23) 


where Pu (Pu = ps) is the liquid pressure in the urethra, and R,(R;,A;) is the total 
resistance of the distal urethra, a function of the resistance and cross-sectional 
areas of the proximal and distal parts of the urethra. 

In the above model, most parameters have both physiological and physical 
meanings. Therefore, it was possible for authors to define a range of their variation 
and measure their actual values. Despite its mathematical and biological naiveté, 
the model results clearly resembled traces obtained from urodynamic studies. 
Several other models of the lower urinary tract have been proposed which are, to a 
certain extent, variations of the Bastiaanssen et al. model. However, all of these 
models lack solid biological and mechanical foundations by design. Therefore, 
they cannot answer important questions related to pathophysiological changes that 
occur in signal transduction mechanisms and biomechanical activity, or predict the 
effects of pharmacological interventions in various diseases. 

Finally, a new class of models have emerged recently, to study the dynamics of 
urinary flow in the bladder and urethra. They are a result of advancements in 
computational, software and imaging technology, MRI with digital three-dimen- 
sional image processing, and image reconstruction, in particular. The approach 
focuses primarily on detailed description of geometry and the application of solid 
and fluid mechanic principles to simulate the lower urinary tract anatomical 
structures and their functions. Although existing models are still biologically and 
biomechanically deficient, i.e., no active forces and deformations are generated in 
the bladder and the urethra during micturition, and no regulatory (parasympathetic 
and sympathetic) inputs are present, they offer original insight into the dynamics of 
bladder-urethra-urine flow. For example, the occurrence of a secondary flow at the 
base of the bladder, which is a result of the action of convective acceleration of 
fluid and viscous friction, has been demonstrated using a computational model by 
Zang et al. (2010). Computational models, when they reach the desired level of 
biological consistency, reliability, and robustness in their design and performance, 
will undeniably become an indispensable clinical tool. They will help solve 
puzzling intricacies of pathophysiological mechanisms of diseases, which are 
inaccessible to investigations in vivo, and will facilitate the development of 
rational pharmacotherapies that do not have deleterious side effects. 

We shall employ a combined reduction—integration systems biology approach 
as the main conceptual strategy to study the urinary bladder. The approach is being 
proven as the most effective and successful scientific method for an in-depth 
investigation of individual biological components and to understand the nonlinear 
multiparametric dynamic interrelationships among them. We shall start our anal- 
ysis with the essentials that are required to understand the concepts involved in 
development of a model of the urinary bladder. Any reader familiar with the 
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subject can proceed directly to Chap. 6. Then we shall study the biomechanics of 
the detrusor fasciculus, followed by the detrusor—intramural ganglion arrange- 
ment, and, finally, shall provide a holistic perspective of the urinary bladder as a 
thin biologically active self-controlled system. 
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Chapter 3 
Geometry of Thin Shells 


Order and simplification are the first steps toward the mastery 
of a subject; the actual enemy is the unknown. 
T. Mann. 


3.1 The Bladder as a Thin Biological Shell 


A solid body bounded by two closely spaced curved surfaces is called a shell. 
Consider a surface S within the shell. Any point M on S can be associated with 
general curvilinear coordinates «!, o? and the unit normal vector m. Let the 
perpendicular distance along m be given by z: —0.5h(a!, o2) < z < 0.5h(a', o2), 
where h is the thickness of the shell. The surface S for which z = 0 is called the 
middle surface. The two surfaces defined by z = +0.5h(a!, 02) are called the faces 
of the shell. Throughout this book we assume that the faces of the shell are smooth 
and contain no singularities, i.e., there are no structural defects or inclusions in the 
wall of the bladder and the thickness h is constant (Fig. 3.1). 

The shell is classified as thin if the max(h/R;) « 1/20, where R; are the radii of 
the curvature of S. However, the estimate is very rough and in various practical 
applications, other geometric and mechanical characteristics should also be 
considered. 

Based on experimental data and clinical observations that: 


1. the thickness of the bladder wall varies from A, = 3.3 + 1.1 mm; the char- 
acteristic radii of the curvature of the middle surface are: R, ~ 8 and 
R ~ 10 cm, respectively; thus max(A/R;) < 1/20 (i = 1,2); 

2. at full fill, the height of the organ measures on average 10 cm, and 
hiL = 4 x 102; 

3. as the bladder empties, the organ changes its shape and size to a "collapsed 
vessel"; 

4. stress-strain curves of uniaxial and biaxial stretching of the bladder tissue in 
in vitro experiments, have shown high degree of nonlinearity and large 
deformability, 


it becomes clear that the bladder can be represented as a thin soft biological 
shell (bioshell). 
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Fig. 3.1 A thin shell 


3.2 Geometry of the Shell 


Let the middle surface S of the bioshell be associated with a set of independent 
parameters «! and o? (Fig. 3.2). Putting «! = const and varying the parameter o? 
we obtain a curve that lies entirely on S. Successively giving o! a series of constant 
values we obtain a family of curves along which only a parameter o? varies. These 
curves are called the o?—coordinate lines. Similarly, setting o? = const we obtain 
the «!—coordinate lines of S. We assume that only one curve of the family passes 
through a point of the given surface. Thus, any point M on S can be treated as a 
cross-intersection of the gl- and «?-curvilinear coordinate lines. 

The position of a point M with respect to the origin O of the reference system is 
defined by the position vector F. Differentiating F with respect to « (i = 1, 2) 
vectors tangent to the «!- and «?-coordinate lines can be found 

naa, n= (3.1) 


Fig. 3.2 Parameterization of 
the surface of a shell 
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The vector m normal to r; and r» is found from 
in = |F, r] and mr, = 0, mr, = 0 (3.2) 


where [71,72] is the vector product. The vectors 71,72, and m are linearly inde- 
pendent and comprise a covariant base {7,,7,m} on S. We introduce the unit 
vectors e; in the direction of r; as 


2 r1 T, - T2 T? 
Ch ———. a C2 St wm? 3.3 
1 Ir. A Ir] A ( ) 
where 
[ri| 2 77 =A = yan, |2| = Toro = A» = yan, Tif. = A1A2 cos 4 = yar. 


(3.4) 


Here y is the angle between coordinate lines, A; are the Lamé parameters. 
Covariant components of the first A (a;,) and second B (b,x) metric tensors on 
S are defined by 


Qik = Tif, bie = —rimy = —rymi; = Tim, (i,k = 1,2). (3.5) 


Mixed components of the second metric tensor are b, = abs. Similarly, covariant 
(cik) and contravariant (c^) components of the discriminant tensor C are given by 


m|, Fi] = ca, |r, 7] = ci, (3.6) 
Hence, 
cum = [na] c*m- [r] (3.7) 
cur” = [m,r] c*ri = [mr] 
it follows that 
Ä=0, Gr Oye Seay cd =en = ya (3.8) 


" 
km __ sm ki o sil. 
Cic = 7. cac? = Ò; = 


Hereforward the summation convention is applied, i.e., when an index variable 
appears twice in a single term, once in an upper (superscript) and once in a lower 
(subscript) position, it implies that we are summing over all of its possible values. 

Vectors 7! are reciprocal to 7; (F = a^,) such that 


rne. rame. rm Sar. F7 = a", (3.9) 


here & is the Kronecker delta 
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Covariant derivatives of a vector or a tensor are found to be 


Via; = Ona; — Tio, Vial = Ox ct! + To, 


VjT* = T* +7" + TT Vila = OTi — TTY, — Ts, 


(3.10) 


where EA (17, = 7 Fix) are the Christoffel symbols of the second kind, V, is the 
sign of covariant differentiation with respect to aik, 0.0%; = Ook. 
The following Gauss-Weingarten equations hold 


TA ET E 
Ty = D} Tj + mb, Mi = 0m = —bir, 


or 
Vin = mby, Vi = mb, Vim = —bkr. (3.11) 


Note that the above formulas are also valid for the deformed surface s. 
(Henceforth, all quantities that refer to the deformed configuration we shall des- 
ignate by an asterisk (*) unless otherwise specified). 

Let u; and u’ be co- and contravariant components of the vector of displacement 


v(a , à?) referred to the deformed coordinate system on 8. Thus 
P—ro (xx), D= uF + mw = ur + ñw, (3.12) 


* * ; 
where F is the position vector of a point M on S, uj, u' are the tangent, and w is the 
normal displacements (deflection), respectively. Covariant differentiation of (3.12) 


" * 
with respect to «' for coordinate vectors r;, we find 


7; = Fel ð + e) + mo, (3.13) 
ej = Vu — whi, oi = Viw + bi, 


The covariant components ¢, of the tensor of deformations E on S are defined 
by 
* LS NS ; 
2£i, = Aik —Gik = Fi Fk —TiTk = eik + ey + a" eijeys + iOR,. (3.14) 
€ik = V ilk = Whig, (i,k = 1,2) 


* 

Here, use is made of Eq. (3.13) for r;. The components £11, £»? describe tangent 
deformations along the coordinate lines a! , x? and £5 is the shear deformation that 
characterizes the change in the angle between them. 


The vector m (m. 5) is found from 


(r'S; + mS3) 
Va, 


d 
m = 


(3.15) 
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where the following notations are introduced: 


S; = ek og — (1 tec e2)ox, S3 = (1 + el)(1 + e3) — eset, 


A = rdc (ce — ee) le) 


The quantities ej, Wi, S; S3 describe the rotations of tangent and normal 
vectors during deformation. 
Covariant components of the tensor of bending deformations on S defined by 


* * 


Seg = Dix — b, (by = —mjr,, bi = — mi f), (3.17) 


where æ; describe changes in curvatures, and «5 is the twist of the surface. 
Substituting E and Mr = 0; in from Eqs. (3.13) and (3.15) into the above, we find 


bik = -J (ex ViS3 + (i + e) ViSj — bi (S3(6} + ei) — «S!)). (3.18) 


Making use of the right-hand sides of equalities 
oxrViS3 + (01 + ef) VS; = —S3Viox — S;V ie], 
(07 + e?) VS; — ojS" = aj ((2 t ei t e)a" — (à; + ej)a"), 
Eq. (3.18) takes the final form 
1 


br = 
k V 


[Ss Vio + Sj Vie] + by ais((2 +e} + ea^ — (ò; ea) |. (3.19) 


3.3 Tensor of Affine Deformation 


The Christoffel symbols of the first kind of the undeformed and deformed middle 
surfaces are defined by 


E. = : (s i +ô; p -6; 2 (i,j,k = 1,2). (3.20) 
Substituting dik = dix + 26% from Eq. (3.14) into (3.20), we get 

La = Ij 4-Ox6jj + Oiej — Ojeis. (3.21) 
The right-hand side of (3.21) can be written in the form 


Oxeij + Õi£jk — Õisik = Pk + ZV nj, 
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where 
Pjik = Vii + Vk&ij — Visi. (3.22) 
and 
Vj£iy = ÕjEik = Mm 22 Din; Vy = aT jit, 


Expression (3.21) can be written in the form 


* 


Tjik = 2 t Pj + 2I nj 


which upon contracting of both sides with the tensor a yields 


* J . ; 
Dy =T} +A] (3.23) 


where AL = 4 Pak Note that although the Christoffel symbols are not tensors 
their difference, AJ, is a tensor of the third valency such that A! _— AL Its tensor is 
called the tensor of affine deformation (Niorden 1954). 

The tensor A, allows to establish a relationship between covariant derivatives 


» * 
of a vector A* with respect to aj, and aik 


V;A* = ViA ASA), | ViA, = ViAy — AĴA; (3.24) 


Similar relationships can be set for tensors of any valency. 
* 
Consider coordinate vectors 7, as tensors of first valency. Their covariant 


derivatives are given by 


* 


here V;—the sign of covariant differentiation in the metric dik . With the help of 
(3.11) we get 


+ ELE. "a 
Tg — Ug T+ big m, mj = —b; rr (x) 
The following is true: 
c Se js 
Vite = bym, | Vim — — bj re 
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from where it follows that 
Vif, = A5 Fe +b, Vim Vim bm. (3.25) 


ak 
Here, use is made of Eq. (3.24). For vectors of the reciprocal basis, r;, we have 


Vur-r-ptsb, yrs. kn 
xk "E^ * oa 
Substituting 7, into the first expression, we find: V; 7^ = bt m. Making use of 
Eq. (3.23), the final formulas are found to be 
x * + " " * xk 
Vit = —Air +bi im, im; = bar, (3.26) 
Multiplying Eqs. (3.25) and (3.26) by 5. we get 
Ab =F Vif, A=- RV. (3.27) 


Contracting both sides of the fist equality with nk,» we have 


k F g 


Paik = Tn ViTk. (3.28) 


To express the tensor P, in terms of displacements we proceed from differ- 
entiating of Eq. (3.13) and making use of (*) for the undeformed surface S. As a 
result we obtain 


" (1) Q) 
Vir, = (0j + ej) Vir, + r, Vie; + m;oy + mV;oy = Wi new, +e) 


(**) 


where 


(1) (2) 
Wi = Vie, — bio, Wy = Vier + biei. 


Substituting the resultant equation into (3.28), we have 


(2) (1) 
Pk — On (wa eu) + Wi (Ans + ens). (3.29) 


On introducing the right-hand side of (**) into the obvious equality c*V; 5 — 0, 
we obtain 
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" (1) (2) 
c* | W3 7, + m| W, +bx ] | =0. 
It is clear that the tensors W}, Hig are symmetrical 
"mU "mo 
c*Wi-0, c“W,=0. (3.30) 

or in expanded form 

c*Viei = c*bio; c'"Vio, = —c'bei, 

kV er, = c"byo;, VR = —c* pie. 


The last two formulas are obtained by contracting the first two with ajs. 
The scalar multiplication of the first equation (3.25) by m, yields bj; — m Vi 5 : 
On substituting m from Eq. (3.15) and Vj 5 from Eq. (**), for the coefficients big 


of the deformed surface 8, we obtain 


* 1 (2) (1) 
j= ag (5s wi |. (3.31) 


3.4 Equations of Continuity of Deformations 


For the surface to retain continuity during deformation, the parameters ¢, and cj; 
(i, k = 1, 2) must satisfy the equations of continuity of deformations. To derive 
them we proceed from covariant differentiation of Eq. (3.25) with respect to o" 


* * * 


* * * * 
ViVi r=m Vn bik Tm, bik +P 5 ikV n pr + r VaPs.ik- 


Substituting expressions for V, r* and m given by Eq. (3.13), we find 


* * * *o ok Xx * * * * 
P = =s s = s = =s 
ViVi r =m Vn bik x bi Das r —AyPsik r! + b, Ps ix m+r VaP; ik 


* * 
Multiplying the above by 7; and m, we obtain two scalar equations 


5 VaV; 5 — big bj AP s jn VaPji (3.32) 


mV, V, 5 =V; by + bi Pj ik- (3:33) 


On substituting V,,P; ix from Eq. (3.32) into the equality 
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Vp (Pir + Pig) = 2VnVkêij 
and contracting the resultant equation with the tensor c'"c/*, we get 
in jk mu S = = 
cin (-2V,Ve8y + bikbyj + ARP sj +T; VanVi n). (3.34) 
In deriving (3.34) use is made of equalities 
ik * * * ik * 
c" ri ViVi Tk = Ti Vac! Vi Ty = 0, 


that follow from the requirements of integrability of c^ V j 5 = 0. The last term in 
Eq. (3.34) can be recast using the Ricci equality given by 


2c" V,V; 5 = eR fs (3.35) 


where R$; are the mixed components of the Riemann-Christoffel tensor of the 
surface defined by 


R? g = barb; — bib. (3.36) 
Thus, we have 
m c" V,Ni T, — 0, (* * x) 


"P * ]1.. " l.. "» 
oe V n = 5€ € ajs = zc” (bubi — Digbnj) + "PR g e. (3*) 


Substituting the right-hand side of Eq. (3.33) into (***), the Gauss-Codazzi 
equations for the deformed surface are found to be 


c” (v. bi, + b] Pin) =0. (3.37) 


Analogously, substituting Eq. (3.34) into (3*), we find a relationship between 
the coefficients of the first and second fundamental forms of the deformed surface 


cic (20, Vat + baby — baby + APs + Rix Gs) =0. (3.38) 
Here, use is made of the obvious fact that 
c" cp ibi = —c"c* bib, 
Applying the Bianchi formular to the last term of (3.38), we find 
Ria = KeniCks, (3.39) 
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where K is the Gaussian curvature of S defined by 
K = bib? — bhp. (3.40) 
Since 


in _ i jk s = Ak S— pam j 
c"cy = —Ój— 2, cce = eese; = —Ópej = —8j, 


then for R* 


ni 


4^ 6j; We have 
c" PRS g Ejs = Kel" F Cni ejs = Kòl cie = 2Kel = 2K(el - 2). (3.41) 


Introducing (3.41) into (3.38) and remembering (4.17) one of the equations of 
continuity of finite deformations is found to be 


xL 1 1 , f ; 
ete (Vave = 2 ike nj = i) = (2Ha'* = b*) ze; = Ka" ej = 0, 


(3.42) 
where H is the mean curvature of S and the following formulas hold: 
c" cp, = 2Ha? — bË, 2H = a* bi, = bi. 


With the help of the Codazzi conditions for the undeformed surface: c"V,bg = 


* 
0, and equalities: A7, = b/ P;ix, two additional equations of compatibility of 


deformations can be obtained from Eqs. (3.37) in the form 


c" V gei = c" (b, — Rin AY (3.43) 
Contracting both sides of (3.43) with the tensor c", we get 


II goa, = cF (Dns — ge, AT. (3.44) 


3.5 Equations of Equilibrium 


Although a thin shell can be treated as a three-dimensional solid the complexity of 
the problem would be reduced significantly if its dimensionality could be reduced 
from three to two. Applying the second Kirchhoff-Love hypothesis that "the 
transverse normal stress is significantly smaller compared to other stresses in the 
shell and thus may be neglected", and recalling that the deformed state of the shell 
is completely defined in terms of deformations and curvatures of its middle sur- 
face, the shell can be regarded as a two-dimensional solid. Furthermore, restricting 
our considerations to the membrane theory of shells, i.e., the lateral forces and 
moments acting upon the element are assumed to equal zero, the stress state of a 
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differential element of the membrane is described entirely by in-plane tangent 
T;(T") and shear T4(T^)(i Æ k) forces per unit length of the element. 

To derive the equations of equilibrium of a shell in terms of the deformed 
configuration, we proceed from the vector equations of equilibrium for a three- 
dimensional solid given by (Galimov 1975) 


dva, Ed 3 yag 0, (3.45) 


On! 
where 
* * 
zi ik iz Š 7 
f =o" rR+o%m (i,k= 1,2). 
fum 
Here o*, o^ are the contravariant components of the stress tensor and the vector 

D^ of external forces acting upon the faces z= +0.5h of the deformed shell, 


* 


respectively, F is the vector of mass forces per unit volume of the deformed 
element, à — dei. i x). Integrating Eq. (3.45) and subsequently the vector product 


of (3.45) by mz over the thickness of the shell, z € [z1, z2], we obtain 


"ik * * 
oT Va +p a 4 ferat-o (3.46) 
E al 


[re ; ei (3.47) 


By introducing Ta and m; given by Eq. (*) into (3.46) and equating to zero the 


coefficients for A and m, we find 
V, ENT 3 (4 ES! 
5. ik rik eie 


P "LT à —Zby)dz = 0 


a (3.48) 


8 73 qe — 
es yira =o. 
Zl 


Here, 


It follows immediately from Eq. (3.47) that the tensor T^ is symmetrical: 
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TU TA. 

Although Eqs. (3.46)-(3.48) are derived under the assumption h = constant, 
they are also valid for shells of variable thickness h = A (a, o). 

Let the components o and in-plane forces 7" be defined by 


. r * : r * . 22 zl es "EE " 
gÉ = g* V 8/8, rik — pik V a/a, TK = J "T (a — isa 
£l 


Hence, the equations of equilibrium (3.48) in terms of the undeformed con- 


figuration take the form 
[rr yari ôi Z bix)dž =0 


V;T* + AET^4 Vert Ea 
a 
bs T* 4+ + jel 33 Eu or dz—0 (3.49) 
2g 


where AL and V;T^ are calculated, with the use of (3.23), as 


Png = Vien + Vieni — Vnêij 


ik _ 1 OyaT^ kij 
VT" Se o FI T; 
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Chapter 4 
Essentials of the Theory of Soft Shells 


Scientists study the world as it is; engineers create the world 
that has never been. 
T. von Karman. 


4.1 Deformation of the Shell 


A class of thin shells that: 


hlL ~ 107° = 1072, L is the characteristic dimension of a shell; 

possess zero flexural rigidity; 

do not withstand compression forces; 

their actual configuration depends entirely on internal/external loads distributed 
over unit surface area; 

5. undergo finite deformations; 

6. their stress-strain states are fully described by in-plane membrane forces per 
unit length; 


TCI 


is called soft shells. Let the middle surface of the undeformed soft shell S be 
parameterized by curvilinear coordinates «!, o? and VM(a!, à?) € S is described 
by a position vector 7(«!, o2). As a result of the action of external and or internal 


* 
loads the shell will deform to attain a new configuration S. We assume that 


"EL * 
deformation is such that VM («!, o?) — M(a! , o?) and is a homeomorphism. Thus, 
the inverse transformation exists. The transformation is defined analytically by 


zs ad (ad, x, (4.1) 


* 


" 
where (a! , 62) and (a, a?) are the coordinates on S and S, respectively. Assuming 


that Eq. (4.1) is continuously differentiable and det(0 o! /O0o^) 40 (i, k = 1, 2), 


then the inverse transformation exists: oj; = o; a! = oi (a, a). 
Let 
* 
oy 
k k 
C; = Da’ C= det(C;), (4.2) 
1 
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* 
cto 9 
io * 7 
Oj 


Assuming that C £ 0, C* # 0, we have 


€ = dei(c*) (4.3) 


GO 100-0 (x 7), 


CC 21, c - [ci]/c. (4.4) 


Here [Cİ] is the cofactor to the element C; of the matrix (CH). 


2 


Deformation of linear elements along the &!, o?—coordinate lines is described 


by stretch ratios 4; (i = 1, 2) and elongations e,; 


ds; — ds; di 
= ay ee 1, 4.5 
É ds; dii a 


* 
where ds;, dsj—lengths of a line element between two infinitely close points on 
* 
S and S, respectively, and are given by 
è * * *. 
ds; = vaida’, dsi = aid o. 


Changes in the angle,y, between coordinate lines and the surface area,ds,q, are 
described by 


(0) * (0 -1 412 
"m 4.6 
Y= %X-L= X-cO anas (4.6) 
ds VE à a sing "E 
ju en ease a9 je (4.7) 


= = 0 
dsi vya aaa sin x 


Making use of Eqs. (4.2), (4.3), for vectors 7;, 5 tangent to coordinate lines on 


S and H we have 
5n- Cn, Ti = R. (4.8) 


* * 
Hence, the unit vectors e; € S, e; € S are found to be 


li Ti k| [di ^k 
= Ci |r| =C 


==) — k 
I| vai dk 
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= Ap. 
e; = Cie, 


55 


(4.9) 


where the following notations are introduced 


(4.10) 


With the help of Eq. (4.10), the scalar and vector products of unit vectors e and 


* 
€, are found to be 


xo 


SE pa * 
6j: Ek :— COS Hj, = 


So PNE. 
DES = MSIN Xik = 


In just the same way, proceeding 
by ex, it can be shown that 

(0) 

COS Yik 


. (0) 
sin Ji 


* 


Cl 


"M Ai Ai 5 (0) 
ej: €; C; Cg = Cj Ck COS Fins (4.11) 
H E Aj An Aj Ea * : (0) 
ej X e, C; Ck = Cj Cy sin yj. 


from the scalar and vector multiplication of e; 


(QU COS Yin; (4.12) 
CC sin y, . 


* 
To calculate the coefficients ed c we proceed from geometric considerations. 


* 


" 
Let vectors 2;, @ at point w(s. 


* 


3 € S be oriented as shown (Fig. 4.1). 


* 
Decomposing e; in the directions of ej, we have 


Ci = MC, C? = CD, Ci = —AB, C2 = MB. (4.13) 
Solving AMCD and AMBA, we find 
M sin(y—y;) + sin 7. 
C = : E) & = : xi 
sin y sin y 
^1 sin(z; - Xo — x) £a sin(a + x2) 
C = > , cy = —_—, (4.14) 
sin y sin y 
Ê = auct = 9 
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Fig. 4.1 Deformation of an element of the soft shell 


* 
Similarly, expanding unit vectors e; along ej, we obtain 


0 0 0 
^i Sin(i + x2) (aom 

1 0 Voy) 0? 
sin y, sin y, 


4 0 0 0 . ,0 0 
^ sin(%ı + X2 — X) a _ sin(y — %) 


C= — Q6 MV, (4.15) 
sin y, sin y, 
sin? 
k 
C — detC; — > 
sin yi 


* * 
Note, that the coefficients e. ck are the functions of y, and Ls while e c 
depend on 7; and 5 and the actual configuration of a shell. 
0 
Let the cut configuration of a soft shell S be different from the undeformed 


0 * 
configuration S. We introduce the coefficients of transformation F; € $ — S by 


ê: = Cty aix / di, e =Ç V ay. / di, (4.16) 


0 : 0 * ; 
where dii, dj are the components of the metric tensor A on § and S, respectively. 


Eliminating e ck from Eq. (4.10), for the coefficients of the cut and deformed 
surfaces we obtain 
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^ 2 A x = 2 
C=O >. =C =. (4.17) 
Ài n 
Similarly, we introduce the coefficients 
ees f6116.8 e 1818, (4.18) 
where 
(4.19) 
(4.20) 
0 
Let E be the tensor of deformation of S(S — S)given by 
E= err, (4.21) 
where 
a — dix 
€ik = : 2 . (4.22) 
Substituting Eq. (4.2) into (4.22) for £i, we find 
(Ait cos 7 — cos x) / dii kk 
€jk = s (4.23) 
2 
It is easy to show that the following relations hold 
eik ccu qd = d cick, 
5 (4.24) 


ex 


* 
= Jon ik || Qjnpuipk 
€ik — Ejn C; Ch, es = g C;C,. 


In the theory of soft thin shells stretch ratios and membrane forces per unit 
length of a differential element are preferred to traditional deformations and 
stresses per unit cross-sectional area of the shell. Thus, dividing (4.23) by the 
surface area ,/aiayy of an element, we get 


Eik (4 COS 7 — cos x) 
čik := : = ; (4.25) 
i Ckk 2 
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& are called the physical components of E. Using Eq. (4.24), for či in terms of the 
deformed configuration, we obtain 


* * x aa 

~ _ z j ^n Viji 

čik = & C; C; , (4.26) 
* * 
Aii Ann 


where the coefficients Ci satisfy Eqs. (4.2), (4.3). 
Making use of Eq. (4.25) in (4.3) and (4.6) for 4; and y, we find 


Ai = 1 + = V1 + 268i, 


~ (0) 
0 2 
y= b cos! = T cos X —. (4.27) 
va + 2580) /(1 + 2852) 


Substituting C^, C^ given by Eq. (4.16) in (4.26), we have 


ae ke 
* 2; 


B. = CC, Bx CC. (4.28) 


* 0 
Finally, formulas for gix in terms of S—configuration of the soft shell takes the 
form 


= Tf 205 - 0, oe. a OS ES ee Oe 0 
ču = m [n sin?(y — %2) + € sin? y; + Wy sin(y — %2) sin y, , 
22 
sin” y 
>o | p. .,,0 0,.,0 0 0... . 0. 0° 0 
£12 SEE |-u sin(y — y5) sin(x1 + X2 — X) + € sin y» sin(y; + 72) 
sin^ y 


- 0 0 0 0 0 
Bi (cost -2y,—3)-— cos cos 71) |; 


e yo 0 0 0 zo -2 0 0 
i 1 [|Énusin^(y, + x5 — x) + én sin (y; + x2) (4.29) 
2 77750 ~. 0 0 o. 0 0 i 
sin? y — Zijo sin(y, + % — y)sin(y, + %2) 


With the help of Eq. (4.25) in (4.28) the physical components can also be 
expressed in terms of stretch ratios and shear angles as 


* o 
* 


d * 0 0. Aj Aan 
Ài Ak COS yg, — COS Hig = (AjAn COS Zin — COS Xin) C; Cy - (4.30) 


0 
Further, on use of Eqs. (4.11), (4.14), (4.30) in terms of S—configuration takes 
the form 


4.1 Deformation of the Shell 59 


T m * 7 ii An 0 

Aj Àk COS Xa, = AjÀn C; C, COS js; (4.31) 
or in expanded form 
* 2.3,0 0 2.50 TED 2,0. 0, , 0112 
Ay = T [a sin'(y — %2) + Ag sin” y +2142 sin(y — y) sin(% — 5) sin y; 


* = . ,0 0 ; 0 0 0 
y= xy; —cos ! I(-4 sin(y — x;)sin(y, + x2 — + 


.0.,0 0 4 0 0 
EE P sin %2 sin(y; + %2) + ArA2(cos(d + 2 y; — x)— 
žo 
—cos y cos 4i) cos(y -)) (Ai 4a sin? | 
7 1 .5,0 0 0 E 
A2 pie A sin? (y, tx—3)- 4: sin? (X; +%)- 
sin’ y (4.32) 


. 0 0 0 0, .,0 | 0472 
—2ÀJacos( —A y) sin(y, + %2 — X) sin(y, + £)| 


Formulars (4.31), (4.32) are preferred in practical applications particularly 
when dealing with finite deformations of shells. 


4.2 Principal Deformations 


* * 
At any point M € S, there exist two mutually orthogonal directions that remain 
orthogonal during deformation and along which the components of E, attain the 
maximum and minimum value. They are called the principal directions. 


0) «x 
To find the orientation of the principal axes, we proceed as follows. Let p 0 


(0 x * 
be the angles of the direction away from the base vectors e; € S$,ej € S, 


respectively. We assume that the cut and undeformed configurations are indis- 


0 
tinguishable S — S. Then, setting Lo = p in the first equation of (4.29), we have 


~~ 1.90 ~ .5,0 0 et hf OS: ee cU retos HO 
é sin’ y = £ sin?(y — p) + 2ép sin(y — p) sing + & sin? p (4.33) 


After simple rearrangements it can be written in the form 
0 7 0 
€ = do + bo cos2 ọ + co sin2 ọ, (4.34) 


where 
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1 l zo " 0 
ao = 2 ju + 8553) — Ecos y|, 
sin? y 
1 L. p - P" 0 0 
bo — Ü 5 nu &52) + (às i cos 7) / cos Z , (4.35) 
sin? y 
E wu 2 0 
€ = — n = y cos 7). 
sin? y 


Differentiating Eq. (4.34) with respect to p and equating the resultant equation 
(0) 
to zero, for the principal axes on the surface S, we find (bọ # 0) 
BS 0 7 . 0 
co 2(čıı cos 7 — Fi) sin z 


tan 2 @ == 7 T (4.36) 
0 £j cos2y — 2é12 cos y + Ex 


Substituting Eq. (4.36) into (4.34), we obtain the principal physical components 
£j,£? of E 


aa uh pO 2 
£12 = ag £ 4j bo +c 


fa m bs 0 
E (či + 8553) — 2j» cos x 


. 70 
2 sin? y 


" "ET. 
1 0 (či — 222) 0 0 
2 511 22 ~~ ae jee T 
+ -sint y — + 8, + E1182 cos? y — £y (£1 + €22) cos Y 


(4.37) 
Henceforth, we assume that maxe, is achieved in the direction of the principal 
axis defined by the angle 0; E p, and min &2—along the axis, defined by the 


angle pə = $ +n/2. Since for the principal directions a = 1/2, from the 
second Eq. (4.29), we find 


Zi = —bosin2 p -F cocos2 Q. (4.38) 


Dividing both sides of Eq. (4.39) by co cos2 9 and using Eq. (4.37), we find 
€12 = 0. Thus, there exist indeed two mutually orthogonal directions at VM (od )€ 


(0) 
S that remain orthogonal throughout deformation. Similar result, i.e. y = 0, can 


be obtained from Eq. (4.13) by setting % = 7/2. 
Substituting Eq. (4.25) in (4.36), (4.37) for the orientation of the principal axes 


0 
on § and the principal stretch ratios, we obtain 
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2 Arde cos(% — y» - 27 cos z] siny 


tan2 9, = = 0 a 3 Q^ (4.39) 
- B+ 2/ Arde cos(y — y) — Aj; cos z] COS Y 
0 0 
Pı = Pı * n/2, 
l fraia 
Aia = [à | 43) [2 A122 cos(% — Y) cos y + 
sin y 


+ (Ci + 12) JA + 4222 cos(2 % — y) cos y — (4.40) 
. 1/27 1/2 
—AA (Ai +4 2) cos(7 — y) cos 7) | 


To find the orientation of the principal axes on the deformed surface 8, consider 
a triangular element on S bounded by the two principal axes and the o, —coor- 
dinate line (Fig. 4.2). Geometric analysis leads to the following obvious equalities 


cos Q, = um ^i ^ cos sing, = E = An 
! ds / A i | ds A (4.41) 
* d : A 0 
iani = = = tanp 
d s; 1 


Fig. 4.2 Principal 
deformations 
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On use of Eqs. (4.39), (4.40) from the above, we find the angles for the principal 


axes Q, and @ = Qi c n/2. 
Finally, substituting Eqs. (4.7) and (4.25) into expressions for the first and 
second invariants of the tensor of deformation E defined by 


* * 
n =e + 2 = č + - 


= T 0 ~~ 
aeu 2612 COS Y + £2 |, 


MEN (4.42) 
IP = ee = ui» - (Bn) = — [68s - (53). 
sin? y 
for the principal stretch ratios and the shear angle, we get 
At + AP =(A1)? + (22)? 
1 -— 0 (4.43) 
a (4 + 4$ — 241% cos(x -1)cos7). 
sin^ y 


E 
Ai Ag sin(x — y) 


0 


AA = yı + uP + 4P = dA. cosy = 5 
sin y 


The last equation is also used to calculate the change of the surface area of S. 


4.3 Membrane Forces 


The stress state of a differential element of the soft shell is described entirely by in- 
plane tangent T;(7?) and shear T&(T*)(i Æ k) forces per unit length of the 
element. To study the equilibrium of the shell, we proceed from consideration of 
triangular elements AMAB and AMCD on S (Fig. 4.3). Analysis of force distri- 
bution in the elements yields 


-MA T B: + MBT "z; + ABT” z; = 0, 
(4.44) 


MD T^ e, 4- CDT"z, — MCT”z; = 0, 


where Cl = MC, C2 = CD, Cl = —AB, C2 = MB (Eq. 4.13). The scalar 
product of Eq. (4.44) and æ yields 


* * * 
lk lip 22 2i Ak pl 
T" = T'C; Cj - T"C; C5, 


* * * 
T* = TOO LT ĉl. 
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Fig. 4.5 Membrane forces in 
the soft shell 


where use is made of Eq. (4.9). Substituting e given by Eq. (4.18) for Cc. we find 
‘i 1 jn Ai À 
pt c? Ce, (4.45) 


On use of Eq. (4.15), the components of the membrane forces are found to be 


B 1 . 0 0 i 0 0 0 
T= 0 D sin? (%1 + y5) + T? sin? (y, + %2 — )*- 
sin y sin y, (4.46) 
2,0. 0 o. 0 0 
oq sin(y, + %2 — x) sin(y, + £)] 
: 1 "A" 0 mad o 0. Ui sd 0 
p =— 7 {-r" sin y; sin(y, + %2) + T? sin(y, + x; — y) sin(y — %) 
sin y sin y, 
0 0 0 0 0 
gru [eos +2 %2 — Y) — cos x cos A ) 
7 1 
TZ = 5 T! sin? p e sin(7 = 1) sin a +7” sin?(% = e) 


sin y sin y, 


Introducing the tensor of membrane forces T 


; l 
T = TYT, = — T'es, (4.47) 
sin y 
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where T" are the physical components of T, and using Eq. (4.9), T'*can be 


expressed in terms of T'*as 


* 


* * 
T = — Teje, = T" sip Ci (dus 


sin y sin y 


* 


sin yi 
Further, making use of Eq. (4.15), we obtain 


* 


pot xx 
Tes TEC (4.48) 
C 


Substituting ĉi given by Eq. (4.14), we get 


1 


E me ; * * ni : * * * 
r= : [F sin? i) +7 sin’(¥, + y — x) 


sin " sin y 
Es . * * * s pe * 
—2 T" sin(y, + x, — y) sin(y — zn 


~ a "A HET: * a : * * NR n * * 
pe e sin 7; sin(y — %2) + T? sin(y, + y; — y) sin(y; + %)+ 


+2 [cos +2 7% —x) — cos X cos A n sin y sin y, 


~ 1 3t * Hai * * * ~ * * 
Oe tu sin 7; 2 T" sin(y, + %2) sin z; +7" sin (y, + 2n 
sin y sin y, 
(4.49) 
Using Eq. (4.17) after simple rearrangements, Eq. (4.48) takes the form 
A142 zik 1 dis Sy Aj kk 
T" — IPC C. 4.50 
n Re ok o j "n ( ) 


Formulas (4.50) are preferred to (4.48) in applications. First, the coefficients 


A kk y ; ? 
C, are used in calculations of both deformations and membrane forces. Second, 


e depend only on parameterization of the initial configuration of the shell. 
Therefore, once calculated they can be used throughout. 
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4.4 Principal Membrane Forces 


* * 

As in the case of principal deformations at any point M € S, there exist two 

mutually orthogonal directions that remain orthogonal throughout deformation and 

along which T attains the extreme values. They are called the principal directions 
and the principal membrane forces, respectively. 

* * 

Assuming that the coordinates di € S and o; € S are related by the angle y, 


then, setting A = 1/2 and Yo = win Eq. (4.46), we find 


i 1 25 ey Bsa) ‘4 
p = T eos V +T” cos? (y — y) -- 2T? cos y cos(y — y) 


sin y 
x 1 * * " * " * " * 
pu (^r cos y sin y +T” sin(y —2 y) +T” cos(y — y) sin(y — in, 
sin y 
* 1 * " * * " * 
qe Te ir sin? y — 2T" sin(y — y) siny +T” sin?(y — in}. 
sin y 
(4.51) 
Equation (4.51) can be written in the form 
x ll * * 
T =a,+b,cos2w+c, sin2y 
x 12 * * 
T -—-b,sin2y-c;cos2y (4.52) 


x22 * * 
T =a,—b,cos2w—c, sin2y 


where the following notations are introduced 


Pae 1 l (rl + T”) — Teos}, 


1 1 * * 
E 5 (T! + T2) 4 (r: 4 T? cos j) esi]. 


ce; =T? + T”? cos 7. 


um * 
Differentiating T" with respect to y and equating the result to zero, we obtain 


12 22 M eid 
erc 2(T +T cos %) sin z 


E i " (4.53) 
bi T 42T! cos ; - T? cos2 y 


tan2y = 


Solving the above for y for the directional angles of the principal axes, we get 
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2 (a + T? cos i) sin 7 


tan2y = T zx 
1 T! + 2T! cosy + T? cos2 x (4.54) 


V, = t n/2. 
Substituting Eq. (4.53) into (4.54), we have 


T” = a + y b7 + e. T? = 0, T? = a — V5 + e 


From Eq. (4.52) the principal membrane forces Tı, T; are found to be 


1 (au +7 x 
T2 = Ê Ege at i ) iori gos 


sin Y 


j: "UI — T2y 4 (T12  TP(TU + T2) cosy + T! T? cos? i}. 
(4.55) 


Thus, at each point of the surface of the soft shell there are two mutually 
orthogonal directions that remain orthogonal throughout deformation. Henceforth, 
we assume that T; > 75, i.e. the maximum stress is in the direction of the principal 


axis defined by the angle y, and the minimum—by the angle y3. 
Analogously to the invariants of the tensor of deformation described by Eq. 
(4.43), we introduce the first and second invariants of T 


a . 1 * 
im =F +D = T"! +T? = —(r" 4 7 4 27" cos 7), 
sin y (4.56) 


it = TT, = TT]? pp — TTZ _ (rey, 


4.5 Equations of Motion in General Curvilinear 
Coordinates 


Let A D and A n be the surface area and mass of a differential element of the soft 
shell in undeformed and deformed configurations. Position of a point M € S at any 
moment of time ¢ is given by vector 7(a!, «, t). Densities of the material in 
undeformed, p, and deformed, p. states are defined by 

Am dm * . Am dm 


ok do Ata do Pg 


4.5 Equations of Motion in General Curvilinear Coordinates 67 
where 


do = vVada!do?, do = daar 


Applying the law of conservation of the mass to Eq. (4.57) we find 


dm = dm = pda Pa tata gal cada 


It follows that 
* a 
a 


Let p;(o, %2, t) be the resultant of the external, p(.,) (o1, oo, t), and internal, 
Poo, 0, t), forces distributed over the outer and inner surfaces of the shell 


ps(01, %2, t) = Pip (%1, %2, f) + pL) (01, %2, t). 


The density of the resultant force per unit area of a deformed element p, is 
defined by 


Ds 


D t) = li x 4.59 
P(%1, oo, t) Ui E (4.59) 
Similarly, we introduce the density of the mass force F(a!, «7, t) by 
" F dF 1dF 
f(a', à, t) = lim — = — = (4.60) 


Am>0Am dm pda’ 


The resultant stress vectors R; acting upon the differential element are found to 
be 


Rı — - (T'a  T"e3) y anda”, Ry = —(T"'e + Te) V anda’, 


— OR * [o * 
(r +da! |) =— (Te + Te) anda? — $ (T'a + T?) dodo! de’, 


a! 


— OR * [o * 
(s dira a) —— (Te, + T?) aj dal — e (Tey + T?) aj, do! de’. 
(4.61) 


Applying the law of conservation of momentum to Eqs. (4.59)-(4.61), for the 
equation of motion of the soft shell we get 
OR; 
dr? Oa! 0o? 


dv? 4- p 4- f p, (4.62) 
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2. = * 
where d7 is acceleration. Substituting R; and p given by Eqs. (4.58), (4.61) into 


(4.62), we get 


dr ð * ð " 

DEP 7 sd [ur + T&)\Van| + E [a + T?ej v an 
+ py à + fpya. 

(4.63) 


Let G;, M p and M '; be the resultant moment vectors acting on the element of the 
shell defined by 


G, =F x Rj, G =7Fx Ry 
— 0G, — 0G, i ] 
aq + Bal à), -(@ + e az) (a + da’ = constant) (4.64) 


M, = (F x pM dda de? My = (F sc Pip) dado 
Assuming the shell is in equilibrium, the sum of the moments vanishes. Hence 
0G, ôG 
Sal 6^ 0o? 
Substituting G;, M, and M; in (4.64), we obtain 


| ( x z ax" (7 x an) + (F x p) + (F x A) 


— (rx R,) da" — (rx Ry) da? = 0. 


da? + M, + M, =0. (4.65) 


Further, on use of Eq. (4.61), we find 


[e x (rz, + T?) + & x (THe + TPS) y aiy a» dala? — F 


OR 1. Oo, | LL 41 _ 
| mE s, d^ P+F = (4.66) 


Since the underlined term equals zero, we have 
(@; x &)T + (e x er?! = 0. (4.67) 


It follows immediately from the above T? = T?!. 
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4.5.1 Remarks 


1. If a soft shell is parameterized along the principal axes then T! = T), 


T? = h, T? = 0, aiz = adj» = O and the equation of motion (4.62) takes 
the simplest form 


dr 0 Fs 0 x n 
pvo = Gui In soe + oi nies + py a+ fpya (4.68) 


2. During the dynamic process of deformation different parts of the soft shell may 
undergo different stress—strain states. The biaxial stress-state occurs when 
1P = 7, + 7T)> 0, i = TiT, > 0, the uniaxial state develops in areas 


where I ie > o, = 0, and the zero stress-state takes place anywhere in the 


shell where im = v) = 0. The uniaxially stressed area (75 = 0) will 
develop wrinkles oriented along the action of the positive principal membrane 
force Tı. The equation of motion for the wrinkled area becomes 


dr 0 f + E 
pva s E In zl + pya + fpva. (4.69) 
1 


To preserve smoothness and continuity of the surface S, the uniaxially stressed 
area is substituted by a “smoothed” surface made out of an array of closely packed 
reinforced fibers. Such approach allows one to use the equations of motion (4.62) 
throughout the deformed surface. 

The governing system of equations of dynamics of the soft shell includes the 
equations of motion, constitutive relations, initial and boundary conditions. 


4.6 Nets 


A special class of soft shell, where discrete reinforced fibers are the main structural 
and weight bearing elements, is called nets. Depending on engineering design and 
practical needs, the fibers may remain discrete or embedded in the connective 
matrix. Although nets have distinct discrete structure, they are modeled as a solid 
continuum. Since the nets have very low resistance to shear forces, then Tj; = 0 
(S" = 0) and the resultant formulas obtained in the previous paragraphs are valid 
in modeling nets. 

Consider a net with the cell structure of a parallelogram. Let the sides of the cell 
be formed by two distinct families of reinforced fibers (Fig. 4.4). 
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Fig. 4.4 A structural 
element of the net formed by 
two distinct types of 
reinforced fibers 


0 
Let the undeformed (S = S) configuration of the net be parameterized by 0, d 
coordinates oriented along the reinforced fibers. For force distribution in the net 


we have to 


1. Find the stretch ratios 2; = A41, 45 = do, using Eq. (4.32). 


2. Calculate forces in reinforced fibers:T], T5, 
3. Find the membrane forces in terms of S—configuration 


: A rA. a0 0 Ai 5. 
pu ou (ri 72 sn (d - 9) + Ti cos (d — V) 
2 


sin y 4» a 

2 1 V 0 0 Vu 

T? = 7 (n z sin? (% V) siny — T; cos — y 
sin y 1 2 

7 À VH 0 A, 0 

T? — z (n ; sm y HT) 7 cos? 2 
sin y 4 a 2 


) cos i). (4.70) 


The principal membrane forces and their directions are found from Eqs. (4.54), 


(4.55) by putting T? = S = 0, T! = T7, T? = T} 


(rt 1) + (rr - T) eermieos ; 


Ti» = 


2sin y 
* TS sin2 y 
tan2 y, ME ia NN 
Ti + T5cos2x 
* * 
V, = *n/2. 


In particular, 


1. If y = z/2then y, = 0, Ti = Tt, T) = T}; 
2. If T! = 0 then, = 0, Ti = 24, T = 0; 


sin y 
3. If T; = 0 then ý; = % Ti = —, T; = 0. 


sin y 


(4.71) 
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4.7 Corollaries of the Fundamental Assumptions 


We conclude this chapter with several corollaries of the fundamental assumptions 
stated at the beginning of this chapter that are specific to thin soft shells. 


1. The zero-flexural rigidity state is natural and unique to thin soft shells in 
contrast to thin elastic shells with finite bending rigidity. 

2. Soft shells do not resist compression forces and thus 7; > 0, 75 > 0 and 

T T 
KD >0, f » o. 

3. Shear membrane forces are significantly smaller compared to stretch forces, 
Ti ~ 10? maxT;. 

4. Areas of the soft shell, where A; x 1 and A» x 1, attain multiple configura- 
tions and are treated as the zero-stressed areas. 

5. Stress states of the soft shell are classified as: (1) biaxial, if T; > 0, T > 0, 
a > 0, 2 > 0), (2) uniaxial, if either T, = 0, T; > 0 or T; > 0, T; = 0, 
ad” > 0, i” — 0) and iii) unstressed, if T, — 0 and T>=0, 
at = 1 = 0). 

6. Constitutive relations for the uniaxial stress—strain state are functions of either 
A, or Ay and empirical mechanical constants cm given by 


Ti = fil, V5 61v Zi) for Ay > 1, Ao< 1, (4 72) 
T, = fo (Ao, V5 C), . Cm, Zi) for ^i < 1, A, > 1. f 


7. Constitutive relations for the in-plane biaxial state, A; > 1, A > 
1 (Tj > 0, T > 0) have the form 


Tı = Fi(Ai, Ao, v5, C]; - - Cg; Zij), 
T? = Fo(Ai, Ao, V, c1,- . Cm, Zij), (4.73) 
y = V(Ai, Ao, V5, Cis. Cm, Zij). 


In general f,(---) Æ F,(--), however, f,(--:) can be defined uniquely if 
F,,(...) is known. 
For example, Corollary 6 of the fundamental assumptions for the nets is given 
by 
IM =T +r = 


(Ti + T3) 20, 


sin% 
IP =T =T eg 
Qu i ee Rone 
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Chapter 5 
Continual Model of the Detrusor 


There is something fascinating about science. One gets such 
wholesale returns of conjecture out of such a trifling investment 


of fact. 
M. Twain. 


5.1 Basic Assumptions 


Histomorphological studies have revealed the exact orientation of smooth muscle 
fibers in the bladder wall. The distinct inner layer has circumferential, and the 
outer layer has longitudinal orientation. Three-dimensional structure analysis of 
the organ suggests smooth muscle fasciculi as functional units of the organ. They 
merge, dichotomize, and intertwine with each other to form an interlacing net- 
work. It is reasonable to assume that the detrusor is a nonlinear heterogeneous 
syncytium and possesses general properties of curvilinear orthotropy with axis of 
anisotropy defined by the orientation of smooth muscle fasciculi. Being hetero- 
geneous, nonlinear, viscoelastic, incompressible composites, detrusor defies sim- 
ple material models. Accounting for these particulars in a constitutive model and 
both experimental evaluations is a great challenge. 

The complete theoretical formulation of a mathematical model of the urinary 
bladder wall tissue is best achieved with application of the principles of thermo- 
dynamics supported by extensive experimentation. The advantage of such 
approach is that it employs generalized quantities that allow tackling specific 
problems related to discrete morphological structure of the tissue and the con- 
tinuum scale of description, which is typically ~ 1 um. We shall base our further 
derivations on the following basic assumptions (Usik 1973): 


1. The biomaterial (the wall of the urinary bladder) is a two-phase, multicom- 
ponent, mechanochemically active, anisotropic medium; phase 1 comprises the 
smooth muscle fasciculi, and phase 2—all remaining structures; 

2. The phase interfaces are semi-permeable; 

3. Active forces of contraction-relaxation produced by smooth muscle are the 
result of intracellular mechanochemical reactions; the reactions run in a large 
number of small loci that are evenly distributed throughout the whole volume 
of the tissue; the sources of chemical reagents are uniformly dispersed within 
the volume of the biocomposite and are ample; 
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4. The biocomposite endows properties of general curvilinear anisotropy and 
viscoelasticity; the viscous properties are due to smooth muscle fiber mechanics 
and the elastic properties depend mainly on the collagen and elastin fibers; 

5. There are no temperature and/or deformation gradients within the tissue; 

6. The biocomposite is incompressible and statistically homogeneous. 


5.2 Model Formulation 


All derivations to follow will be obtained for the averaged parameters. Here we 
adopt the following notation: the quantities obtained by averaging over the volume 
of a particular phase are contained in the angle brackets, and those free of brackets 
are attained by averaging over the entire volume. 

Let p be the mean density of the tissue. The partial density of the Cth substrate 
(¢ = 1,n) in the f phase (f = 1, 2) is defined as 


pt = mb ["V, 
where mf is the mass of the Cth substrate, mY is the total elementary volume of the 


£ —B 
tissue "V = ares "V . The mass and the effective concentrations of substrates 
are 


g 
ct = pf /p, (h) = mE" V of. (5.1) 


Assuming p = (pP) — const, we have 
(0°) = > (5.2) 


Setting f = 1 we find 


"S 
where y is the porosity of phase fj (n — "V /"V ). It is easy to show that 


n=) cac. (5.3) 
c1 


Hence the mass e and the effective < ) concentrations are interrelated by c = 


n(c j The sum of all concentrations er in the medium equals one (5 7 , e =1). 
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Change in the concentration of constituents in different phases is due to the 
exchange of the matter among phases, external fluxes, chemical reactions, and 
diffusion. Since chemical reactions run only in phase 2 and the substrates move at 
the same velocity, there is no diffusion within phases. Hence, the equation of the 
conservation of mass of the (th substrate in the medium is 


dc! 
pol UT ar = 08+ Yov (5.4) 


Here o! is the velocity of influx of the ¢th substrate into the phase o, v;jJ; is the 


rate of Cth formation in the jth chemical reaction (j = 1,r). The quantity vg is 
related to the molecular mass M; of the substrate ¢ and is analogous to the stoi- 
chiometric coefficient in the jth reaction. vz; takes positive values if the substrate is 
formed and becomes negative if the substrate disassociates. Since the mass of 
reacting components is conserved in each chemical reaction, we have 


n 
, Vg =0. 
el 


Assume that there is a flux of of the matter into: (1) phase 1 from phase 2 or the 
“distributed” external sources and (2) phase 2 from phases 1 only. Hence, we have 


O; = —Or + OF, OF = Or, (5.5) 


where Q* is the flux of distributed sources, Qç is the exchange flux between phases. 
Applying the incompressibility condition to (5.5), we have: 7 4 Of — 0. Addi- 
. n 
tionally, let also 5 /* ., Qg = 0. 

Assuming that the effective concentration of substrates remains constant 


throughout and neglecting the convective transport of matter within phases, then 
with the help of Eqs. (5.5) -(5.4), we obtain 


Oc} : Oc? r 
Pa =- Ob ose T QUE v (5.6) 
j=l 


The equations of continuity and the conservation of momentum for the tissue 
treated as a three-dimensional solid in a fixed Cartesian coordinate system is given 
by 


Ou; Ôo 


Pep = Ox; + A Sn 


Here u; are the components of the displacement vector, f; is the mass force, and 
gj; is the stress tensor for the whole medium. 
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Let UP, sP), c; be the free energy, 
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entropy and stresses of each phase. The 


following equalities are an extension of the assumption v 


gets 


U 


oy = elo") + 


+(1- c')(s*) 
eg +(1- eU 


(1—c!)(o),. 


The Gibbs relation for each phase is defined by 


1 
c'(U'} = Up (ct, T) t 2p Piintütim (5.8) 
1 
(1 = c!) (U?) = US ef, T) + 5 Yimin (5.9) 
d(U*) == (0!) de} — (s')aT + D> (ut )d(c!) (5.10) 
AU?) = (0?) dc} - ()ar +) (jale) (5.11) 


ç=1 


(ul) = o(U*)/0( cf), (sf) = 


where ut is the chemical potential of the Cth substrate in the fj (f = 1,2) phase, 


e(u^y/oT, 


ub = wade cj; is the elastic (phase 1) and Aj; is the viscous (phase 2) part 
of deformation ie = =g + Ab ij); E ijkl, Yijxı are the fourth rank tensors, T is tem- 
perature. Making use of the identities 


pol B B 
Oc P(UPY [ce = (ut )+ + (UP) — 2:004) 
Equations (5.10) and (5.11) can be written in the form 


d(c'(U')) = Tela!) dal = c! (s! dT + 5 ude; 


> 3 (5.12) 


d(c2(U2)) = Soe) dsj - 2(s2)aT + Y dad. 


ç=1 
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Assuming that the mass sources are present only in phase 1, the general heat 
flux and the second law of thermodynamics for the tissue are described by 


1 b Aor 
dU — 7 cjdej — sdT — dq’ + i 3 acl Qtdt, (5.13) 
e / = 1 Q; 
Tds = dq’ + dd +5 TS! d, (5.14) 


ç=1 


T- as u OF 
real = sro’ 
SS Te (OAL), 34.8 . 


Here U — S cP (UP), s = Tai cP(sP), and St is the partial entropy of the 
tissue. On use of Eqs. (5.12) in (5.13), (5.14), the equation of the balance of 
entropy of the detrusor takes the form 


ds ds . q] R 
p ; Qt = -divz (5.15) 
dt ra Oc? ~> T T 
d TN ME 
Re Bade. +> (i -4)9: - 34; ply) 
ç=1 j=l 
Aj = $ vak (5.17) 
(= 


Here R is the dissipative function, Aj is the affinity constant of the jth chemical 
reaction. 
Let the generalized thermodynamic forces and the heat flux, g, be defined by 


how 1 dA; (ui - 22) Aj 
EI LA LE MEE x (5.18) 
pdq) = —divgdt. 


Then, the thermodynamic fluxes q, Qf, Qç, J; and stresses 0; can be expressed 
as linear functions of the generalized thermodynamic forces 
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oT 
qi = —Ny—, 
j 
1 
Oy = Ejuéij, 
P dA. xr 
Tu = Bini Ge - 2 Pubs + So Mu =e) Sis 
a=1 5.19 
Jy = Dy d ee up (Hy — 1), 


Qa = MagN, — > “Lapp + 2’ aß (uj - iis). 


= (v) (=) 
Here "/,5.5, (m = 1,3) are scalars, Bj, Dij; Maij, Nij are the parameters of ten- 


sorial nature that satisfy the Onsager reciprocal relations 
Biu = By, Da -D pij Nj; = Nij, M uj E -Mj lap py =" lhah, 


and of, (B = 1,2) are the components of the total stress tensor, ox, of the 
biocomposite 


On = E que; + EAE 
where Ex Eig are the tensors of elastic and viscous characteristics (Eija = EX) 
of the detrusor, i.e., collagen and elastin fibers of the stroma and smooth muscle 
elements, respectively. The term Xg- D puAg is the biological input—“biofac- 
Or", Z, and accounts for biological processes, e.g., electromechanical, chemical, 
remodeling, aging, etc., in the tissue. With the help of (5.17) and remembering that 


(ut) = a(UP) (ef), it is found to be 


nx oU? 
Zu == S Dm, = » Y Paras Sex ac? 2 è (5.20) 


p=1 ß=1 a=1 


Thus, the system, including equations of conservation of momentum (5.7) with 
ob, defined by (5.19); conservation of mass of the reacting components (5.4) 
together with (5.16) and Q, defined by (5.19); heat influx which is obtained on 
substituting s = —OU/OT into Eq. (5.16) and the explicit expression for q, 
describes the continual mechanical model of the detrusor. However, it does not 
guarantee the required relationships between in-plane forces and deformations in 
the soft shell. To establish the missing link, recall the second Kirchhoff—Love 
hypothesis, which states that the normal stress o33 is significantly smaller 
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compared to oj (i,j = 1,2). Then, the terms containing 033 can be eliminated. 
In general, the end formulas are very bulky and are not given here. In applications 
though, depending on a specific tissue, the formulas can be simplified to a certain 
degree and even take a compact form. 


5.3 Biofactor Zy 


The fundamental function of the bladder—filling and voiding—is closely related 
to electromechanical wave processes and the coordinated propagation of the waves 
of contraction-relaxation in the organ. 

Consider smooth muscle syncytia to be an electrically excitable biological 
medium, (Plonsey and Barr 1984). Applying Ohm's law we have 


j= - (auia, +39 9T (5.21) 
cS 8il m ei 8i2 Ax e2 j, ‘ 

_ . OY. . OP, _ 

Jo=- (2. & 4 + 2o2 ca) ; (5.22) 


where J;, J, are the intracellular(i) and extracellular (o) currents, V;, V, are the 
scalar electrical potentials, 2;;, ĝoj (j = 1,2) are the conductivities, and e;,e» are 
the unit vectors in the directions of «1, o» coordinate lines. Both cellular spaces are 
coupled through the transmembrane current J; and potential Vm as 


Imi = —divJ; = divJ,, (5.23) 
Vin = Yi — Yo. (5.24) 
Substituting Eqs. (5.21) and (5.22) into (5.23), we get 


Ow; Ow; 
Im = ĉn ——-61- $0 ——-62, 5.25 
1 dE €i t 8n oo e ( ) 
PF, A, 
deeem um debo m. 5.26 
1 Bol Bag Oi + Bo Bop e2 ( ) 


Equating Eqs. (5.25) and (5.26), we find 


Ow, OY, ev, ev, 
pup e Eae dio. =, ĝo . 5.27 
(in + 81) Gaz ^ (8i2 + Ea) Qo Eol ag? + 802 G3 ( ) 


Solving Eq. (5.27) for V;, we obtain 
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AE Bol a3 o2 y 
~ an Sil T 8o1 w ên F o2 oY? 
/ $ 2 1 $ 
+ hoe( (x - x)' +(x- Y) ) exar. 
where the following substitutions are used: X = o / gii ^ goi, Y = %/ 
vên + ĉo2. Here the integration variables are primed, and the unprimed variables 


indicate the space point (o, o5) at which V; is evaluated. The reverse substitutions 
of X and Y give 


eV 
br -x ll (ea ote ox? + Bo ze) 
, 2 
+ llog (a, — a) 4 (% T 2;) i oe — 
gil + ĉo 812 + £02 VIT + Bo1)(Bi2 + B02) 


Introducing Eq. (5.28) into (5.25), after some algebra we obtain 


2 
oda 0 Vm 
In = o 7 > + D 7 
1 2xG + ju) ja) TIC 2 + 802 at) 


2 1\2 1\2\ 2 
a, —o Ly — q ay — a Uy — q ie act 
x (ec aad) (B cs) +e) da, doo, 


(5.28) 


here 
Gi = $1 + Bo1, Go = Ên + Bor, 
= y GiG, Hy = 10/8; Hy = Bo2/8i2- 
Substituting Eq. (5.28) into (5.25), we find the contribution of an e-neighbor- 


hood of CA =0,% = 0) to mı. Using the transformations given by X = 


a / Vn, Y = */ /80, we find 


Vign Vm ~ E YS 
4xG (Ê aya t 82 gY? da 


x* y? P 
x V? do — + - )je da,. 
/ | (am G2/ dis 


Here V? is the Laplace operator. Applying the divergence theorem and 
performing the gradient operation, the integral in Eq. (5.30) is converted into a line 
integral 


Im = 


(5.30) 
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- üdC , (5.31) 


J (2X 8i1/G1)@ + (2Y 82/G2)e» 
(X285/G1) + (Y?852/G2) 


where dC’ is an element of the e-contour. The result of integration yields 


Ov, O^ VN (ên 2(], — iia) Gi 
In = ĉo 7 Bo 7 : E = t a — 32 
: (sn Get ee SP) (B+ a em m 2 pus 


To simulate the excitation and propagation pattern in the anisotropic smooth 
muscle syncytium we employ the Hodgkin-Huxley formalism described by 


OV, 
Cm 27 = —(Ini + In? zs Lon), 


where Cm is the membrane capacitance, Jion is the total ion current through the 
membrane. Substituting expressions for /mı and In. given by (5.29), (5.32), we 
obtain 


OVin — I mm Hy ^ 0^ V, 
Cop = 2nG(1 + fy) 4- all (m4 bol By ae EC oy? 
2 2 (X2 (122 
- (a — a) (o5 — as) / (a — 94) ! (a5 — o) dz, das 
Gi G5 Gi G2 


zo P Va o 8 Val fr ie jh G 
(2 19 t 802 /2 2) (&. t T Ha) tan ! ) xz Tion; 
Oot; 02 GP m(l- Hy )(1 + fy) G2 


where Jion is the function depending on the type and ion channel properties of the 
biological tissue. 

In the case of electrical isotropy, 1, = M = Ji, the integral in Eq. (5.33) 
vanishes and we get 


2y 2y 
Ca 3E --g : " (i A e da a t =) i, (5.34) 

Finally, the constitutive relations of mechano-chemically active electrogenic 
detrusor include equations of the continual mechanical model as described above, 
Eqs. (5.33) and/or (5.34). The system is closed by formulating the free energy, ion 
currents, initial and boundary conditions, and the function Z; = Zi (Vm. His ij; ĉoj). 
It is noteworthy that the closed system of equations describes the development of 
forces in absence of active strains and vice versa, a condition which is unique to all 
biological materials. 

Models are generally evaluated for the degree of parameters and constants 
involved as well as for their accurate and meaningful experimental determination. 
The constitutive model as described above includes 6 +n + r+ 'h(n? +r? + 
n + r) independents scalar quantities. To obtain their experimental estimates will 
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always remain a great challenge. Therefore, simplified models and descriptors are 
required to serve specific needs of the investigator. While phenomenological 
constitutive models are able to fit the experimental data with a high degree of 
accuracy, they are limited in that they do not give insight into the underlying cause 
to the particulars of mechanical behavior. Fine molecular and structure-based 
models help avoid such ambiguities and are able to reveal the intricacies of 
functions of tissues. However, they are beyond the scope of this book. 


5.4 Special Cases 


Constitutive relations for soft biological tissues are usually obtained along struc- 
turally preferred directions that are defined by the orientation of reinforced smooth 
muscle, collagen and elastin fibers, and thus, make them easy to use in 
calculations. 

However, if constitutive relations are obtained in the directions different from 
the actual parameterization of the shell, the task then is to calculate membrane 
forces in the principal directions. Consider two typical situations: 

Case 1. Constitutive relations are given by Eq. (4.73). Then, 


1. From Eqs. (5.39) and (5.40), we calculate the principal deformations A,, ^» 
and the angle qj, 

2. Using Eq. (4.73), we compute the principal membrane forces T; and T, and the 
angle y, 


3. Finally, setting 7, = 2/2, x; = y, T! = T;, T? = T>, T? = 0 in Eq. (4.49), 
we find 


T! = n sin?(y —W+tTh cos? (y — V) / sin y, 
T? = (T, sin y sin(y — y) — T cos y cos(y — V) / sin y, (5.35) 
T? = n sin? y + T> cos? y) / sin x. 


Case 2. Constitutive relations are formulated for the orientation of reinforced 
fibers, superscript (7), 


Ti = Fi (21, 45, y, c1. Cm, Zij), 
T5 = FAA}, 45, Y, €1, + - -Cm Zij), (5.36) 
v= S'(A, 45, y, C1,.. Cm, Zi): 
* 
Let di € S be an auxiliary orthogonal coordinate system oriented with respect to 


a set of reinforced fibers by y. Then, 
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1. petting y1 = 8/2, x4 = Y, X1 = n/2 — y" in Eq. (5.32), where 41 := 41, 45 := 


p for the stretch ratios and the shear angle y", we have 
0 0 
i= G sin?(7 — V) + 23 sin? y 
o., 0 0 gs 0 
--24443 cos y sin( y — V/) sin v) /sin x 


1 1 0 
y = sin! (— (- z4 sin 2(y 7 y) 
A142 sin? X (5.37) 
1 2... 0 ‘ 0 . ,0 0 
+ 5^ sin 2 y +2122 cos y sin(y — V) 
o 0 0 
As = (2 cos?(y — y) + 45 cos? y 
o 0 o 0 I2 0 
—244142 cos(y — W) cos y cos v) /sin x 


2. Using Eq. (4.73) we find T], T5 and S", 
3. The angle y is found from Eq. (5.32) by putting h V, % 0, 4 = À, 


4 r 
M = A, 


* E 1 . 40 0 . Oo. 0 
V = cos ee A, sin(y — V) + 42 cos y sin y | |, (5.38) 
| sin y 


4. Setting y, = 1/2 — y", y, = V, I := pu, D i -= 72, S = -72 in Eq. (4.49), 
we obtain 


Ti! -mn sin 2(y — V) + Ti cos? (x = pay) 
—28" cos(7 = +y") sin(y — j o sin 7 cos 7”, 

pm =fr sin usin(%— V) — T; cos(W +7") cos( — V) (5.39) 
+S bin —2 y +y") — cos y sin r| n sin 7 cos 7”, 


= * * * * $ 
pe = sin? y —25" sin( —") sin y +T} cos? (y - sin y cos y”. 
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Formulas (5.39) can be written in more concise form if we introduce gen- 
eralized forces defined by N* = pa A. 


N! 1 T” A5 NZ Ti A NT S" 
, 2 25 , 
Then, Eq. (4.50) takes the form 
ik 1 a A xi Axk 
N" = 5, N" C, C,- (5.40) 
Substituting n given by Eqs. (4.16) and (4.18), we find 
N" ={w" si? ({— Jo) +N? sin? (J +p - Y) 
-2 N” sin(7, + % — 7) sin(7 — ia) n sin % sin A 
N”? ={W" sin sin(7 ia) 4a? 2 sin(%ı +2 — z) sin(%ı + %) 
(5.41) 


+N? [eos (i +2 %2 — 7) = cos % cos A ! sin 7 sin 7, 
N? ELA sin? 2 -2N" sin (ži + ia) sin Xo 
HS sin? (4 + ia) n sin x sin 7, 
Putting " = n/2, % I J in Eq. (5.41), for the membrane forces in terms of the 


0 
undeformed surface S(S = S), we have 


. À A . 0 0 k AT 0 0 
quit! (n T sin? (Z — Y) + T} z cos^(y — y) 
1 2 


pao 40 9 . 0 
—25S" sin2(y — y) | / sin y, 
` 0 0 gr 0 0 
T? = (rsm G- y) sin y —7; Z+ cos? (7 — V) cos y (5.42) 
1 2 


"UNE 0 
+S" sin(y — 2y) | / sin y, 


m A A A, 
T? = a (n 7 sin 2y ud J +S" sin2 i) / sin 2. 
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Formulas (5.42) depend only on parameterization of S and the axes of anisotropy 
and are less computationally demanding compared to Eqs. (5.38), (5.39). 
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Chapter 6 
A Model of the Detrusor Fasciculus 


One never notices what has been done; one can only see what 
remains to be done. 
M. Curie. 


6.1 Formulation of the Model 


Consider a detrusor fasciculus embedded in the extracellular matrix of connective 
tissue. Assume that (see Chap. | for details): 


1. 


2. 


Ww 


Smooth muscle cells are arranged in a soft fiber thread with properties of an 
electromechanical biological continuum; 

The myofiber has nonlinear viscoelastic characteristics; inactive smooth muscle 
cells, collagen, elastin fibers of the extracellular matrix define the “passive”, 
T?(A, c;), and intracellular contractile proteins describe the “active” compo- 


nent, T? (a, zie , [Cart]. a), of its total force, T" 


7 = TUA cj) + T(4, zm [Ca?*], eh: (6.1) 


where / is the stretch ratio, c; are empirical material constants, ze is the 
“biofactor”, and [Ca?*] is the concentration of free cytosolic calcium; 


. Contractions of the myofiber are isometric and deformations are finite; 
. The fasciculus possesses cable and myogenic electrical properties; myogenic 


electrical events are a result of activity of the fast (T-type) and slow (L-type) 
inward Ca", BKc,, voltage-dependent K}, and leak Cl” currents; the fasciculus 
is in the silent state and the transformation to a firing state is a result of 
depolarization and/or stretch deformation that alters the conductance for L- and 
T-type Ca?* channels, while the stretch affects permeability of L-type channels; 


. A smooth muscle cell or a group of cells within the fasciculus has intrinsic 


pacemaker properties; additionally, an a priori defined “pacemaker” may 
provide an excitation to the fiber. 


Let the fasciculus of a length L be referred to a local Lagrange coordinate 


system «. Its equation of motion is given by 
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Qo 0 
~= T, (O<a<L 6.2 
pa =T, OSaSL) (6.2) 
where p is density, v is the velocity, and the meaning of other parameters are as 
described above. Following the working assumptions, the total force T' can be 
decomposed as 


o(A — 1) 


T =k, 
Ot 


T (129), [Ca], «) + T". o), (6.3) 


ie mn? 


where the viscoelastic term has been added to Eq. (6.1). Here ky is viscosity of 
tissue. Substituting the above into Eq. (6.3) we obtain 


O 08 f, 0-1 
Ps =~ (x ( : ip (4, 24), [ca], «) + 7°(2, «)) (6.4) 


The force-stretch ratio relationship satisfies Eq. (6.1) that is used to fit the 
experimental results on uniaxial stretching of linear strips of the tissue 


p [jelexpe(4—1)-1], 42 1.0, 
iiis n otherwise. (6.5) 


Although the exact data on the sequence of mechanical events in the longitu- 
dinal and circular muscle syncytia during contractions of the detrusor are not 
available, it is reasonable to assume that both muscle syncytia contract simulta- 
neously. Such coordination results in generation of the strong active forces T?, that 
are needed to effectively expel urine from the bladder. This assumption, together 
with the fine fibrillar structure of smooth muscle, suggests that 72; are produced 
only in the preferred directions, either longitudinal or circumferential, and as such 
can be characterized in full by uniaxial tests. Further, constructive modeling 
requires formulation of the excitation-contraction coupling phenomenon that 
provides a link between electrical events and the muscle mechanics. The key 
player in the dynamics of the transformation of microscopic properties of electrical 
excitatory events into macroscopic contractions is free intracellular Ca”. 
Although many details of this process have become well established within the 
past decade, a convincing experimental correlation between ultrastructural chan- 
ges and force development has not been clarified yet. 

As a first approximation of the excitation-contraction phenomenon we adopt 
the active force—intracellular Ca^* relationship given by 


0, [Ca2*] < 0.1 uM 

c; + cs[Ca2*]* + co [Ca2*? + eio [Ca2^ ? + c11 [C2], 
0.1 < [Ca2*] < 1uM 
max T“, [Ca2*] > 1uM 


T= 


ol. 


(6.6) 
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where [Ca?*] is the intracellular concentration of calcium ions. 
The system of equations for the oscillatory activity of the membrane potential V is 


4C,,— = p» : (6.7) 


where, 4 is the numerical parameter, C,, is the smooth muscle cell membrane 
capacitance, and J; is the sum of the respective ion currents 


I. SUL (V — Voa); 

Ra =8tažca (V ^ Vea) 

Ii =8xiit"(V — Vi), 
i _8ca [C8] (V — Vaa) 
up 0.5 + [Ca] 

la =8a(V 7 Va). 


Here Voa, Vki, Vc, are the reversal potentials, and gh,, g$, gk, 2c, 4o 8c] are 


(6.8) 


T 


the maximal conductances for the ion currents, m,h,n and Xe, are dynamic 
variables described by 


Xm 


dm + Bm 


dh . "m 
A —&,(1— h) — B,h, 


da x: 
an =4 (128) — Bh, (6.9) 
y dxc, 1 = 
T Cs = Xe 
* dr ~ exp(—0.15(V + 50)) “© 
acaz] 


ee 2+ 
A dt =P caca (Voa - V) ES [Ca; ib 
where the activation à, and deactivation P, O =M, h, n) parameters of ion 


channels satisfy the empirical relations 


~ |. 0.1(50-¥) u d) 2 
Ln = exp(5-0.1V)-1" Bn = 4 exp = — 
A 25-0.1V z 
3, = 007 ea, hoic; (6.10) 
~ | 0.01(55-V) = xp V) 
On = ae(55-017)-1? B, = 0.125 exp. 


Here V = (127V + 8265) /105, Txca is the time constant, gc, is the parameter 
referring to the dynamics of Ca** channels, fi is a numerical constant. 
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The evolution of L- and T-type Ca**-channels depends on the wave of depo- 
larization, V*, and is defined by 


gealt) =[0(V) + (A(t) — 1)](maxge,), (6.11) 
NO =(A(t) xw 1)g6. 
where 
5 1, fr V>VS 
— 1 — P 
AG 2 OE) { 0, otherwise 
Here V5 is the threshold value for V*. 
The propagation of the wave of excitation V? is described by 
OVS dm 0 ovs 
C = A Ina + Ig2 + f 6.12 
" Oy Pag (005) Una + tka + Te) oY) 


where dm is the diameter, R, is the specific resistance of the fasciculus, and 


Iya =gyañ h (V5 = Vna) 
Ik —gxofi* (V° — Vra) (6.13) 
Ia =8a(V* — Va). 

Here 24,,2y5,2c, are the maximal conductances, and Vha; Vko; Vc, are the 


reversal potentials of Na*, K$, and Cl” membrane currents, respectively. The 
dynamics of the variables ñ, h, ñ are described by 


dt = mi m) Be 
dh. re 
T =ô, (1 — Â) — Bh (6.14) 
dà |, TEM 
dt —á,(1 £z f) = p, 


with the activation &, and deactivation By (y = M, h, ñ) parameters given by 


0.005(V—V,,) = (VV) 
Xm = exp0.1(V—Vn)—1? Bn => 0.2 exp 38 ; 
K —(Va+V" ^ 
y = 0.014 exp Ust"). f, = pU (6.15) 
0.006 (V*—V,) 


EE ELLA EA HEC f — _ ys 
"^ exp0.1(V8=Va)—1’ fj, = 0.75 exp(Vn = V^). 

Here Vm, Vn, Vn are the reversal potentials for activation and inactivation of 
Na* and Kj? ion currents. 

In the following numerical experiments, we assume that at the initial moment of 
time the functional unit is in unexcitable state 
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V*(a,0) = 0, v («,0) =0,|Ca*| = [ca] (6.16) 


L Lh > ~ ~ __ x00 
h= ho: N= Ny, Xo, = ea 


It is activated by a series of discharges of externally applied action potentials 


0 
V*(0, t) = D e E V*(0, t) = V(t), (6.17) 


The ends of the myofiber are clamped and remain unexcitable throughout 
V*(0, t) = V*(L,t) = 0, v(0, 1) = v(L, t) = 0. (6.18) 


Equations (6.4)-(6.15), complemented by initial and boundary conditions 
(6.16—6.18) constitute the mathematical formulation of the model of the electro- 
mechanical activity of the fasciculus. It describes: 


. Self-oscillatory behavior and myoelectrical activity of the fasciculus; 

. Generation and propagation of the wave of depolarization along the myofiber; 

. Generation of action potentials (spikes); 

. Dynamics of the cytosolic Ca? transients; 

. Active and passive force generation; 

. Deformation of the fasciculus and the following excitation of the cell mem- 
brane with subsequent contractions. 


NNnBWNY Re 


The governing system was solved numerically using ABS Technologies® 
computational platform. It employed a hybrid finite difference scheme and finite 
element method of second-order accuracy, with respect to spatial and time vari- 
ables. The parameters and constants used in simulations were derived from the 
published literature. The values, which could not be found in the literature, were 
adjusted during the experiments in order to mimic closely the behavior of the 
biological prototype. 


6.2 Physiological Condition 


The resting membrane potential of the unexcited fasciculus is V' = —45 mV. 
Continuous fluctuations at low rate and amplitudes of the T-type Ca?* —0.006 + 
0.01 nA, the BKc, —0.05 — 0.11 nA, and the small chloride —0.04 nA currents 
result in oscillations of the membrane potential V which resembles slow waves. 
Their frequency, v = 0.5 Hz, and the amplitude, V — 6 + 10 mV, vary in time. 
The maximum rate of depolarization is dV/dt = 10 mV/s and of repolarization is 
dV/dt = — 7.2 mV/s (Fig. 6.1). 

Slow waves induce the flux of Ca? ions inside the cell at an average rate of 
0.04 uM/s. There is a short 0.65 s time delay in the intracellular calcium transients 
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as compared to the wave of depolarization. Free cytosolic calcium at max 
[Ca?*] = 0.42 uM activates the contractile protein system with the production of 
active contractions, max7* =1.2 mN/cm. They follow in-phase and time the 
dynamics of calcium oscillations and are normally preceded by slow waves. 

Electrical stimulation of the detrusor myofiber triggers the generation of high 
frequency and magnitude ion currents: J$, —1.1, IL, =0.038, Z1 =7.6, 1c, y —1.5, 
and la — 0.4 (nA), respectively. Action potentials of amplitudes 35 — 40 mV ata 
frequency ranging between 4 and 6 Hz occur as single discharges or bursts of 4—10 
spikes. Each spike consists of rapid depolarizing, repolarizing followed by after- 
hyperpolarizing (AHP) phases. The rates of depolarization and repolarization are 
dV /dt = 7.5— 8 mV/ms with the half-duration of a spike ~6 ms. Amplitudes of 
fast AHPs vary between 8 and 12 mV (see a fast timescale of Fig. 6.1). A con- 
comitant rise in intracellular calcium to 0.45 uM leads to the development of active 
force, max7" —1.4 mN/cm. Upon the termination of electrical discharges the 
myofiber returns to its unexcited state. 


6.3 Effects of Pharmacological and Extracellular Ion 
Changes on Electromechanical Activity of the Detrusor 


6.3.1 Changes in Kj 


A four-fold increase in the concentration of extracellular potassium lastingly 
depolarizes the detrusor cell membrane, V' — —17 mV, and abolishes its oscil- 
latory activity (Fig. 6.2). It affects the transient dynamics of the fast Ca** as well 
as BKc, and K} channels. They remain open and active throughout the time of 
depolarization. The strength of T-Ca?* current is E, = 0.028, Tey = 0.03, and 
doce — 0.06 (nA), respectively. There is a significant drop in the free intracellular 
calcium level, [Ca?*] = 0.35 uM, and an associated reduction in active force of 
the myofiber, T* = 0.9 mN/cm. 

The following gradual decrease in [Kj ] hyperpolarizes the detrusor membrane: 
V" = —25, —36, —38 (mV). The intensity of the outward potassium current 
increases to 1.2 nA with a concomitant attenuation of the T-Ca?*, i, = 0.03 nA, 
and the calcium-activated potassium, Ts. Se = 0.02 nA, currents. The fast and slow 
Ca** channels regain their dynamics and detrusor starts generating slow waves of a 
small amplitude, V = 3 mV, and a frequency, v = 0.13 Hz. There is a steady rise in 
[Ca?*] to ~5.0 uM along with the intensity of contractions, maxT* = 1.59 mN/cm. 

Interestingly, an electrical excitation does not lead to a production of spikes. 
Only after [Kj ] returns to its physiological values does the fasciculus regain 
normal electromechanical activity. 
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Fig. 6.1 The dynamics of ion currents, electrical pattern, intracellular Ca?* and total force in the 
detrusor fasciculus following excitation of the smooth muscle membrane 


6.3.2 L- and T-type Ca?* Channel Antagonists 


Treatment of the fasciculus with a selective T-type Ca?* channel blocker, mibe- 
fradil, produces dose-dependent reduction of the fast inward calcium current, the 
frequency of spikes, calcium transients, and the intensity of contractions (Fig. 6.3). 
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Fig. 6.2 Electromechanical response of the detrusor to excess extracellular K* 


The drug at “low” dose decreases the amplitude of the fast T-type Ca?* current to 
0.036 nA. The frequency of slow waves and action potentials diminish to 0.2 and 
1.5 Hz, respectively. The detrusor produces bursts of 2—4 actions potentials at a 
time. 

The main influx and rise in [Ca?*] occurs during the times of stimulation. The 
maximum concentration of intracellular calcium reaches 0.45 uM and the active 
force —1.41 mN/cm. At "high" concentrations mibefradil further reduces the 
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frequency of spikes, however, it does not affect the pattern and strength of phasic 
contractions. 

Nifedipine, a nonselective T- and L-type Ca** channels antagonist, hyperpo- 
larizes the cell membrane and shift the resting membrane potential to V' = 
— 52 mV. There is an increase in the amplitude, V = 14 mV, and the frequency, 
v = 0.7 Hz of slow waves (Fig. 6.4). Electrical stimulation in the presence of 
nifedipine induces high amplitude action potentials, V = 38 mV, although, at a 
lower frequency, v = 3.5 Hz, compared to the physiological norm. No changes are 
seen in the dynamics of j A decrease in the amplitude of the fast inward calcium 
current J£, = 0.025 nA results in a decrease in the level of free cytosolic calcium, 
[Ca?*] = 0.16 uM. The myofiber continues to generate phasic contractions of 
0.86 mN/cm. 

Application of verapamil, a selective L-type Ca^* channels blocker, inhibits 
spontaneous and induced myoelectrical activity in the detrusor. The intracellular 
calcium level decreases to 0.01 uM and no contractions are produced. 


6.3.3 BKc, Channel Agonist/Antagonist 


Application of forskolin (FSK), a Ca” -activated K* channel agonist, dose- 
dependently hyperpolarizes detrusor: V = —62, —68 (mV) and abolishes slow 
wave oscillatory activity (Fig. 6.5). Although the conductivity of fast and slow 
Ca?* channels is unaltered, there is a fall in intracellular calcium concentration, 
[Ca?*] = 0.034 uM. As a result, the fasciculus remains relaxed throughout. 

Treatment of the electrically stimulated detrusor with charybdotoxin (CTX), a 
selective maxi Ca?*-activated K* channel blocker, exerts an excitatory effect on its 
myoelectrical activity. Charybdotoxin does not affect the resting membrane 
potential but abolishes slow waves and AHPs. There is an increase in the ampli- 
tude, V = 50 + 55 mV and the frequency, v = 8 Hz, of action potentials. The 
level of free intracellular calcium rises to a maximum 0.48 uM and the myofiber 
generates tonic-type contractions, max7? = 1.5 mN/cm (Fig. 6.6). Similar elec- 
tromechanical effects are observed after application of iberiotoxin. 

Concomitant addition of nifedipine disrupts the high frequency bursting activity 
in detrusor. The cell membrane becomes hyperpolarized V' = —52 mV and the 
fasciculus generates steady slow waves of amplitude V = 20 mV and a constant 
frequency, v = 0.85 Hz. Smooth muscle produces bursts of 7—9 action potentials, 
V=50 + 55 mV, at a frequency, v = 2.5 Hz. Free intracellular calcium of 
0.17 uM triggers the active force development of intensity max7^ = 0.9 mN/cm. 
The contractility pattern also changes from a tonic to phasic-type. 
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Fig. 6.3 Dose-dependent effects of mifebradil on electromechanical activity and the dynamics of 
selective ion currents of the detrusor fasciculus 


6.3.4 K* Channel Agonist/Antagonist 


Lemakalim (LEM), a selective K*-channel agonist, has a strong hyperpolarizing 
effect on smooth muscle cells of the bladder. The resting membrane potential 
reduces by 26 mV (V' = —70 mV). LEM eliminates completely slow waves and 
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Fig. 6.4 Effects of nifedipine and verapamil on electromechanical activity and the dynamics of 
selective ion currents of the detrusor fasciculus 


action potentials. The concentration of free [Ca?*] ~ 0 and the detrusor fails to 
contract (Fig. 6.7). 

Addition of tetraethylammonium chloride (TEA), a nonselective voltage-gated 
K*-channel antagonist, increases the amplitude, V = 55 + 60 mV, and the fre- 
quency, v = 20 Hz, of action potentials (Fig. 6.8). These effects are dose-dependent. 
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There is an increase in the influx of calcium ions inside the cell, [Ca?*] = 0.47 uM, 
and a rise in the tension of contractions, T* = 1.46 mN/cm. 

Detrusor continues to generate spontaneously spikes of the amplitude V = 28 
+ 34 mV and the frequency, v = 3.5 Hz, in the presence of TEA even after 
withdrawal of external electrical stimuli. The smooth muscle fasciculus evokes 
long-lasting bursts of regular phasic contractions of a maximum 1.38 mN/cm. 

4-aminopyridine (4-AP), a selective K^-channel antagonist, depolarizes the 
membrane, V" = —42 mV, but does not affect slow wave activity. There is an 
increase in the amplitude and frequency of spikes, V — 47 mV, v — 7.1 Hz, 
without significant changes in the dynamics of their generation (Fig. 6.9). The 
intensity of Tk is reduced by 25 %, max/k — 5.6 nA, and the effect is dose- 
dependent. 

The maximum influx of extracellular calcium is not altered. Noticeable changes 
are seen in the amplitude of oscillations in free Car, The detrusor smooth muscle 
contracts at a relatively constant strength of 1.44 mN/cm. 
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Fig. 6.5 Dose-dependent effect of forskolin on electromechanical activity of the detrusor 
fasciculus 


6.3 Effects of Pharmacological and Extracellular Ion Changes 99 


moo CTX 
EN Nifedipine 

9.0 
_~ -14.0 
> 
E 
= -37.0 
> 

-60.0 LLLLLI LEPLLLI 
0.0 25 50 75 100 


1.5 
z m 
x 143 
z 
ES = 
e) 
v 0.7 
Ó 3 
2 LLELELEEELEEEEEEETELLTELLL ILLI 0.3 
0.0 25 50 75 100 
Time (s) 


Fig. 6.6 Selective and conjoint effects of charybdotoxin and nifedipine on electromechanical 
activity of the fasciculus 
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Fig. 6.7 Electromechanical response of the detrusor to application of lemakalim 


6.3.5 Ca?*-ATPase Inhibitors 


Two agents have been described as selective inhibitors of the Ca?+-ATPase: 
cyclopiazonic acid (CPA), a mycotoxin from Aspergillus and Penicillium, and 
thapsigargin, a natural compound. CPA easily passes into the cytoplasm through 
the plasma membrane and reduces Ca?*-ATPase activity while the effects of 
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Fig. 6.8 Electromechanical response of the detrusor to tetraethylammonium chloride 


thapsigargin are more complicated. At concentrations «1 uM, it prevents the 
filling of the Ca?* stores, while at higher concentrations it interacts with voltage- 
dependent L-type Ca? channels in the plasma membrane. 

The overall effect of CPA, an increase in the magnitude of membrane ionic 
currents and the electromechanical activity of detrusor, is time-dependent 
(Fig. 6.10). Thus, the inward slow and fast calcium currents reach J$, = 2.0 nA 
and IL, = 0.024 nA, respectively, while the dynamics of the transient outward Tk 
and Iu currents is not affected. 

CPA depolarizes the detrusor muscle membrane by 8 mV and abolishes slow 
waves. The detrusor muscle shifts to a new stable excitable state in which high 
amplitude, 52 mV, and high frequency, 30 Hz, action potentials are generated. 
The amount of free intracellular calcium is slightly reduced, however, is reduced, 
max [Ca?*] = 0.51 uM and as a result, the intensity of phasic contractions falls to 
T* = 1.2 mN/cm. 
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Fig. 6.9 Dose-dependent effects of 4-aminopyridine on electromechanical activity of the 
detrusor fasciculus 


Thapsigargin causes an increase, compared to the norm, in the frequency of là, 
and IE currents and their amplitudes. The dynamics of the currents changes to 
bursting. The amplitude of the influx of Ca?* through the L-type Ca?* channel is 
0.18 nA and that through the T-type Ca?* channel is 0.04 nA. Also, there is 
amplification of the outward Tę and Ic, y currents. 

Thapsigargin has little effect on muscle resting membrane potential but it 
increases the frequency, v — 0.57 Hz, of self-oscillatory activity. Continuous 
bursts of action potentials, amplitude 38 — 55 (mV) and frequency, — 19 Hz, are 
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Fig. 6.10 Effects of cyclopiazonic acid and thapsigargin on electromechanical activity and the 
dynamics of selective ion currents of the detrusor fasciculus 


generated. There is a rise in the concentration of intracellular calcium, max 


ca? 
ulus 


*] = 0.48 uM, with greater amplitude in its variation, 0.21 uM. The fascic- 
produces high frequency phasic contractions, max7? — 1.49 mN/cm. 


The response of the detrusor myofiber pretreated with CPA is sensitive to a 


conc 


omitant addition of nifedipine and verapamil. Nifedipine has shown to have a 
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Fig. 6.11 Selective and conjoint effects of cyclopiazonic acid, nifedipine, and verapamil on 
electromechanical activity of the detrusor fasciculus 


detrimental effect on 73, and IL, currents, oscillatory and spiking activity of the 
urinary smooth muscle. The frequency of action potentials reduces significantly, 
v = 0.03 Hz while their amplitude remains unchanged, V = 52 mV (Fig. 6.11). 
The drug has a positive effect on the amplitude and the frequency of slow waves. A 
decrease in the amplitude of the inward calcium flux results in a decrease in the 
level of free cytosolic calcium, [Ca2*] — 0.28 uM. The myofiber sustains its 
phasic contractility and generates the active force of 0.9 mN/cm. 

Subsequent washout of nifedipine and application of verapamil completely 
abolishes electromechanical responses in detrusor. The fasciculus becomes hy- 
perpolarized, V = 58 mV, and remains relaxed, 7? = 0. 

The detrusor smooth muscle cell fails to respond to CPA in the environment 
free of extracellular calcium, [Ca?*] = 0. It becomes hyperpolarized, V = 52 mV, 
and no active forces are produced (Fig. 6.12). An increase in the concentration of 
pu has an excitatory effect on the fasciculus with the development of strong 
slow calcium I = — 1.1 + —32 (nA) and fast calcium I = — 0.0014 + 
—0.045 (nA) ion currents, and the generation of high frequency and amplitude 
spikes, maxV — 60 mV. The amount of free intracellular calcium also rises, max 
[Ca?*] = 0.37 uM, and as a result the active force of contraction of the intensity 
1.3 mN/cm develops. 
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Fig. 6.12 Conjoint effects of cyclopiazonic acid and excess extracellular Ca?* and Kt ions on 
electromechanical activity and the dynamics of selective ion currents of the detrusor fasciculus 
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Concomitant application of CPA and elevation of the concentration of extra- 
cellular potassium ions leads to stable depolarization of the detrusor by ~21 mV 


(V = —29 mV) and development of a long-lasting tonic-type contraction, 
T* = 1.8 mN/cm. 


Chapter 7 
The Intrinsic Regulatory Pathways 


There is something fascinating about science. One gets such 
wholesale returns of conjecture out of such a trifling investment 


of fact. 
M. Twain. 


7.1 Models of Electrical Activity of Neurons 


Electrical activity in neurons is ubiquitous and is essential to maintain vital 
functions of the body. It is manifested by a variety of patterns of discharges (action 
potentials, spikes) and the specificity of transduction mechanisms that enable the 
system to integrate and to coordinate the overall behaviors in space and time. 

At the level of mathematical modeling of neuronal activity two types of 
complexity must be dealt with: the interplay of ion channel dynamics that underlie 
the excitability, and the neuronal morphology that allows neurons to communicate 
signals among them. Despite of significant advancements in simulation and a large 
number of proposed mathematical models of excitable media the only biologically 
plausible model remains the Hodgkin and Huxley (H-H) model (1952). It contains 
a current balance equation that satisfies Kirchoff’s law 


dV ANM. 
"ar - 2 N= Vi) + Tappi (t), va) 


and relaxation equations for the n ionic conductances 


dx; E Xi (V) — Xi 


= , i-ln, 
df T(V) (7.2) 
dy; L Yioo(V) —»i AT 
4^ wy 


Here C is the membrane capacitance, V is the membrane potential, V; is the 
Nernst potential for the ith ion, g; is the maximal conductance of the channel for 
the different ion-selective channels, x? (xis), Y (vis) are the actual (steady) state 
activation and inactivation variables, respectively, Tx, and c, are the relaxation 
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time constants, lappi is the applied current. The nonlinear system (7.1), (7.2) 
accurately reproduces the rich dynamic behavior in various excitable tissues 
ranging from tonic bursting to chaos, and from stationary to traveling wave phe- 
nomena. Such versatility and reliability is a result of mathematical formulation 
that incorporates explicitly/implicitly intrinsic properties of the cell membrane. 
Dynamic systems analysis of (7.1), (7.2) has provided insights into how cell 
activity is shaped by individual parameters. Thus, it has been demonstrated that a 
qualitative change from large amplitude, stable periodic behavior to small 
amplitude and unstable periodic oscillations—a subcritical Hopf bifurcation—can 
be achieved by an increase in [a2]. a reduction in membrane maximal K* 
conductance, or by a shift of the Nernst potential for Kt in the depolarizing 
direction (Holden and Yoda 1981). 

Several other models, which are simplifications of the H-H model to mimic 
membrane potential and ion current dynamics (Fitzhugh 1955; Kepler and Marder 
1993; Wilson 1999; Gerstner and Kistler 2002; Izhikevich 2003) have demon- 
strated their robustness in describing complex electrical events in large populations 
of neurons and neural fields. However, their applications in "electronic pharma- 
cological" studies are severely limited because of constructive biological insuffi- 
ciencies (Miura 2002). 

A distinctive feature of neurons is their morphology. Much of the research 
recently has been directed at understanding the impact of anatomical architecture 
on signal integration within a neuron and neuronal ensembles. Originally, a single, 
electrically equivalent cable model was used to simulate the generation and 
propagation of voltage waves and current changes in the axon and dendrites 
(Hodgkin and Huxley 1952; Schierwagen 2009) 


dV d; ev n 
dt 2R. = my (V — Vi) +h 7. 
dr 2R, 0€ ds Yi (V — Vi) + lappi (t), (1.3) 


where Rm is the membrane resistance, d; is the cross-sectional diameter of the 
nerve fiber, x is the Lagrange coordinate (0 € x € L), and the meaning of other 
parameters and constants is as described above. However, it could not answer 
questions related to the neuronal branching and synaptic transmission. As an 
extension, compartmental models have been developed to explore the role of 
evolving dendritic morphology and to embrace biological mechanisms at molec- 
ular and cellular levels that support the temporal and spatial stability of signaling 
in diverse neuronal topologies (Graham and van Ooyen 2006). 

Neurohormonal modulation and electro-chemo-mechanical coupling in 
detrusor involve a cascade of electrical and chemical processes including gen- 
eration and propagation of the electrical signal along axons and dendrites, 
synthesis, storage, stimulation, release, diffusion, and binding of various sub- 
strates to specific receptors with activation of intracellular second messenger 
systems and the generation of a variety of physiological responses. Qualitative 
analysis and quantitative evaluation of each and every step experimentally is 
very difficult, and sometimes practically impossible. Therefore, different classes 
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of mathematical models of neurohormonal modulation and synaptic neuro- 
transmission, ranging from the most comprehensive "integrated" —microphysi- 
ological—to “reductionist” —deterministic—have been proposed to study 
intricacies of the processes of neuroendocrine regulations. With a microphysio- 
logical approach, attempts to reproduce reality in great detail lead to mathe- 
matically challenging and computationally demanding tasks. In contrast, 
deterministic models aim to capture accurately phenomenological behavior of the 
system. They not only provide macroscopic explanation of complex biophysical 
processes but are general enough to offer a coherent description of essential 
biochemical reactions within the unified framework. These models are inherently 
flexible and can accommodate spatiotemporal and structural interactions into a 
tractable representation. 


7.2 A Model of Neuronal Assembly in the Bladder 


Bearing in mind the paucity of qualitative and quantitative data about morpho- 
logical arrangements, electrical and pharmacological properties, the receptive field 
distribution and projections in the major pelvic and intramural ganglia of the 
human urinary bladder, consider a hypothetical neuromuscular module as shown in 
Fig. 7.1. 

It includes the sensory fiber, the postganglionic parasympathetic motor neuron, 
the postganglionic sympathetic nerve fiber, and the detrusor muscle. They are 
linked together by multiple synapses which: (1) provide the major postganglionic 
sympathetic and parasympathetic inflows to the detrusor, (2) pass on afferent and 
regulatory signals from the stretch receptors, the postganglionic pelvic and the 
hypogastric nerve fibers to the motor neuron, and (3) inhibit the activity of the 
postganglionic pelvic nerve fiber presynaptically. The sequence of events that it 
describes includes: 


1. mechanical stretch of the mucosa triggers a depolarization in the free nerve 
endings of the sensory afferent fiber; the excitatory signal passes along the 
unmyelinated nerve fiber to the somatic cholinergic synapse on the motor 
neuron; 

2. the motor neuron receives postganglionic parasympathetic (excitatory) and 
sympathetic (inhibitory) inputs from the major pelvic ganglion; it is assumed 
that the strength of these inputs is known a priori; the parasympathetic axonal 
terminal is under the inhibitory control of the sympathetic outflow; 

3. discharge of the motor neuron and generation of the wave of depolarization that 
travels along the unmyelinated axon to reach the detrusor muscle; additionally, 
the detrusor receives a direct inhibitory signal from the postganglionic hypo- 
gastric nerve. 


Although it has been confirmed experimentally that multiple transmitters are 
colocalized and coreleased in the process of neurotransmission at different 
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synapses at different levels, we shall focus our attention in the future on the 

analysis of the two major substrates only, i.e., acetylcholine and adrenaline. 
Assume that mechanoreceptors convert mechanical stretch stimuli, e(t), to the 

dendritic receptor potential, Va, via activation of stretch-activated Na* ion channels 


dV, & - s 
Cy = = —(Ina + Fx + Jci) + (Vi — Va) /Rm (7.4) 


dV. 
Tm a7 Vi — k(V; — Va) /Ra 


where C, is the dendritic membrane capacitance, V; is the nerve ending potential, 
Tm, k are the membrane time and numerical constants, respectively, Tua, Ik, la are 
the ion current 


Ína = gu h(Va — Vya), Ik = gk (Va — Vk), la = Ba(Va — Va). — (7.5) 


Here @na, £k, £ci are the maximal conductances of Na*, K* and Cl” channels, 
m, n, h are probabilities of opening of these channels, Vya, Vg, Vy are the reversal 
potentials of the respective currents. The activation and deactivation of the 
channels is described by 

dy * * cm n], 
Dc 0 = 9") — By) "=A, fh) (7.6) 
where the rates of switching the channels from a closed to an open state, %,, and 


reverse, f,, are 


&, = 0.221 exp(e(t) +0.01Va) _ — f, = 4.5exp(—Va/18) 
&, = 0.048exp(—V4/36), — fj, = 40.12/(1 + exp(3.4 — 0.2Va) (7.7) 
&, = 0.33 exp(1.1 — 0.1Va), Ba = 0.185 exp( —Va/18). 


The dynamics of the soma action potential at the motor neuron, V,m, is a result 
of the integrated function of five ion channels including voltage-dependent 
Ca?* , Ca?*-activated K*, and the H-H type Na*, K* and Cl” channels 


Fig. 7.1 Scheme of the one- Sensory Motor 
dimensional neuromuscular neuron neuron 
module used in numerical 
experiments to simulate 


Postganglionic 


responses by the detrusor. pelvic nerve 
Here ———» is the Sensory fiber 
excitatory, and ———- is the receptors 


inhibitory synapse 


Detrusor 
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Cam m — (Ic n Ica_k F Ina F Îk + la) + Text; (7.8) 


where Cpm is the membrane capacitance, i...) are ionic currents carried through 
different ion channels, and 7.4, is the external membrane current (input). The 
equations for currents are 


T Baz ^ 

I a — Vam Vi a 

B 1+¢[Ca?*] ( ca) 

5 BCa-K foc [Ca?*] 7 

Ica-K = Vim — Vea— 

A E E (1.9) 


Ina = anal A( Vam = Vya)8c1(Vam = Vk) 
Ix = êk (Vam — Vk) 
la = &c1(Vam — Va). 


Here Vc are the reversal potentials for the respective currents, 2(.. are the 


maxi-mal conductances, ( is the parameter for calcium inhibition of Ca?* chan- 
nels, z, 94, , Mñ, fi, h are dynamic variables of the ion channels given by 


z= (ts — 2)/ ts Gi = (foo — Â) /Ta 
Gh — Ag (hoo — h)/t, fh = fus (Vas) (7.10) 
Poo = (1 + exp 0.45(Vam + 57)) |. 
Here 
Joo = Byoo Ty ($ — i, h, it) 
Ty = 1/(&yos t Pu) 
; 0.12(Vam + 27) A 
— ; ass = 4A exp(Vam + 47)/25 
Ae = T oxp Vaa 425/89 P cep Von 147) 
" 1 
à, = 0.09 exp(—(Vam + 47)/27), Bi. = (7.11) 
OG oc exp( ( =F )/ ), Pi. 1 + exp(—(Vam Em 33)/8) 
P 0.012(Vam + 12) d 
foo = ; Weg = 0.12 — Vum 26)/87 
As =T expl (Van + 13/8) ^ expUs Vase par 
1 
Zoo 


~ 1+ exp(—0.15(Vam + 42)) 
Dynamics of the free intracellular calcium yield 


d[Ca?*] 0.234 x 10-4z(Vam — Vea) 


24 
d LEO" 0.003[Ca?*]. (7.12) 
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The model of the detrusor fasciculus, initial and boundary conditions are given 
by Eqs. (6.4)-(6.15), (6.16)-(6.18) (see Chap. 6 for details). It is assumed that the 
motor neuron provides an ionotropic input to the L-type Ca?™ channel via a 
cholinergic 43 synapse. The combined system of equations comprises the math- 
ematical formulation of the neural control of the functional unit of the urinary 
bladder. 


7.3 Response of the Detrusor to Stretching 


Distension of the urinary bladder with urine, ~300 ml, causes the stretch of 
mechanoreceptors in the mucosa and evokes a receptor action potential Vy = 
82.5 mV of a duration 30 s at the free nerve afferent sensory fibers. A characteristic 
feature of the potential is the long lasting “step” with the constant amplitude of 
approximately 45 mV. The level of depolarization is sufficient to develop and to 
sustain the propagation of the wave of excitation, Vf = 69 mV along the unmy- 
elinated dendrites. It reaches the dendrosomatic synapse on the motor neuron and 
activates electrochemical signal transduction mechanisms with the generation of 
the fast excitatory postsynaptic potential (EPSP), Vsyn = 87.1 mV (The details of 
the multistep process are discussed in the following chapter). 

Depolarization of the membrane of the motor neuron triggers ion currents 
of intensities max Ic, = —1.2, Ina = —2, Ica_k = 1.6, Jk = 1.2 — 1.8, and Îcı = 
1.5(nA), respectively. Their strength and dynamics are defined by the activation and 
deactivation properties. As a result the soma of the neuron produces tonic bursts of 
high amplitude spikes, Vam = 72 mV, at a frequency ~3 Hz. They propagate along 
the axon toward the neuromuscular cholinergic synapse on the L-type Ca?* channel 
on the detrusor membrane. Their activation causes an influx of extracellular calcium 
inside the cell and, as a result, the generation of high frequency and amplitude 
oscillations of the potential V = 50-60 mV, v = 4-6 Hz. A rise in intracellular 
calcium [Car] = 0.48 uM and the activation of contractile proteins leads to the 
development of force, maxT — 1.7 mN/cm (Figs. 7.2, 7.3). 

Electrical activity in the neuronal pathway ceases upon voiding and the detrusor 
returns to its resting state. 

The response of the smooth muscle fasciculus to simultaneous deformation and 
spontaneous discharges of action potentials does not differ significantly from that 
described in previous sections. The detrusor responds with phasic contractions of 
varying amplitude 0.4—0.6 mN/cm and a frequency 0.4 Hz (Fig. 7.4). 


7.3 Response of the Detrusor to Stretching 


Fig. 7.2 Normal 
electromechanical activity of 
the ganglion-detrusor in 
response to stretching. Traces 
from top to bottom indicate: 
depolarization wave 
dynamics on the soma of the 
motor neuron, the detrusor, 
intracellular calcium changes, 
and total force 


Fig. 7.3 The dynamics of 
ion currents and the action 
potential on the soma of the 
motor neuron 
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7.4 Pharmacology of the Neuronal Pathway 


7.4.1 Effect of Iberiotoxin and Charybdotoxin 


Iberiotoxin and charybdotoxin are specific inhibitors of high conductance Ca?*- 
activated K* channels. During the numerical simulation the mechanoreceptor field 
has been excited by a single stretch ¢(2.0) = 0.75. Application of CTX and ibe- 
riotoxin does not affect the activity of free sensory nerve fibers. However, there are 
significant changes in the pattern of firing of the motor neuron. The soma of the 
neuron responds to the electrical signals from the free nerve endings by generating 
ion currents of intensities: Ica = 6.9, tes = 0.016, Ix = 0.15 and iq = 0.8(nA), 
respectively. Action potentials of amplitudes 94 mV at a frequency of 3.1 Hz 
occur as continuous discharges (Fig. 7.5). 

There is a concomitant increase in the frequency and amplitude of action 
potentials, without effect on their configuration, produced by the detrusor. Cyto- 
solic calcium transients result in strong and prolonged contractions of the smooth 
muscle fasciculus. 
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Fig. 7.5 Effect of Ca?*-K 
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7.4.2 Effect of w-CgTX 


c-Conotoxin is a selective, long-lasting blocker of N-type Ca?* channels. It 
completely blocks electrical activity at the level of the soma of the effector neuron 
which remains hyperpolarized throughout. No action potentials are produced. As a 
result, no electromechanical activity is detected in the detrusor. 


7.4.3 Effect of TTX 


The mechanism of action of TTX is a decrease in the permeability of Na* and K* 
channels. It is also suggested that in neuronal cells treated with TTX, sodium ions 
can enter the cell via voltage-dependent Ca?* channels. Therefore, it is of par- 
ticular interest to analyze the influence of the intensity of mechanical stimuli on 
the electrical response at the level of free nerve endings. 

The results of numerical simulations have shown that treatment of the detrusor 
fasciculus with TTX does not affect electrical activity of mechanoreceptors. 
Generated dendritic potentials of amplitude 84 mV initiate the traveling wave of 
excitation that propagates along the afferent fibers toward the dorsal root ganglion. 

In contrast, the motor neuron appears to be very sensitive to the toxin. TTX 
abolishes its electrical activity and, as a result, the detrusor remains relaxed 
throughout. 


7.4.4 Effect of Purinoceptor Agonists 


There is experimental evidence that two distinct adenosine subtype receptors, A; 
and A3, are present on intramural neurons. An A, receptor agonist—2-chloro-N°- 
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cyclopentyladenosine (CCPA) exerts its effect by increasing the permeability of 
voltage-gated K* and BK, channels. Treatment of the functional unit with CCPA 
results in a continuous bursting of spikes by the motor neuron and in intense phasic 
contractions of the detrusor. 

Application of selective A; and A; receptor antagonists-NF449 and A-317491, 
respectively, decreases the input resistance, and increases the conductance of BKc, 
channels. The compounds have no effect on electrical activity of the mechanore- 
ceptors. However, they hyperpolarized the motor neuron, which remains unexcited 
throughout, and no mechanical activity by the smooth muscle fasciculus is recorded. 


7.4.5 Effect of Protein Kinase C Activator 


A proposed mechanism of action of phorbol 12,13-dibutyrate (PDBu) on neurons is 
a simultaneous increase in the permeability of voltage-dependent Ca?* and BKca 
channels. Treatment of the functional unit of the bladder with PDBu at low con- 
centration results in a regular bursting of the soma of the motor neuron. APs of 
amplitude 82 mV and a frequency of 17 Hz are produced. An increase in the con- 
centration of PDBu leads to the shift in the pattern of discharges to a beating mode at 
the maximal frequency 20 Hz while the amplitude of APs remains unchanged. 
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Chapter 8 
The Synaptic Transmission 


Science does not know its debt to imagination. 
R. W. Emerson. 


8.1 A Model of the Synapse 


Let a synapse be an open three-compartmental system. Compartment 1 comprises 
presynaptic elements where synthesis and storage of a neurotransmitter or a hor- 
mone occurs. For example, in case of cholinergic or adrenergic synapses it cor- 
responds morphologically to a nerve terminal of the unmyelinated axon. 

Ligands, i.e., acetylcholine, adrenaline are released upon neural stimulation by 
exocytosis to the synaptic cleft and bloodstream. The common sequence of events 
involved in the dynamics of their transduction includes: 


1. depolarization of the nerve terminal or cell membrane; 

2. influx of extracellular calcium through voltage-gated Ca?" channels; 
3. binding of free cytosolic Ca?*to transmitter-containing vesicles; 

4. release of vesicular/granular stored ligand, L,, into the synaptic cleft. 


Propagation of the wave of depolarization, Vf, along the nerve fiber and the 
nerve terminal is described by the modified Hodgkin-Huxley system of equations 
(Miftahof et al. 2009) 


oV 1 0 ovt l 

f Ë f f 

AE d? d Gg. 44 

m Qr = Raa (4 (a) e) Una + Ik +e) 
ad, I <a X I5 — I3 

ds (a) = S S S 
2d; a> D — Lẹ, t» 0, « € (0, D). 


(8.1) 


Here C is the specific capacitance of the nerve fiber, R! is the membrane 


resistance, d; is the cross sectional diameter of the terminal, « is the Lagrange 
coordinate, and L*, Lj are the length of the axon and the terminal, respectively. The 
above ion currents are defined by 


= ghamphy(V! — Va) 
Ik = gnghj(Vf — Vk) (8.2) 
A = gal vf - Vi); 


R. N. Miftahof and H. G. Nam, Biomechanics of the Human Urinary Bladder, 117 
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where gui gh, an are the maximal conductances for Na‘, K*, and C17 currents, 


respectively, Vis Vie Vy are the equilibrium potentials for the respective ion 
currents, and my,ny,hry are the state variables that are calculated from 


dy/dt —a,(1—y)—B,y y= (mp, ny, hy). (8.3) 


The activation, «,, and deactivation, f,, parameters satisfy the following 
empirical relations 


| QAT(2.5 - VP) 7 j 
“mf exp(2.5 — 0.1V/)' Bms = 4 Texp(-V /18) 
ong = 0.07Texp(—0.05V/^), Bhs =T/(1 + exp(3 — 0.1 V/)) (8.4) 
0.17 (10 — V^) 


nf — z ; n — 0.125 T —0.125 V ). 
"f = BD -0dV)-1' Pas expt ) 


Here Tis temperature. 
The cytosolic calcium turnover is given by 


d[Ca?* 
1 — ge WOC] — lca, (8.5) 
where gs is the conductivity of the voltage-gated Ca?^* channel at the synaptic 


end, ky is the intracellular buffer system constant. 
Finally, the release of a stored fraction of the ligand, L,, is described by the 
state diagram (Fig. 8.1) 


d[X;] 


= —ko[Xi][X2], (8.6) 


Here X, := [Ca?*], X» :— [L,] are concentrations of the cytosolic calcium and 
vesicular-stored ligand, respectively, ko is the rate constant of association of 
Ca?* with calcium-dependent centers on the vesicles, kı is the diffusion constant. 

At t = 0 the nerve axon and the synapse are assumed to be at the resting state 
and the concentrations of reacting components are known 


my(0) = mro, hj(0) = hyo, nj(0) = nro; 


Vf (0, a) = 0, L,(0) = Lyo. pug 


Fig. 8.1 The state diagram k,Ca* 
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The synapse is excited at the free end (x = 0) by the electric impulse of an 
0 
amplitude Vf and a duration 1^, and the presynaptic terminal end (x = L) remains 


unexcited throughout. 


0 
f d 
V(0,) 2 4 V» O«t«t Wit) =0 (8.8) 
0, 27 


Compartment 3—the postsynaptic membrane—contains a complex network of 
surface membrane proteins which modulate the transmission. It involves the fol- 
lowing processes: 


1. binding of free ligand in the cleft, L,, to the G protein-coupled receptor, R, and 
its conformational change, L, - R; 

2. active configuration of the receptor, R*, (the Le - R* reactive complex is able to 
produce a biological effect, i.e., the postsynaptic potential, Vsyn ); 

3. binding of Le - R* to G protein, the formation of the Le - R* - G —complex and 
the initiation of guanosine diphosphate/guanosine triphosphate (GDP/GTP) 
exchange; 

4. the dissociation of G protein « and fy subunits with the subsequent release of 
Gacı protein that interacts with downstream effector pathways; 

5. enzymatic, E, clearance of the excess of L, in the synaptic cleft through the 
formation of intermediate complexes, Le - E and ES, and the final metabolite, S. 


They are described by the following state diagram (Fig. 8.2). 


Fig. 8.2 The state diagram 
of a neurotransmitter 
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Here k+; are the forward (+i) and backward (—i) rate constants of chemical 
reactions, C;(j = 1,4) are the intermediate complexes, and the meaning of other 
parameters as described above. 

Assuming that: 


1. the distribution of reactive substrates is uniform throughout and no chemical 
gradients exist; 

2. the total enzyme concentration does not change over time, Ey = constant; 

3. the total ligand concentration is much larger that the total enzyme 

concentration; 

4. no product is present at the beginning of the reaction, and 

5. the maximum rate of chemical reaction occurs when the enzyme is saturated, 
i.e., all enzyme molecules are tied up with a substrate, 


then all reactions are of the first order and satisfy the Michaelis-Menten kinetics. 
Hence, the system of equations for the ligand conversion is given by 


dX/dr = AX(t), (8.9) 


where the vector X(t) = (xj) G — 1,20) has the components 


Xj = [C277], X; = [L], X; = [L Cat], X4 = [L], 
Xs = [R], X6 = [R*], X; = [Le R'], X; = [Le R], 
* * 
Xo = |Le- R* -G], Xio=[R*-G], Xu = [A - E], Xn = [G], 
Xi5 = [Ci]; Xu = [Co], X15 = [C3]; X16 = [C4], 
A 
Xy = [G], Xis = [Gaci], X19 = [ES], X29 = [S]. 
The matrix A(aj) (i,j = 1,20) has the nonzero elements 
an = ge V (O) [Cag], an = —ko[Xi]. 413 = ko 
an = —ko|Xi], az = k-o 
az = ko[Xi], a33 = —(k-o + kı), 
ag = ki, 444 = ka ([Eo] — [Xs]) + ko [Xs] + k-s[X6] + k-11 [X10], 
ası = kg, dag = k6, as = kii, 
daa = ko, 
ass = —(k.s + ks) [Xa], ass = ks[Xi2], asg = k_6, 
des = k-s, aso = —k- [Xi] — (ks + Kioio[Xis])[Xi2] 
ası = kg[Xi2], 6,10 = kio, 6,16 = ki7, 
aw = k_s|Xa], an = —(k-7 + ks + ko)kig[X12][Xis], 
an = ky, ok12 [X40], 4744 = ki, 
ass = ke|X4], 1X12]; agg = —k_¢ — ky, 
49; = kok o[X12][Xis], doy = —k_g — ku — kı2[X40], 4549 = k-n [X4], 
4106 = k-io[Xi2], aio = kr, 41010 = — (kio + K-11) [Xa] — Kis [Xao] 
ana = b([Xa] — [Xs]), ania = —(k2 + ks), anas = ks, 
ai26 = k-s + k-s[Xa], ang — kn ang =k», 41240 = kio 
41242 (ks + K-iofio[Xia]) [Xo] — (k-7 + ks + KokioXis]) [Xs], 
4139 = kin[Xao], 43,13 = —kı4 [X39], 1314 = k-14, 
41443 = Kia [Xso], aja = —(k14 + kie), 
415,10 = Kis[Xao], 1545 = —kis[X30], 415,16 = Kis, 
41645 = Kis[Xso], 41616 = — (kis + k-15), 
41744 = kis; 41746 = kı7, 41717 = —kis; 
ag = kis, 41818 = —kio, 
aiga = ks, 419,49 = —(k-3 + ka), 41920 = k-4([Xai] — [Xs]), 


40,19 = ka, 42920 = —k_4([Xai] — [Xs]) 
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Here [X39], [X40], [X4i1]are given concentrations of GTP, GDP, and E enzymes, 
respectively. 

The generation of the excitatory postsynaptic potential, V.y,, is given by 
(Miftahof et al. 2009) 


dV, 
dt 


Vs n 
+ Van (7 Q[Xo] + Ry!) = > (8.10) 


Cp 
where C, is the capacitance of the postsynaptic membrane, A, is the resistance of 
the synaptic structures, Q is the empirical constant, Vsyno is the resting postsyn- 
aptic potential. 

Given the concentrations of reacting components and the state of the synapse 


X(0) = Xo, Ven —0 (8.11) 


Equations (8.1, 8.11) provide mathematical formulation of the dynamics of the 
common pathway of neurotransmission. 

The quantitative assessment of the velocities of reactions shows that the rates, 
ki», kı3, of exchange of G protein for GDP at L.. R* -G and R*-G sites, 
respectively, are significantly smaller compared to the rates of other reactions. 
Hence, the system (Figs. 8.2 and 8.3) can be viewed as a combination of rapid 
equilibrium segments interconnected through slow, rate-limiting steps. The char- 
acteristic feature of such system is that it attains the steady state when the rapid 
segment has already reached quasi-equilibrium. This fact allows us to simplify 
(8.9) as follows: let Ji and h be the fractional concentration factors of the rapid 
segment “product” -substrates 


f-IXg/Y. h = [X10]/] (8.12) 
where Y is 


[Y] = [Xo] + [X10] + [X12]. 


Define the association constant as 
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[X6][X12] x [Xa] [X10] " [X7][X10] T. [Xa][Xs] 


K = , = E 3 = 3 = , 
^w ^ pm qe t gx 
after simple algebra, for fis we have 
; _ [Xa][Xe]/KiK2 _ [X4] [Xs] 
f- - 
Coq BE E] Ki Ko + [Xo] Dra] + K2) 
142 1 
(8.13) 
pee LC NY K2|X6] 
2 1 q Ex] | Xel KioKi, + [X6]([X4] + Ka) ` 
KiK Ki 


Then, concentration distribution equations for the main reactants can be 
obtained in the form (for details see King and Altman 1956; Cha 1968) 


(D)(Y]/[Y ]) = kiakiskiski7kiski9[X30]” 


= kiokiskizkigkiof | [X30](k-14 + Ki) 


(8.14) 


(231/171) 
(y (pat) = kokukiskiskiskiof, [Xs] 
( ) = ki kakiokigkiof [X39] (Kk 15 + ki7) 


(D)([Xig]/[Y]) = kukiakiskiekiskiofo[X30] 


(D) (ps/(Y']) = klaskiskisfa eol" (ofi + kufa) 
where 


[Y ] = [Y] + [X13] + [Xia] + [5] + [X16] + [X7]; 
and D is the sum of all the values on the right side of Eqs. (7.13). Note that Eqs. 
(8.14) are algebraic equations. The initial velocity equation for Y complexes for- 


mation is given by 


d[Y]/dt = kigki9[Xq7]. (8.15) 


Finally, substituting X,7 from Eq. (8.14) and making use of Eqs. (8.13) we get 
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[X4][X6] BA 
2 | ho TER 
Lay | kyakiskiokizkigk19 [X39] eee | 
[Y] dt kyakiskioki7kiski9 [X39] + ' 
[Xa] [Xo] 
kiokiskiskiskio[Xso] (k-14 + kie) P e E 
Kj Ko Ki 
[Xa] [Xo] 
+kkiakiski7kigkio[X30]- S N al 
KiK Ki 
kia [X4] [Xo] +k [Xs] 
+kakiskioki7[X30]” cara E 
1+ KıK2 tK 
p] 
2 K 
ykiakiskigkigkio[Xao] TEM 
Ki, ^ Ki 
[Xe] 
K 
+ki1kışskiokiskio [X39] (k-15 + kı7) rie] 
Kk ^ Ki 
or, after some algebraic rearrangements, in the form 
1 d[Y] E kiakiskikizkigkio[Xa9] [i2 [X4] + i1 K2] 
[Y] dt [X4] [X30] KiaKisKis [Ki6Kigkio + KizKigkio + ki2ki6|+ 


+[X39]kiakiski6K2[kizkigkio(1 + x) + kitki7 + kitkigkiol+ 
[X4] [Aiokiskizkigkio(k—14 + K16)] + KiikiakickigkioKo(k-15 + kı7) 
(8.16) 


Both of the mathematical formulations, i.e., the system of differential equations 
(8.9) or a simplified model given by (8.16), supplemented with initial and 
boundary conditions, provide a detailed description of electrochemical coupling at 
the myometrial synapse. In practice, the preference for either of them depends on 
the intended application and is entirely the researcher’s choice. 


8.2 cAMP-Dependent Pathway 


Stimulation of cholinergic 45 and adrenergic f-receptors leads to production of 
cAMP—the process that is controlled by a number of adenylate cyclase (AC) 
enzymes. They are dually regulated by a family of G proteins, forskolin, and other 
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class-specific substrates. A diversity and expression of ACs isoforms provides a 
mechanism for integrating positively or negatively the responses to various 
transmitters. Isoforms V and VI of ACs have high affinity to Gfy subunits of 
inhibitory proteins (Gao and Gilman 1991). Their functionality is affected by 
phosphorylation with PKC. Types V and VI are inhibited directly by low levels of 
Ca?* (Yoshimura and Cooper 1992; Premont et al. 1992; Cali et al. 1994). 
Adenylyl cyclase isoforms I, III, and VIII are shown to be upregulated by a 
calmodulin dependent protein kinase in response to the elevation in Cat, while 
the isoforms II, IV, and IX are stimulated by in the presence of Ga, subunits 
(Hanoune et al. 1997). Although the exact presence and distribution of the specific 
types of AC isoforms in the wall of the human urinary bladder are not known, it is 
obvious that the same stimulus may trigger different physiological responses in the 
detrusor depending not only on the type of receptors involved but also on the type 
of adenylyl cyclase to which they are coupled. 

The proposed structure of ACs consists of a short amino terminal region and 
two cytoplasmic domains. The latter are separated by two extremely hydrophobic 
domains which take the form of six transmembrane helices. The catalytic core of 
the enzyme consists of a pseudosymmetric heterodimer composed of two highly 
conserved portions of the cytoplasmic domains. It binds one molecule of Ga, 
which, in turn, catalyzes the conversion of ATP into cyclic adenosine mono- 
phosphate. cAMP influences a wide range of physiological effects including: (1) 
the increase in Ca?* channel conductance, (2) the activation of protein kinase A 
(PKA) enzymes, and (3) the cytokine production. 

Protein kinase A is a holoenzyme that consists of two regulatory and two 
catalytic subunits. Binding of cAMP to the two binding sites on the regulatory 
subunits causes the release of the catalytic subunits and the transfer of ATP 
terminal phosphates to myosin light chain kinase. The result is a decrease in the 
affinity of MLCK for the calcium-calmodulin complex and smooth muscle 
relaxation. Additionally, PKA may promote relaxation by inhibiting phospholipase 
C, intracellular Ca?* entry, and by activating BKca channels and calcium pumps. 
PKA activity is controlled entirely by cAMP. Under low levels of cAMP, it 
remains intact and catalytically inactive. 

The level of cAMP is regulated both by the activity of adenylyl cyclase and by 
phosphodiesterases (PDEs) that degrade it to 5'-AMP. At least 11 families of PDE 
isoenzymes are identified. Their hydrolytic activity is determined by the catalytic 
domain and conserved areas of amino acids specific to each family. Thus, PDE 4, 
7, and 8 hydrolyze cAMP; PDE 5, 6, and 9 hydrolyze cGMP, and the rest dem- 
onstrate dual specificity to cAMP and cGMP (Wheeler et al. 2005). In humans, 
PDE 1, 2, 3, and 5 enzymes specific for cGMP and PDE 4, highly specific for 
CAMP have been identified (Truss et al. 1995). Another mechanism that reduces 
the production of cAMP is the activation of Go; — 4, proteins which directly 
inhibits adenylyl cyclase through the MAPK signaling cascade. 

The state diagram of the cAMP-dependent pathway as described above is 
shown in Fig. 8.3. 
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The system of equations for the cAMP pathway activation is given by 
dX/dt = BX(t), (8.17) 
where X(t) = (x) (i = 21,26) has the components 


Xn = [AC], X72 = [ATP], X23 = [cAMP], 
X3, = [PKA], X25 = [PKA*], X26 = [5 — AMP], 


and the square matrix B(b;;) (i,j = 6,6)has the nonzero elements 


bii — — kiok20[Xis], b, = — ka [X27], 
b; =k, [Xo7], b33 = — kas|[Xos]o, 
bag = — ko[X23], 

bs4 =k22 [X23], bss = — koa, 

bes —kos[Xos]o; bos = — los, 


The active form of adenylyl cyclase, X27 :— AC", is obtained from 
[Xo7] (0) = Xz] — (X21) (0) (8.18) 


where [X1], = const is the initial concentration of the enzyme. We also assume 
that the level of phosphodiesterase enzyme, |Xos]; := PDE, remains constant 
throughout. 

The initial conditions provide concentrations of the reacting components 


X(0) — Xo (8.19) 


8.3 PLC-Dependent Pathway 


Activation of 43—receptors results in downstream stimulation of the intracellular 
phospholipase C (PLC)-protein kinase C (PKC) pathway. Four f), two y, four 6, 
and € isoforms of PLC enzymes have been isolated from the human detrusor. 
PLC/ members are triggered by Ca?* , but are differently regulated by G proteins. 
PLCf1 and PLCf4 are sensitive to Gag;;, whereas PLCf2 and PLC/3 can be 
activated by Gag}; and Gfy subunits. Without exception, their activation leads to 
the break down of inositide-4,5-biphosphate (PIP5) and the generation of second 
messenger molecules—inositol-1,4,5-triphosphate (IP3) and 1,2-diacylglycerol 
(DAG) (Frasier et al. 2008). IP3 is a highly soluble structure and it quickly diffuses 
through the cytosol toward the sarcoplasmic reticulum. Here it binds to Rip; 
surface receptors and triggers the mobilization of stored Ca?*. 
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Diacylglycerol consists of two fatty acid chains covalently bonded to a glycerol 
molecule. Compared to IP3, it is hydrophobic and therefore remains bound to the 
myometrial plasma membrane. DAG has a number of functions including the 
activation of protein kinase C, a subfamily of TRPC cation channels, and pros- 
taglandin production. DAG together with acyl-CoA is converted to triacylglycerol 
by the addition of a third fatty acid to its molecule. The reaction is catalyzed by 
two distinct isoforms of diglyceride acyltransferases. 

Protein kinase C comprises a family of eleven isoenzymes that are divided into 
three subfamilies—conventional, novel, and atypical—based on their second 
messenger requirements for activation. Conventional PKCs (o, Pı, f», y) require 
free Cx, DAG, and a phospholipid; novel PKCs (6, €, 0, u) need only DAG, and 
atypical PKCs (C, 1, 2) require none of them for activation. 

PKCs consist of a variable regulatory and a highly conserved catalytic domain, 
tethered together by a hinge region. The regulatory domain contains two subre- 
gions, namely C1 and C2. The C1 subregion has a binding site for DAG and 
phorbol esters and C2—acts as a Ca** sensor and is functional only in PKCs o, pı, 
f» and y. Upon their activation, kinases are translocated to the plasma membrane 
by RACK proteins where they remain active for a long period of time. The effect is 
attributed to the property of diacyglycerol per se. The enzymes play important 
roles in several signal transduction cascades. In detrusor, they phosphorylates 
MLCK and thus catalyze the contraction response. 

The simplified state diagram of the PLC pathway is outlined in Fig. 8.4. 

The corresponding system of equations is given by 


dX/dt = CX(t) + Co (8.20) 
where X(t) = (xj) ( — 29,34) has the components 


X», = [PIP;], X39 = [IP3], X31 = [DAG], 
X», = [DAGT], X33 = [PKC], X34 = [Rips]. 


The matrix C(cj;) (i,j = 6,6) has the nonzero elements 
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c11 = — kigkog[Xig], 


C21 =— Ci, cn = — ky, 

c31 —C21; c33 = — kag[X32] — kz9k30|X1] [X33], 
C44 = — ky [X31], 

55 = — kook30[X1] [X31], 

C66 = — kxj|Xao] — k-27, 


and the vector-column is: Co = (0, k 27 [X37], 0,0, 0, k-27 [X3;])T. 

Concentrations of the active form of protein kinase C,X35 :— PKC*, and the 
activated IP;—Ryp3—complex, X36 :— /P3 - Rips, can be obtained from the alge- 
braic relations 


(8.21) 


Here [Xa;]y. [X3g],are the initial concentrations of the IP5—receptor on the 
endoplasmic reticulum and protein kinase C enzyme, respectively. 

The dynamics of intracellular calcium release from the stores is described by 
(Miftahof et al. 2009) 


k3[X1] 
[Xi]? — ka’ 


d[Xi] 
dt 


— ko (à — h [X Bal (Cai) a (8.22) 


where ki (i= 0,4) are the kinetic parameters related to the release of sarcoplasmic 
Cat. 

SR 

Provided that initial concentrations of reactive substrates are known, the system 
of Eqs. (8.20-8.22) models the PLC pathway dynamics. 

The basic hypotheses in the model rely on data about axon cable characteristics, 
synapse geometry, electrical, and chemical mechanisms underlying signal transfer. 
The electrochemical response of the synapse to external stimulation can be pre- 
dicted accurately, assuming the accuracy of input data. In terms of the required 
number of parameters necessary to specify, the synapse model is quite complex. 
In calculations, we used the values of constants obtained from direct measure- 
ments made in in vitro experiments on different neurons. Some of the chemical 
kinetic parameters we have been unable to estimate from literature and they have 
been adjusted during the numerical simulation to fit experimental results well. 
Results of numerical simulations of the neurotransmission at a synapse are dis- 
cussed in the following paragraphs. 


128 8 The Synaptic Transmission 


8.4 Physiological Response 


The depolarization of the presynaptic membrane triggers a short-term influx of 
calcium ions into the terminal through Ca?* voltage-dependent channels. In 
Ata = 0.75 ms after the beginning of depolarization in the presynaptic zone the 
concentration of free cytosolic calcium quickly increases to 5.55 uM and reaches 
the max[Ca?*] = 19.4 uM in the next 0.65 ms. Some of the ions are immediately 
absorbed by the intracellular buffer system; others diffuse toward the vesicles 
where they bind to the active centers. Their activation triggers the release of the 
neurotransmitter, i.e., ACh, or NA,. This process depends on the dynamics of 
Car and is most rapid for 1.41 < tg < 2.13 ms. Thus at the beginning, the free 
fraction of ACh; rises at a rate of 0.5 mM/ms with a max d[ACh;]/dt. = 1.6 mM/ 
ms and of NA; at d(NA;]/drt. = 0.66 uM/ms and max d[NA;]/dr. = 2.87 uM/ms, 
respectively. During the whole cycle ~ 10 96 of stored vesicular ACh and ~ 2.1 96 
of NA are released. 

The neurotransmitter exit into the synaptic cleft starts as soon as 
[ACh;] 2 10.6 mM and [NA;] — 86.3 uM. Only half of free acetylcholine, 
5.38 mM, diffuses from the presynaptic terminal while the amount of NA, that 
reaches the cleft is comparatively higher to the norm: max[NA,] = 86.3 uM. 

The main part of ACh, reacts with the postsynaptic receptors with formation of 
(ACh-R)-complex, 0.11 mM: 99.2 % of receptors become active. In the case of 
the adrenergic synapse only 35.2 uM of NA, binds to 3-adrenoceptors, 
max[NA,-ARg] = 8.53 uM: 8.8 % of the whole number of available receptors 
become active. Another part of freely available NA, binds to presynaptic 
a-adrenoceptors. Their activation causes the autoinhibition of the neurotransmitter 
release. 

Highly reactive neurotransmitter-receptor complexes lead to the generation of 
the fast excitatory (EPSP), in the case of the cholinergic, and the fast inhibitory 
(IPSP), for the adrenergic synapse, postsynaptic potentials. The EPSP/IPSPs 
quickly attain their maxima: Vsyn = 87.1 mV and —76.2 mV, respectively 
(Fig. 8.5). The (ACh-R) and (NA,-ARg) complexes are very unstable and rapidly 
dissociate, releasing a receptive zone. As a result, the postsynaptic membrane 
returns to the resting state. 

The postsynaptic acetylcholine undergoes further fission by acetylcholine 
esterase enzyme (E): max[ACh,-E] = 0.47 mM is formed. The complex quickly 
dissociates into enzyme and choline which is reabsorbed by the nerve terminal and 
is drawn into a new cycle of acetylcholine synthesis. 

The effects of released NA are attenuated and terminated rapidly by reuptake of 
a large proportion into the nerve terminal (uptake-1 mechanism) where it is further 
taken up and stored in the dense-cored synaptic vesicles, while its free extraneu- 
ronal fraction in the cytoplasm runs the gamut of mitochondrial catechol-O- 
methyltransferase. 
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Fig. 8.5 The dynamics of the cytosolic Ca?*, vesicular and free acetylcholine fractions, 
concentrations and excitatory postsynaptic potential development in the case of physiological 
norm (1), after application of divalent cations (2), and botulinum toxin (3) 


8.5 Variations in Synaptic Neurotransmission 


The processes responsible for synaptic signal transduction can be affected by 
altering the intensity of stimulation, the depolarization of the presynaptic nerve 
terminal, the dynamics of voltage-dependent Ca?™ channels, the release of trans- 
mitters from the presynaptic terminal into the cleft, the interaction of neuro- 
transmitters with receptors, and the activity of enzymes. 
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8.5.1 Chloride Salts of Divalent Cations 


The chloride salts of Mg?*, Co?*, Cd^* , and Ni?* displace free cytosolic Ca2^ in 
the synaptic terminal without changes in the intensity of spontaneous transmitter 
release from the store. Thus, an increase in absorption of calcium by the cell buffer 
systems causes a rapid decrease in [Ca?"] = 5.4 uM (Fig. 8.5), and subsequently a 
reduction in the vesicular ACh or NA release. Most of the neurotransmitter remains 
intact; e.g., max[ACh,] = 79.5 mM, and only a small amount of its free fraction 
readily diffuses into the cleft. As a result, low amplitude EPSP/IPSPs are produced. 


8.5.2 Change in [Ca2'] 


A decrease in extracellular calcium reduces the rate of Ca?* ion flow into the 
presynaptic nerve terminal. Thus, a two-fold decrease in the concentration of C 
results in d[Ca?*]/dt = 10 mM/ms, and a ten-fold decrease — 1.5 mM/ms. Con- 
sequently, the amount of the free calcium available for binding with the active 
sites on the vesicles reduces to 9.2 and 2 (uM), respectively. 

The amount of discharged ACh, and NA, decreases. Thus, a two-fold decrease in 
Qu does not affect the dynamics of the release, diffusion, binding to and activation 
of receptors with the generation of the postsynaptic electrical signal. However, a 
tenfold decrease in Cg leads to a total block of synaptic signal transmission. 


8.5.3 p-Bungarotoxin 


p-BuTX and its derivatives suppress the diffusion of the free fraction of ACh into 
the synaptic cleft. The toxin at a “high” concentration causes a decrease in the 
diffusion velocity of ACh, from the terminal and its content in the cleft to 3.9 mM. 
This amount is insufficient to maintain the signal transduction in the synapse. In 
the case of a “low” dose of f-bungarotoxin the synapse generates 
Vsyn = 41.09 mV. Suffice it to note that the intensities of voltage-dependent Ca?” 
ion influx and ACh release from the vesicles remain unchanged and coincide with 
the data observed for the case of physiological norm. 


8.5.4 Botulinum Toxin 


Botulinum toxin blocks the active centers of Ca?* ions binding on the vesicles. 
Addition of the toxin suppresses the release of vesicular neurotransmitter. As a 
result blockage of the nerve-pulse transmission is observed. 
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8.5.5 Tetrodotoxin 


Artificial application of TTX in the vicinity of the presynaptic terminal has a 
detrimental effect on the inward sodium flux and, as a result, completely blocks the 
action potential generation and propagation in and through the affected zone. The 
nerve terminal remains in an unexcited state, and Ca?* ion flow is abolished. 
Intracellular content of free calcium remains at a low level, 0.1 mM, throughout. 
This amount is insufficient to sustain the process of the release of neurotransmitter 
form the stores. Subsequently, no synaptic signal transduction takes place. 


8.5.6 High-Frequency Stimulation 


High-frequency stimulation of the synapse by a train of 3, 10, and 20 impulses 
causes an increase in the amount of released neurotransmitter. Thus, in the case of 
action of 20 impulses, 2.2 mM of ACh, and 72.77 uM of NA, are released from 
the stores. The dynamics of exocytosis of ACh has its maximum of 0.5 mM/ms 
which is observed at the beginning of the process, followed by a sharp decrease of 
acetylcholine release to 0.001 mM/ms, while NA is released at a steady rate of 
1.67 uM/ms. The concentration of the free fractions of neurotransmitters quickly 
reaches its maximum, 10.6 mM and 99.1 uM, respectively, and they start to dif- 
fuse into the synaptic cleft. Results of simulations clearly show a quantum release 
of neurotransmitters. 

After stimulation by a series of twenty impulses max[ACh,] = 6.3 mM and 
max [NA,] — 66.2 uM reach the postsynaptic membrane where they bind to 
receptors. During the whole process max[ACh-R] = 1.224 mM and 
max[NA,-ARg] = 10.3 uM are formed. As a result, EPSP of an amplitude 
Vsyn = 89.4 mV and IPSP, V,,, = —81.2 mV are generated. 

Suffice it to note, that the dynamics of enzymatic and uptake-1,2 neurotrans- 
mitter utilization does not depend on the frequency of excitation. 


8.5.7 Inhibition of Neuronal Uptake-1 Mechanism 


Cocaine, imipramine, and sympathomimetic amines, e.g., metaraminol, dexam- 
phetamine, phenylethyl-amine, isoprenaline, and their derivatives, methylpheni- 
date, amantadine, chlorpromazine, and guanethidine, which are structurally 
analogous to the native NA, inhibit the neuronal uptake-1 mechanism in the 
adrenergic synapse. 

Application of cocaine causes an increase in the amount of diffused NA into the 
synaptic cleft, max[NA,] = 115.3 uM, without any affect on the time-character- 
istics of the process. The maximum concentration, 54.35 uM, of neurotransmitter 
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becomes available at the postsynaptic membrane. Approximately 13 % of the total 
number of f/;-adrenoceptors are activated with the max[NAp-ARg] = 12.2 uM. 

A dynamic equilibrium exists between neuronal and extraneuronal uptake 
mechanisms, i.e., when neuronal uptake-1 is inhibited more NA is utilized by 
extraneuronal uptake-2 and vice versa. Inside the cells, noradrenaline is exposed to 
the action of intracellular COMT enzyme. In the presence of cocaine, the amount 
of inactivated neurotransmitter increases, max[NA,-COMT] = 44.7 uM with a 
concomitant elevation of consumed enzyme, 92.3 uM. Finally, IPSP of the 
amplitude, V;,,, = —82.1 mV is produced. 
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Chapter 9 
Pharmacology of Detrusor Activity 


If (an experimentalist) had a needle to find in a haystack, he 
would proceed at once with the diligence of the bee to examine 
straw after straw until he found the object of his search.... A 
little theory and calculation would have saved him ninety per 
cent of his labor. 

N. Tesla. 


9.1 Classes of Drugs 


The majority of pharmacological agents used in clinical practice acts to alter the 
processes responsible for -transmission, by facilitating or inhibiting: (1) release, 
(2) enzymatic degradation of the neurotransmitter or modulator, (3) function of 
specific postsynaptic receptors, (4) second messenger system, or (5) intracellular 
regulatory pathways. For example, N-type calcium ion channel blockers—deriva- 
tives of w-conopeptides—interfere with the dynamics of cytosolicCa?*in the pre- 
synaptic nerve terminal. The decreased Qu concentration prevents activation of 
calmodulin protein and movement of vesicles containing neurotransmitters toward 
to presynaptic membrane. Chemical agents that facilitate cholinergic and adren- 
ergic neurotransmission can inhibit true and pseudo-acetylcholinesterase, mono- 
amine oxidase, and catechol-O-methyltransferase enzymes in the synaptic cleft. 

There are more than 20 families of receptors that are present in the plasma 
membrane, altogether representing over 1,000 proteins of the receptorome (Stra- 
chan et al. 2006). Transmembrane and intracellular receptors are being used as 
drug targets and they have a wide array of potential ligands. However, to date only 
a small percentage of the receptorome has been characterized. Abilities of a ligand 
to react with a receptor depend on its specificity, affinity, and efficacy. Selectivity 
is determined by the chemical structure and is related to physicochemical asso- 
ciation of the drug with a recognition (orthosteric) site on the receptor. The 
probability at which the ligand occupies the recognition site is referred to as 
affinity, and the degree at which the drug produces the physiological effect is 
defined as efficacy. Because different receptors are expressed on human detrusor, it 
is evident that disparate ligands acting alone or conjointly may elicit similar 
responses. 

A given receptor may contain one or more binding sites for various ligands and 
can be linked to different second messenger systems. The interaction of a drug 
with the site that is topographically distinct from the orthosteric site is called 
allosteric (Monod et al. 1965). The essential features of an allosteric drug—receptor 
interaction are: (1) the binding sites are not overlapping, (2) interactions are 
reciprocal in nature, and (3) the effect of an allosteric modulator could be either 
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positive or negative with respect to association and/or function of the orthosteric 
ligand. Thus, binding of a drug to the receptor changes its conformational state 
from the original tense to the relaxed form and hence, either facilitates or inhibits 
linking of the endogenous transmitter. 

Drugs that act at receptors are broadly divided into agonists and antagonists. 
Ligands that interact with the orthosteric site of a receptor and trigger the maxi- 
mum response are called full agonists. Related to them structurally are partial 
agonists, however, with lower biological efficacy. They are regarded as ligands 
with both agonistic and antagonistic effects, i.e., in the case of conjoint application 
of a full and partial agonist, the latter competes for receptor association and causes 
a net decrease in its activation (Kenakin 2004). In practice, partial agonists either 
induce or blunt a physiological effect, depending on whether inadequate or 
excessive amount of endogenous transmitter is present, respectively. 

Receptors that exhibit intrinsic basal activity and may initiate biological effects 
in the absence of a bound ligand are called constitutively active. Their function is 
blocked by application of inverse agonists—drugs that not only inhibit the asso- 
ciation of an agonist with the receptor but also interfere with its activity. This 
pharmacokinetic characteristic distinguishes them from true competitive antago- 
nist. Many drugs that have been previously classified as antagonists are being 
reclassified as inverse agonists. 

A class of drugs that have selectivity and affinity, but no efficacy for their 
cognate receptor are called antagonists. Antagonists that interact reversibly at the 
active site are known as competitive. Once bound, they block further association of 
an agonist with the receptor and thus prevent the development of a biological 
response. Ligands that react allosterically are called non-competitive. They stop 
conformational changes in the receptor necessary for its activation. A subtype of 
non-competitive ligands that require agonist-receptor binding prior to their asso- 
ciation with a separate allosteric site is called uncompetitive. Their characteristic 
property is related to the effective block of higher, rather than lower, agonist 
concentrations. 

Ligands that affect second messenger system function are classified according 
to the enzyme they act on. Thus, there are competitive selective and nonselective 
cyclic adenosine monophosphate, protein kinase A and C, phosphodiesterase, 
diacylglycerol, Ca**- adenosine triphosphatase, etc., activators and inhibitors. 
Because a single enzyme is often involved in multiple regulatory pathways, drugs 
of this category have a narrow, therapeutic index, and many side effects. Despite 
their pharmacokinetic and clinical limitations, they are widely used in laboratory 
research and represent a new approach in the management of inflammation- 
induced preterm delivery. 

All drugs, depending on the stability of a drug-acceptor complex that is being 
formed, show reversible or irreversible interaction. Reversible ligands have strong 
chemical affinity to a natural transmitter or modulator and normally form an 
unstable complex which quickly dissociates into a drug and a “receptor”. In 
contrast, irreversible drugs are often chemically unrelated to the endogenous 
transmitter and covalently bind to the target creating a stable complex. 
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9.2 Current Pharmacological Therapies of Bladder 
Dysfunction 


The control of detrusor excitability has important therapeutic implications since its 
malfunction may lead to bladder overactivity and urinary incontinence or urinary 
retention. These syndromes could be a result of pathological changes in the ana- 
tomical structure of the bladder and surrounding organs (e.g., benign prostatic 
hyperplasia, bladder calculus, developmental abnormality), the central and 
peripheral nervous system that are involved in control of the bladder (e.g., brain 
and spinal cord lesions, sacral nerve damage), the detrusor (e.g., hyperreflexia, 
instability or areflexia, intramuscular fibrosis, ischemia), coordinated function of 
the external sphincter and the detrusor muscle, urethral sphincter bradykinesia, or 
be a complication of multiple sclerosis, Parkinson disease, dementia, spinal cord 
injury, and myelomeningocele among others. The conditions cause a health, 
hygienic, and social concern to the patient. Thus, patients with neurogenic bladder 
are highly susceptible to urinary tract infections, and those with spinal cord 
injury—to stone formation in the bladder and kidneys and bladder cancer. 
Unfortunately, no adequate pharmacological modalities are available at the 
moment to treat them. 

For example, overactive bladder symptoms are experiences by approximately 
17 % of both men and women over the age 40 years and are becoming even more 
common—33 %—with increasing age (Stewart et al. 2003). Apart from the 
incalculable costs associated with the physical and emotional suffering inflicted by 
overactive bladder, especially in the wet patient, the disorder exacts a heavy 
economic toll on individuals, their families, and society. An estimated total eco- 
nomic cost of the disorder, which includes both direct—the diagnosis and treat- 
ment—and the indirect costs—lost wages and productivity, is $12.6 billion. The 
mainstay therapy of overactive bladder rests on the use of ji and u3 receptor 
antagonists—oxybutynin, darifenacin, fesoterodine, propiverine, solifenacin, tol- 
terodine, and trospium. Although the drugs offer a significant improvement in the 
management of symptoms their overall efficacy, tolerability, and toxicity are still 
less than optimal. Potential novel formulations such as selective %,p-adrenoceptor, 
serotonin 5-HT;4, and 5-HT;, neurokinin NK; and NK» receptor antagonists 
(SP-saponin, aprepirant), /3-adrenoceptors (terbutaline, clenbuterol), and potas- 
sium channel agonists (retigabine, forskolin), vanilloids (capsaicin, resiniferatox- 
in) botulinum toxin are currently under investigation (Sval@ et al. 2012; 
D' Agostino et al. 2006; Hristov et al. 2008; Apostolidis et al. 2009; Sacco et al. 
2008; Staskin et al. 2011). Despite the large amount of new biological insights, 
very few drugs with mechanism of action other than antimuscarinics have passed 
as yet the proof-of-concept stage. 

Successful drug discovery requires deep understanding of the mechanisms of 
diseases, the full biological context of the drug target, biochemical mechanisms of 
drug action. It should involve a multilevel conceptual framework which would 
allow the integration and variation of parameters and constants within the 
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biological system with high precision. Systems (computational) biology approach, 
as a thorough, quantitative and qualitative, interrogation of biological processes 
within the physiological milieu in which they function, provides a new paradigm 
to study the combined behavior of interacting components through the integration 
of experimental, mathematical, and computational methods. 

The development of drugs for clinical use calls for standard test models in 
which a large number of substances can be analyzed simultaneously to determine 
whether, and how, drugs contribute to the overall physiological response. Math- 
ematical modeling and simulations in the form of virtual laboratories offer a 
uniform platform for bridging gaps in our understanding of intricate biological 
mechanisms and generation of new hypotheses that are not realizable within a real 
laboratory. Verified experimentally, they serve as predictors in the iterative pro- 
cess of identification, validation, and evaluation of drug targets, and supply a 
valuable insight into questions of drug development, efficacy, and safety. 


9.3 Model of Competitive Antagonist Action 


Let La,be the competitive reversible antagonist. A part of the state diagram of 
signal transmission in the myometrium that describes ligand action is shown in 
Fig. 9.1, while other reactions in the general cycle remain unchanged. 

Assuming that the reactions of association/dissociation of the antagonists (Lat) 
with the receptor (R) and the drug-receptor complex (La; - R) formation satisfy the 
Michaelis—Menten kinetics, the governing system of equations is 


dXA /dt == DaXat(t) + Coat: (9.1) 


The matrix Da; is the extension of the matrix D (Eq. 8.10) 


Fig. 9.1 The state diagram 
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where the modified elements and new elements are given by 


5,21 = —ka|Xs], 4522 = k-31 + k32, 
435,35 = —k3i[Xs], 35,36 = k-31, 
d36,35 = k31 [Xs], | d3636 = —(k-31 + k32). 


New components of the vector Xa are define: Xy :— Lar, X43 :— Lar: RX43 : 
= La R. To close the system, it should be complemented by initial values for La: 
and La, - R. 

In case of partial chemical equilibrium, which could be achieved after pro- 
longed treatment of the detrusor with the antagonist, the dynamics of the receptor, 
drug, and drug-receptor complex conversions can be described by, 


Dese] = — ka |Xs] [X42] + k-31 [X43] (9:2) 


ae =k) [X5] [X42] — k-31 [X44]. 


Here [X45], [X43] are equilibrium concentrations of the drug and the bounded 
complex, respectively. Summation of equations for [X5], [X43] yields 


d|Xs] | d[X43] 


= 0 or [Xs](t) + [Xa] (t) = const, (9.3) 
from where letting the concentration of total available receptors [X5])= constant, 


we get 


[Xs](r) = [Xs]o— [Xas] (0). (9.4) 


Since d(L,, - R)/dt = 0, substituting (9.4) into the second equation of (9.2) and 
after simple algebra, we obtain 


[X43] (t) mmm wp? 


(9.5) 
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Here K* = kzı /k-31is the Michaelis-Menten equilibrium constant. 
Finally, the dynamics of receptors in presence of a competitive antagonist and a 
corresponding endogenous transmitter (L,) is given by 


[Xs] (2) = [Xs]o—[Xas] (0) — [Xo] (0) — [Xs] (4). (9.6) 


Here [Xs](t) is the concentration of constitutively active receptors, and [Xs] (t) 
is the equilibrium concentration of the L, - R— complex. 

Substitution of Eq. (9.6) into (9.1) allows some simplifications in the governing 
system of equations. We have left it as an exercise to the diligent reader. 


9.4 Model of Allosteric Interaction 


Two broad conceptual views underlie the majority of studies of allosterism. The 
first, was developed initially in the field of enzymology, is based on the assumption 
that proteins possess more than one binding site that can react successively with 
more than one ligand. The second considers allosterism as the ability of receptors 
to undergo changes that eventually yield an alteration in affinity of the orthosteric 
sites for endogenous transmitters (Monod et al. 1965; Koshland et al. 1966). 

Respectively, two types of mathematical models, i.e., concerted and sequential, 
along with their various expansions and modifications, have been proposed to 
simulate cooperative binding. The concerted model assumes that: 


]. Enzyme (receptor) subunits in equilibrium attain the identical—tensed or 
relaxed—conformation, 

2. It is affected by an allosteric effector, 

3. A conformational change in one subunit is conferred equally to all other 
subunits. 


In contrast, the sequential model does not requires the satisfaction of conditions 
(1), (3) but instead dictates an induced fit binding of the ligand with subsequent 
molding of the target, instead. 

Consider a modified part of the general state diagram of allosteric ligand- 
receptor interaction (Fig. 9.2). Positive non-competitive allosteric mechanism 
assumes binding of the ligand Lanto the receptor in the R' conformation. 

The Lay - RR— complex further associates with the endogenous transmitter 
Land produces the active complex— Le - RR. Lan. In case on uncompetitive 
positive allosteric mechanism, though, the transmitter Le binds first to the RT— 
receptor, changes its configuration to Le - R form and only then the ligand Lay 
occupies the allosteric site. The Le - RÈ - Lan— complex reacts with the G protein 
system and enters the cascade of chemical transformations as described above. 

Comparison of the unperturbed and current schemes gives: 
X; := R™,X6 := Lay RÈ, X7 := Lp R®- Lan , Xg := Lo: RT, Assuming that all 
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reactions satisfy the Michaelis-Menten kinetics, the system of equations for 
allosteric interaction is given by 


dXan/dt = DauXAn(t) + Coan. (9.7) 


Here 

Xan(t) = (X1,...,X34,Xa4)" , Coat = (0,...C30,0,...0,C34,0)" , Xas := Lan. 
The matrix Day contains the modified matrix D (Eq. 8.10) and the additional 

new elements 


d635 = k31[Xs], 

4735 = k32[Xs], 

ag35 = —(k_6 + kx2) [Xs], 

d35 4 = ks[Xo] + ks|Xs], das; = k-7, dass = —k3ı [X5] (k-6 + k32) [Xs]. 


In contrast to competitive antagonists which cause a theoretically limitless 
rightward shift of the dose-occupancy and dose-effect curves for endogenous 
transmitter, allosteric ligands attain a limit which is defined by the binding factor. 
Thus, allosteric agonists, applied conjointly with agonists and endogenous trans- 
mitters, can enhance their spatial and temporal selectivity at given receptors and 
guarantee a required level of safety. 

Recently, a precursor of a drug, THRX-160209, an allosteric modulator of u2 
receptor has been investigated for its potential properties to treat overactive 
bladder. Experiments using radioligand binding assays to assess cooperativity and 
affinity suggested that the compound binds in a multivalent manner to the u2 
receptor, simultaneously occupying the orthosteric site and a spatially distinct 
allosteric site. 


140 9 Pharmacology of Detrusor Activity 


9.5 Allosteric Modulation of Competitive Agonist/ 
Antagonist Action 


One of the intriguing pharmacological properties of allosteric drugs is their 
potential ability to alter selectivity, affinity, and efficacy of bound and non-bound 
competitive agonists/antagonists by enhancing or inhibiting their cooperativity at 
receptor sites. During the last decade, effects of different allosteric compounds 
have been studied extensively in vivo and in vitro. For example, it was found that 
gallamine diminishes the affinity of bound acetylcholine and inhibits its negative 
ionotropic and chronotropic effects in the myocardium, while alcuronium exerts 
positive allosteric modulation on the affinity of ACh (Stockton et al. 1983; Tucek 
et al. 1990). Studies with the human adenosine A, receptor revealed the diverse 
regulatory effects of the allosteric modulator (PD81, 723) on the affinity of a 
partial agonist (LUF5831), a full agonist (N°-cyclopentyladenosine (CPA)), and an 
inverse agonist/antagonist (8-cyclopentyl-1,3-dipropylxanthine-DPCPX)), for the 
receptor. Results demonstrated that it increased the affinity of CPA, slightly 
decreased the affinity of LUF5831 and significantly reduced the affinity of DPCPX 
(Heitman et al. 2006). Thus, therapeutically, allosteric ligands are capable of 
modifying signals carried by the exogenous and/or endogenous ligands in the 
system. 

Let the detrusor muscle be exposed simultaneously to a competitive agonist 
(antagonist), Lag(ant); allosteric ligand, Lay, and the endogenous transmitter, Le. 
The proposed state diagram of their interactions is shown in Fig. 9.3 

It combines noncompetitive and uncompetitive allosteric mechanisms of action 
which involve binding of Lay to the receptor RT, formation of the Lay - RF— 
complex followed by binding of the agonist (antagonist) L4,(,,; and the trans- 
mitter L. to it, and an inverse sequence, i.e., binding of L4,(,) and L. to the 
receptor RT first with the subsequent addition of the ligand Lay. As a result of both 
processes, the Le - LAg(ant) - RF - Lay active (inactive) complex is produced. 

The governing system of equations for positive noncompetitive allosteric 
modulation of the agonist (antagonist) action is 


dX am/dt = DamXam(t) + Co. Am (9.8) 
Fig. 9.3 The state diagram Ko( Lag) 
of the effects of a competitive k 
agonist/antagonist and AE. 
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The vector of reacting components is Xam(t) = (ace CP X42; -< TD 


where X; := RT, Xo :— Lag: RÈ, X7:= Lo- RP- Lan: Lag Q Xs := 
Lo- R" -Lag , Xy := Lag, Xa := Lag: R", Xyu := Lan , Xas := Lanc RÈ- 
Lag , X46 := Lan - RF. Le. The meaning of other components is as described in 
Eq. (8.24) The vector of constant concentrations of substrates is 
Co,Am = (0,...C30, 0, ...0, C34,0, . . .0)7. 

The matrix DA, is obtained from the general matrix D where the following 
elements are adjusted and new elements are introduced 


d44 =k2( [E0] — |Xs]) + ke[Xs] + (k- + ks + ks) [Xe] 
Hk ai[Xio] + (k-8 + kg) [X45], 

d442 =k34|X6] + k-6[Xa6], a444 = kai [Xs], 
ds,42 —ks|Xe], as44 = —kai[Xs], 

a4 —kg[Xas] — ke[Xe], ass = kai [X44], 

ass = — k_3[X4] — (ks + k-iokio[Xis]) [X12] — ks ([X42] + [X4]), 
6,42 —k.e|Xas] — ksa[Xo], 
4745 =k_g[X4] + ka[X42]; 
d4»4 =(ke + ks)|Xas], | daos = k31[Xaa], 
di»; —k.6,  da242 = (ks + k34)[Xo] + (k-6 + k34) [X46], 
djs —ks|Xao],  daa.aa = ks [Xo] — kai [Xs], 
dis =k34[X42], das; = kg, — dasas = — (k-s + kg) [Xa], 
d46,6 —ke|Xa], dae; = ka, — daea4e = —(k-6 + ka4) [Xap]. 


In case of uncompetitive positive allosteric modulation the vector of reacting 
components is Xam(t) = (X1,...,X34, X42, .. X47)", where X47 := Le- RT and 
the meaning of other components is as described above. The vector of constant 
concentrations of substrates remains unchanged. The matrix Dam has new 
elements 


44 —k»([Eo] [Xs]) + ks [Xs] + (k-8 + ks + kc) [Xo] + 


d442 =k_6[Xa7], a444 = k-6 [Xe], 

54 =k_6[X43], ass = (Kk. + ko) [Xa] — Ka [X42], 
5,42 =k_6[Xa7], 47,44 = k32[X], 

ags = — k-e — k-7 — (k-6 + k32)[Xa4], agaz = Kas[Xa], | 8843 = ke[Xa4], 
dya = — k33[Xs5] + (k-6 + k33) [X47], | daos = k-e[Xa], da244 = k-6 [Xs], 
daz a2 =k3|Xs], | daas = —(k-6 + ke)|Xa], | das44 = ko [Xs], 
day —k.3,  daaar = k33[Xa7],  dagag = —(k-6 + kz) [Xe] + ko [Xa]. 


Initial concentrations of the reacting components close the system. 
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9.6 Effect of Cholinergic Antagonists/Agonists 


The addition of ACh-antagonists causes a deactivation of chemical processes on 
the postsynaptic membrane. These effects are dose dependent. According to the 
numerical results, the compound in concentrations of 0.13-0.26 mM causes a 
complete blockade in excitation transmission, V,,, = 0; while at lower concen- 
trations 0.026 mM, a decrease in the number of choline receptors which are 
available for interaction is [R] = 17 uM. Subsequently, the concentration of 
(ACh-R) complex is reduced to 0.016 mM, which, however, is enough to produce 
a depolarization of the postsynaptic membrane with the generation of EPSP, 
Vsyn = 69.6 mV. 

In the case of a partial chemical equilibrium the compounds in a concentration 
of 0.26 mM, which is two times higher than the whole number of ACh receptors 
on the postsynaptic membrane, does not cause the blockade of the synaptic 
transmission. The maximum of [ACh-R] = 1.87 uM is formed, which leads to the 
generation of EPSP of an amplitude of 53.2 mV (Fig. 9.4). 

In addition to the drugs that completely or partially block synaptic transmission, 
there is a large group of chemical agents that facilitate cholinergic transmission. 
Their primary mechanism of action is to inhibit the activity of acetylcholinesterase 
enzyme in the synaptic cleft, and to increase the effectiveness of ACh release from 
the nerve terminals. 

Strictly reversible anticholinesterase agents, such as galantamine, ambenonium, 
endorphonium and aceclidine, are not hydrolyzed by cholinesterase and relatively 
quickly dissociate intact from the enzyme. Others, such as neostigmine and its 
derivatives, remain bound to the esoteric site of enzyme for a long time and 
undergo the reaction of hydrolysis, just as ACh does, but at a very slow rate. 
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Fig. 9.4 Receptor availability on the postsynaptic membrane and excitatory postsynaptic 
potential development in the case of physiological norm (7), in the presence of a cholinergic 
antagonist (2), and in the case of chemical equilibrium (3) 
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The treatment of the nerve terminal by cholinesterase inhibitors in concentra- 
tion 0.1 mM causes a significant decrease in the active enzyme content in the 
synaptic cleft [E] = 3.5 uM. The presence of the inhibitor leads to an expected 
accumulation of free fraction of ACh in perfusate: max[ACh,] = 5.42mM. The 
amount of inactivated neurotransmitter in the cleft is [AChE] = 0.22 mM 
(Fig. 9.5). Consequently, 0.124 mM of (ACh-R) complex is developed on the 
postsynaptic membrane which induces the generation of EPSP of an amplitude 
89.08 mV. In the case of the addition of 0.05 mM, Vsyn = 88mV is registered. 

Interestingly, in the case of a partial chemical equilibrium between the inhibitor 
and the enzyme, the application of 0.15 mM of the drug, which is three times 
higher than the concentration of acetylcholinesterase, does not abolishes the 
enzyme activity. The concentration of free enzyme falls to 1.3 uM. A slight 
increase in ACh content in the cleft, max[ACh,] = 5.44mM, and a significant 
decrease in inactivated acetylcholine, [AChE] = 1.22 uM, are observed. The 
amplitude of EPSP attains 89.1 mV. 


9.7 Inhibition of Catechol-O-methyltransferase 


Catechol, pyrogallol, quercin, rutin, substances described as flavonoids of the 
vitamin P group, tropolones, and a group of the benzoic acid derivatives are potent 
COMT inhibitors. Their addition causes a significant increase in the concentration 
of NA; on the postsynaptic membrane. The response is nonlinear and depends on 
the concentration of compound added. Thus, in the case of treatment of the 
adrenergic neuron with 0.15 mM of a drug, [NA,] = 58.2uM is formed, and after 
application of 0.2 mM of inhibitor, max[NA,] = 80.5 uM is produced (Fig. 9.6). 
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Fig. 9.5 Changes in acetylcholinesterase enzyme and inactivated acetylcholine in the synaptic 
cleft as a result of application of cholinesterase inhibitors in concentrations 0.1 mM (2) and 
0.15 mM (3). The curve (1) corresponds to the case of physiological norm 
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Fig. 9.6 Changes in noradrenaline and noradrenaline-receptor complex on the postsynaptic 
membrane in the presence of catechol-O-methyltransferase inhibitors. The curves (1), (2) refer to 
the concentrations of 0.13 and 0.2 mM of the inhibitor present in the state of chemical 
equilibrium, and the curves (3), (4) correspond to addition of 0.15 and 0.2 mM of the inhibitor, 
respectively 


As a result [NAp — ARg] = 15.5 and 19.8 (uM) are formed and 18.8 versus 
13.8 uM of the f/;-adrenoceptors are bound, respectively. The amount of utilized 
NA by COMT also decreases to 25.2 uM. After treatment with a “high” dose of 
inhibitor, 0.2 mM, only 20.04 uM are produced. These affect the level of hyper- 
polarization of the postsynaptic membrane: V,,, = —83.3mV (—78.7 mV) are 
recorded. 

In the case of a partial chemical equilibrium between the inhibitor and the 
enzyme the addition of IK in concentrations of 0.13 and 0.2 mM causes an 
increase in [NA,] to 119.3 and 135.5 (uM), respectively. Consequently, the con- 
centration of bound NA, to adrenoceptors on the postsynaptic membrane increa- 
ses: max[NA, — ARg] = 24.5 and 26.98 (uM). The blockage of 81-87 96 of 
COMT leads to a decrease in the intensity of degradation of NA,. Only 9.59 and 
13.0 uM of the (NA, — COMT) complex is formed. As a result, the amplitude of 
IPSP increases: Vsy, = —85.6 and — 86.1(mV). 


9.8 Effect of f: Adrenoceptor Antagonists 


Drugs that block the adrenoceptors for sympathomimetic amines are classified on 
the basis of their chemical structure: haloalkylamines (phenoxybenzamine), imi- 
dazolines (phentolamine, tolazoline), phenylalkylamines (thymoxamine). The 
mechanism of their action is rather complex and not completely understood. It is 
supposed that the antagonism is competitive and slowly reversible due to the high 
strength of the covalent bond formed with the receptor. As well as acting on 
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postsynaptic receptors, they can affect adrenergic transmission preventing the 
uptake of released NA into neuronal and extraneuronal sites, and inhibiting 
adrenoceptors on the presynaptic membrane. The second mechanism is more 
important since the compounds that inhibit NA uptake cause only a slight increase 
in NA release. 

The addition of antagonist in rising concentrations blocks 55.4—67.5 % of the 
receptors on the postsynaptic membrane, respectively. The concentration of free 
f3-adrenoceptors also decreases to 44.6-33.51 (uM) and, as a result, the rate and 
the maximum concentration of (NA, — ARg) complex development falls. Thus, 
max[NA, — ARg] = 4.81uM is formed (Fig. 9.7). The inhibitor at a “high” dose 
causes further decrease in the amount of active complex, 3.49 uM. The above 
changes affect the postsynaptic response: IPSPs of the amplitude V;,, = 
— 11.42 and — 66.24mV are produced. 

After the partial chemical equilibrium is achieved, an added drug in concen- 
trations 0.2 and 0.3 (mM) blocks 69.2—75 % of f;-receptors. This is associated 
with a decrease in formation of [NA, — ARg] 3.12 and 2.36 uM. As a result, a 
higher level of hyperpolarization is attained: Vsyn = —64.3 and — 58.86(mV). 
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Chapter 10 
Human Urinary Bladder as a Soft 
Biological Shell 


Start doing what’s necessary, then what’s possible, and suddenly 
you are doing the impossible. 
St. Francis of Assisi 


10.1 Basic Assumptions 


In Chap. 3 we have demonstrated that the urinary bladder satisfies all the 
hypotheses and assumptions of the theory of thin soft shells. Furthermore, the wall 
of the organ, its geometry, and stress-strain distribution within it can be approx- 
imated, with a sufficient degree of accuracy, by the middle surface of the organ. 
Contractions in the bladder are triggered by underlying electrical events in the 
detrusor. The basic anatomical and morphological data can be summarized as 
follows: 


1. the wall of the bladder is composed of two interspersed muscle layers— 
longitudinal and circumferential—with fasciculi running orthogonally within 
them; 

2. the detrusor has the properties of the electrically isotropic syncytium; multiple 
gap junctions provide the low electrical resistance gating among smooth muscle 
cells; 

3. the electrical activity represents the integrated function of ion channels: volt- 
age-dependent Ca?* channels of L-and T-types, large Ca^* -activated K* 
channels (BKc,), potential sensitive K* channels, and CI" channels; the 
properties of the channels are neurohormonally modulated; this effect is 
assumed to be mainly chronotropic with an increase in the time of permeability 
for specific ions; 

4. although pacemaker cells per se have not been found in the organ, there are 
groups of cells that are responsible for generation of electrical signals and 
coordinated rhythmic contractions of labor; 

5. contractions are dependent of external and internal control; the coordination is 
provided by multiple neurotransmitters and hormones; 

6. the key player in the process of electromechanical coupling is free cytosolic 
calcium; the active force 7" - Ca?" -activity relationship is characteristic for the 
smooth muscle and satisfies Eq. 6.6; 
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7. the total force T;.; generated by the muscle layers is the result of deformation of 
its passive elements, T” (2e, 4), cj), i.e., collagen and elastin fibers, and active 


contraction-relaxation of myofibrils, T° (2e, z&, [Ca2*], ci) 


Toy = TP (2e, 21, ci) + T° (444, ZG), [Ca?*], ci), (10.1) 
where the meaning of parameters and constants as described above; 


8. electrochemical and chemoelectrical coupling in the system is guaranteed by 
separate and/or conjoint activation of multiple intracellular pathways, and 
ligand-dependent transmembrane ion channels. 


The above data, 1-8, have been incorporated in a mathematical model of the 
human urinary bladder as discussed below. 


10.2 Model of the Bladder 


0 
Let S be the cut configuration of the bioshell that describes the anatomical shape of 
the organ with accuracy of bending in absence of loads. Assume that when the 
organ is empty of urine its cut configuration coincides with the undeformed, S, 


0 
configuration (S = S) and: 


1. stretch ratios 4;; = 1.0,(i, j = 1, 2); 
2. in-plane total forces T;;(A;;) = 0 throughout the bioshell; 
3. intravesicular pressure p = 0. 


The bladder attains the deformed S state with 4;; > 1.0, Ty > 0, V 7 0 and 
p > 0 during the filling phase. 
Let the middle surface of the bladder be associated with a cylindrical coordinate 


system {r, 9, z} 
r = r(0,05,t), 
9 = (%1, %,¢), (10.2) 
d z(a, 05, f); 


that is related to the orthogonal Cartesian coordinates (xi, x2, x3] as 
Xj = rcosQ, X2 = rsin, x3 = Z, 
The position vector 7 of point M(r, g,z) € S is given by 


T = rk, + zk,, (10.3) 
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where 

k, = i cosọ +i,sing, k, = —i, sing +i, cos q. (10.4) 
Differentiating Eq. (10.3) with respect to o; with the help of (10.4), we obtain 


1 


or Or, _ ok Oz z 


"em Og Oa, | Ou Oa ° (10.5) 
wea $ D "S. | 
r 
aj Oxi Qo; ^' 
where 
Ok, 0 (i i gi ) 
= =x i, cos i, sin 
Qa; Qa; 1 p 2 p (10 6) 
0 Zoe s 0g - l 
= =< (-i sing + i, cos p) = T 
Projections of 7; in the direction of the r, , and z axes are given by 
H Or 0p |. Oz 
= — yf. = = i => —, 
" — Qa "^ — Qa; * — Qaj 
Hence, 
ri = Tikk, + Tipk, + Tick, . (10.7) 
Decomposing e; along the base (kk b, kets we find 
e; = liki + ligko + licks, (10.8) 
here the direction cosines l;(i = 1,2,3; j = r,@,z) are given by 


2122 = lyb, + hole + lil = cosy = at,/s/ai ah 
= ly, + lile + lila: = 0, (10.9) 


£s = E +Ê, Ael. (n= 1,2,3), 


a 
1 
w 

l 


Expanding the resultant of internal uterine pressure, p, and the mass force, f, in 
the direction of unit vectors @,, we get 


p = epi + &2P2 + @3p3, (10.10) 
f = kf. + Kf + kfz. (10.11) 


The vector of acceleration a(a,, do, az) in cylindrical coordinates is given by 


a, = — 


2 2 p 10.12 
de "ar + 2r a (10.12) 


dr (de do dr dọ d?z 
dg = —r 
m dt? dt dt’ ^ d? 
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Substituting Eqs. (10.8), (10.10—10.12) in (4.63), the equations of motion of 
the soft shell in cylindrical coordinates take the form 


dr (dp\*\ o - 
= Tlie + Ty 
oval $3 (2) Oo í 1H € 2r) an2 
0 T T * vi 
dp ex; (Th, F Tl) aiil + (pilir + pol) Va 
"D (bolz = PANTE + pf. a, 


BU eue cs ; 
eve apt? B 0n [rns ERSTE 


ð if * * 
4 Bap [o s T? loo) ail + (pilio anv 


+ pa(hzb, — lirl2z) dii ax + pf, v'a, 


Taz z : 9 Hm A 1 
pva =x | (Tih Tits) an| + az (rts Th) | 


Qo 


+ (Pilz pala) a+ P3 (liso = holo) V an 45) + pfa. 


(10.13) 
In case of pacemaker activity, V,, its dynamics is given by 
dV, 
ACh T Sol laxi); (10.14) 
Í 


where Ch is the specific smooth muscle cell membrane capacitance, J; are 


respective ion currents (j = L- and T-type Ca**, BKc,, and voltage-dependent 
Na‘, K*, and CI), Lexy i) = Vi /Rp is the external membrane current, and R, is the 
input cellular resistance. The currents are defined by 


ris &ca(i)zca (Vp — Vea) 8Ca-K(i)Poo(Vp — Vea-x) 
"1+ 8ca[Ca?*] C 0.5 + [Ca] 
Iya = 8Na(i) Mga (Vp = Vna), Ik = grong (Vp = Vk), 
Iq = Zera (Vp = Va); 
where Ve, Vc, k; Vna: Vio Ver and 86,5. Sc. (y Sway Eka Sa ate the 


reversal potentials and the maximal conductances of the respective ion currents, Jca 
is the numerical parameter, and the turnover of free intracellular calcium yields 


Ica-K = 
(10.15) 


d[Ca?*] | 0.2zca(Vp — Vea) 
dt 1+ Dca[Ca?*] 


0.3[Ca?*] (10.16) 
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The dynamic variables of the ion channels zca, p,,. Mya, Aya, and ny are given 
by 
dzc,/ dt = (Zoo = ZCa)/Tz, 
dhya/dt = An (Pis — hna)/Th, 


die) di 3L to fi) s VO 
mya = ma(Vj), Poo = [1 + exp0.15(V, + 47)] |. 
In the above m,,, 1, ,n,,,z,, are calculated as 
yoo 
Yoo ~= 73 (y = m,h,n), 
yoo d Byes (10.18) 
Zoo = [1 + exp(—0.15(V + 42))]7", 
where 
0.12(V, + 27 
moo = ( P ) Pma = 4exp(—(Vp aE 47) /15), 


1 — exp(- (V + 27) /8) ' 
hoo = 0.07 exp(—(Vp + 47)/17), Broo = [1 + exp(- (V, + 22)/8)] , 
|. 0.012(Vp + 12) 

1 — exp(- (V, + 12) /8)’ 


noo 


Broo = 0.125 exp(—(Vp + 20) /67). 


The propagation of the electrical wave of depolarization, V; ,, along the detrusor 


gl 
is given by 
ov: My 0 90. oV? 0 80.4 oV? 
j= obi) fe (tes, qiiis 
Ot 1+ py, (Oa, X Ae Ooty Ox \ À O0 


Here £o) have the meaning of maximal intracellular and interstitial space 


conductivities, %1, 0» are Lagrange coordinates, and the meaning of other param- 
eters as described above. 
The total ion current Jon is given by 


lion = nañ (V7, = Vna) T gk (V3, DS Vk) ap &a(Vz, = Vo), (10.20) 


where 2y,,x.£c; represent maximal conductances, Vna; Vx, Vc;-reversal poten- 
tials of Na*, K^, and Cl” currents. The dynamics of the probability variables, the 
activation and deactivation parameters of the ion gates are given by Eqs. (6.14) 
and (6.15), respectively. 

Neurohormonal regulatory processes in the detrusor are described by 


dX/dr = DX(r) + Co, (10.21) 
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where the matrix D and vectors of reacting substrates X and C, are given by Eqs. 
(8.9), (8.17), and (8.20) and the algebraic relationships (8.18), (8.21). 

The following anatomically and physiologically justifiable initial and boundary 
conditions assume that: 


9. the initial undeformed configuration of the urinary bladder is defined by the 
intravesicular pressure, p, and the organ is at rest, i.e., electrically quiescent; 


10. the excitation of known intensity, V, and duration, d , may be triggered by 
electrical discharges in the a priori defined pacemaker regions; 

11. concentrations of the reacting substrates are known; 

12. the bladder sphincter is rigidly fixed throughout filling. 


Thus, 
0, O«r«n 

Att —0: V,—4p > Ver = Vin Ve, = 0, 
Vp du 


0 
[Ca?*] = [ca |, X(0) = Xo, (0) = Co, 
r(&1, 2) = ro(&1, 02), P(%1,%2) = Po(%1, %2), z(&1, %2) = Zola, 92); 


and the dynamic variables of ion channels involved are 


m — Mo, h= hy, n = næ, ZCa = ZCaoo; hya = hwasc; nk = nKo, 


h— hy, it = ño, Xca = Xcasc, ed 
fort >0 
r(%1, 02), 0 = ro(0, %2) = ro(L, %2), 
olai, 02), Lo, = Po(0, %2) = Po(L, #2), (10.23) 


z(01, %2)y, 0,2 = 20(0, %2) = zo(L, %2), 
Ve (0, %2) = Ve (L, %2). 


To close the problem, the above equations should be complemented by con- 
stitutive relations. Finally, the mathematical model describes: 


1. the urinary bladder as a soft electromechanically active biological shell; 
2. the dynamics of electromechanical processes and neurohormonal regulation of 
contractions in the organ. 
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10.3 Bladder Mechanics 
10.3.1 Filling Stage 


The results of numerical experiments of the study of the statics and dynamics of 
changes in shapes, force-stretch ratio distributions, electromechanical activity, and 
pharmacological modulation of the urinary bladder are presented in the following 
paragraphs. It must be recognized that the model contains numerous parameters 
and constants that have not yet been evaluated experimentally. For example, no 
information is available currently on mechanical constants c; of the human bladder 
under biaxial loading, chemical reactions k;,, and electrical properties of the 
detrusor as an electroconductive medium. To overcome the problem of missing 
data, during simulations, the wall of the organ is regarded as a curvilinear 
anisotropic nonlinear viscoelastic material with constitutive parameters derived to 
resemble similar soft human tissues. Thus, the uniaxial passive force-stretch ratio 
T? (A, cj) relationship assumes to show “pseudoelastic” behavior; the active force 
function exhibits nonlinear dependence on the intracellular calcium concentration; 
and the biaxial in-plane T(4;, 4;) function is a generalized form of the uniaxial 
constitutive relationship. Other parameters and constants are adjusted to resemble 
the physiological and diseased states of the organ. Therefore, the results of 
modeling are not aimed to achieve an accurate quantitative representation, but 
instead offer a qualitative assessment of complex processes in the bladder during 
filling and voiding. 

Anatomical shapes of the human urinary bladder are highly variable. An 
example of a three-dimensional CT image reconstruction of the filled urinary 
bladder is shown in Fig. 10.1. It can be imagined that the fundus is located in the 
lower pelvis and is significantly restricted to movements, while the body and the 
apex undergo extensive deformations and changes in configuration. 

The cut shape of the bioshell corresponds to the electrically quiescent organ 
empty of urine. It can attain many forms simultaneously and depend on the status 
of the tissue per se and surrounding organs. The initial undeformed configuration 
of the bladder though is determined by the amount of urine it contains at the given 
moment of time and the intravesicular pressure, p. In simulations we have assumed 
that at f = 0 the bioshell is loaded by p = 1.0 kPa. Although urine, as a second 
phase, is not present in mathematical formulation of the problem, its mechanical 
effect has been reproduced numerically by a gradual increase (in case of filling) or 
decrease (in case of voiding) in p. Such a modeling approach allowed us to follow 
the dynamics of changes in shapes and stress-strain fields in the organ during 
different stages of its function. 

During filling of the bladder, assume that the detrusor remains actively relaxed, 
i.e., no excitatory or inhibitory signals are generated in the tissue. Analysis of the 
total force distribution during this stage indicates that at the beginning, 
p — 1.2 KPa, excessive tensions develop in the longitudinal direction in the body 
and the apex of the organ (Fig. 10.2). A high level of circumferential force is 
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Fig. 10.1 3D MRI image 
reconstruction of the filled 
human urinary bladder 


E RA. 


observed in the region of the flexure that appears along the posterior wall of the 
bladder. It extends bilaterally and upward to the body of the organ. Low-intensity 
longitudinal and circumferential stretch forces are seen in the fundus. With a 
gradual increase in p to 2.7 and 6.5 kPa, respectively, maximum T; and T, values 
are produced in the body and the apex of the bioshell which undergo even biaxial 
tension. A persistent zone of high degree T, is recorded at the flexure. Throughout 
the process of filling all regions of the organ experience biaxial stretching 
with); > 1.0. 

Coronal, sagittal, and transverse views of the three commonly observed ana- 
tomical shapes of the organ for p = 4.0 kPa are shown in Fig. 10.3. Attained 
undeformed geometries demonstrate a high degree of variability and asymmetry. It 
can be expected that the pattern of total force and strain distribution differs 
significantly. For example, in the bladder of the given shape (c), the excessive in- 
plane forces 7|,7, develop along the posterior wall of the body and the apex of the 
organ, while the fundal region remains almost free of tension. Note a smooth 
gradient in 7,1, i.e., the change in intensity extends from the fundus to the apex of 
the bioshell (Figs. 10.4, 10.5). 


10.3.2 Voiding Stage 


The distension of the bladder with urine and activation of intramurally located 
mechanoreceptors lead to “micro-contractions” of the organ. The mechanism is 
attributed to self-excitatory, pacemaker-type reactions which can be detected 
throughout the entire detrusor. 
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(a) i 3 2 
(c) 


Coronal Sagittal Transverse 


Fig. 10.2 Dynamics of changes in the longitudinal, Tj, and circumferential, T., tension fields in 
the bladder during filling 


Let multiple *pacemakers" be randomly distributed over the surface of the 


bioshell. They discharge impulses of average amplitude, V, — 70-80 mV, and 
duration, d = 10 ms, that incite waves of depolarization V; and V., in the lon- 
gitudinal and circular smooth muscle layers of the detrusor. They quickly spread 
over a long distance in the longitudinal direction and over a short distance cir- 
cumferentially (Fig. 10.6). The fronts of the waves may collide with generation of 
a single solitary wave Vi. However, the strength of excitation and poor electrical 
connectivity within the muscle tissue prevent the propagation of Vo; over the 
organ. Only an area that includes the posterior wall of the body, the fundus and the 
apex appear depolarized. 

The activation of L- and T-type Ca?* channels, the influx and rise in intracel- 
lular calcium concentration lead to contractions of the detrusor. Extremely 
localized, patchy areas, mainly encircling the *pacemaker" zones, are subjected to 
weak active forces Te p Which are consistent with experimental observations 
described as “micro-contractions”. 
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Fig. 10.3 Coronal, sagittal, and transverse views of actual configurations of different anatomical 
shapes of the human urinary bladder at p = 1 kPa 


A chain of signal transduction events that involve the pontine micturition 
center, the lumbosacral nuclei, and parasympathetic ganglionic motor neurons is 
marked by widespread long lasting, intense, tonic-type contractions of the 
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Fig. 10.4 Dynamics of changes in axial, &, and circumferential, £<, deformation fields in the 
bladder during filling 


detrusor. The general topology of intercellular arrangements, electrochemical 
coupling and mechanical processes, along with mathematical aspects of the model 
formulation have been discussed earlier. With these numerical experiments, we 
have aimed to reproduce and analyze the following sequence of physiological and 
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Fig. 10.5 Circumferential, A., (top) and axial, Aj, (bottom row) stretch ratio fields in the fully 
filled bladder of the initial shape (c) as in Fig. 10.2. Oblique views of the bladder are shown 


mechanical events that occur during voiding: deformation of free nerve endings of 
mechanoreceptors; generation of APs and their propagation along the unmyeli- 
nated fibers; AP generation at the soma of neurons; electrochemical coupling at 
synapses and production of fast excitatory and inhibitory postsynaptic potentials; 
ion channel dynamics and spread of the wave of depolarization within the detrusor; 
and active force generation and configuration changes of the bladder—bioshell. 

A mechanical stretch of the bladder with urine, € = 0.75, evokes single receptor 
potentials of maximum amplitude, 82.5 mV, and duration —2.3 ms at the 
mechanosensitive afferents of the pelvic nerve. They trigger the development of 
action potentials (APs) of an amplitude 69 mV in the unmyelinated nerve fibers 
that propagate toward the somas of neurons in the lumbosacral nuclei. There they 
initiate high frequency, 8-10 Hz, APs of the amplitude 81.3 mV. The further 
efferent spread of the electrical excitation with subsequent activation of the gan- 
glionic motor neurons and muscarinic cholinergic synapses on smooth muscle 
results in generation of fEPSPs and waves of depolarization V;, Ve. The quanti- 
tative analysis of the processes has been discussed in detail in Chap. 8. 

The waves V., extend rapidly along the muscular syncytium throughout the 
organ. The bladder undergoes a uniform electrical excitation of an average 
intensity of 80 mV. This level of depolarization is strong enough to produce a 
rapid influx of Ca” inside cells. 

A rise in the free cytosolic calcium ion concentration leads to activation of the 
cascade of mechanical reactions with production of the active forces of contrac- 
tion, Te, (Fig. 10.7). Their pattern of propagation resembles the dynamics of 
spread of electrical waves. Weak in strength localized contractions quickly pro- 
liferate throughout different areas of the bladder. The intensity increases in time 
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Fig. 10.6 Dynamics of microcontractions in the fully distended urinary bladder 


and the whole organ becomes uniformly contracted—a condition that is essential 
for effective emptying of the organ. Note that there are no changes in the actual 
configuration of the bioshell despite strong active forces being present. The 
intravesicular pressure rises to ~ 10 kPa. 
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RCC 


Fig. 10.7 Dynamics of contractions of the detrusor as a result of activation of the micturition 
regulatory mechanisms 


The relaxation of the internal urinary sphincter and emptying of the bladder, 
which has been achieved in the model by decreasing p values, results in gradual 
changes in shape and the reduction in dimensions of the organ (Fig. 10.8). The 
bioshell attains homeomorphic configurations during the process. 
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Fig. 10.8 Changes in configurations of the bladder during voiding 


10.3.3 Pharmacological Modulation of Voiding 


There are a limited number of therapeutics available currently to treat bladder 
dysmotility. The leading edge of therapy relies on selective and nonselective 
muscarinic 44.3 receptor antagonists, e.g. oxybutynin, darifenacin, fesoterodine, 
etc., and, to a lesser extent, f; 4 adrenergic receptor agonists, e.g., terbutaline, 
clenbuterol. We have studied numerically their effects on the fully distended 
urinary organ. 

The application of cholinergic antagonists has a negative effect on the dynamics 
of signal transduction processes both at the neuro-neuronal synapses and at the 
detrusor, by significantly reducing the amplitude of fEPSPs and APs, respectively, 
and inhibiting depolarization of the muscle membrane. Adrenergic receptor ago- 
nists facilitate generation of fIPSPs postsynaptically, and thus inhibit the spread of 
excitation along the parasympathetic pathways and depolarization of the detrusor. 
As a result, the wall of the bladder remains relaxed or produces low-intensity 
active forces of contraction Tj, ~ 0. Remarkably, the strength of total in-plane 
forces, T, and T,, increase throughout the bioshell while its shape does not change. 

These results are not unexpected and could have been easily predicted from the 
reductionist analysis of the neuronal arrangements and the detrusor that has been 
done separately. The mysteries of bladder function are concealed in the multi- 
functional integration of its morphostructural elements. Although the extension of 
the current model to simulate: (1) co-localization and co-transmission by multiple 
neurotransmitters, (2) receptor polymodality, (3) receptor distribution, (4) 
mechanical heterogeneity, and (5) detailed neuromorphological (intraganglionic) 
ensembles is beyond the scope of this book, it is important to recognize their 
implications on the functionality of the organ. 


Chapter 11 
Challenges in Human Urinary Bladder 
Mechanics 


You will never do anything in this world without courage. 
Aristotle 


11.1 Urine-Bladder Interaction 


The study of the coupling between the motion of the bladder wall and urine leads 
to a mathematical and numerical problem whose complexity is enormous, but 
within the realm of today’s computational capabilities. The models proposed are 
based on simplified assumptions which allow one to capture only a few aspects of 
the physics of events and utilize special structural features, such as geometry, 
symmetry, and periodicity. Currently, there are no analytical results proving the 
well-posedness of the mathematical problem of urine-bladder interaction. 

Throughout the book, we have focused our attention on the mechanics of the 
wall of the urinary bladder, regulatory mechanisms and pharmacological effects of 
drugs on its contractility. Although we have assumed in our modeling and 
numerical simulations that urine is implicitly present, it has been broadly 
embodied by a single parameter-intravesicular pressure p. Such approximation has 
shown to be adequate in the analysis of shape changes and stress-strain distribution 
in the organ under different stages of filling and voiding. However, the explicit 
inclusion of urine as a second phase becomes essential in problems related to the 
quantification, assessment, diagnosis, and establishment of the efficacious treat- 
ment of patients with abnormal urodynamics. In the following paragraphs, we 
provide a brief coupled urine-bladder mathematical formulation of the problem to 
reflect the major physiological and pathological phenomena observed clinically 
and experimentally. 

A medium where the stress vector p, at any point VM € S coincides with the 
direction of the normal 7 to the surface S is called an ideal fluid. The tensor surface 
in that case is a sphere and p; = p» = ps, i.e., the principal components of the 
stress tensor are equal. We shall denote p; — p» — p3 — —p, where p is the 
(mechanical) pressure. 

A fluid for which the stress tensor is given by 


p! = —pg! +7! (i,j = 1,2,3), (11.1) 
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where, 
p> P (Pus Ty, se adis 
A = fÜ (èag, 81, X1, = ibus 


is called a viscous fluid. Here, p, is the density of fluid, T is temperature, t(tË) is 
the shear stress exerted by the fluid, £5 are covariant components of the stress and 
of the rate of deformation tensors, respectively, g are contravariant components of 
the metric tensor, and 7,,..., y, are physicochemical parameters, e.g., the fraction 
of solutes in urine. 

Assume that f? can be expanded in a Taylor series with respect to é,g and the 
fluid satisfies the condition: p = 0— é,; = 0. Then, for tË, neglecting nonlinear 
terms in the expansion, we obtain 


tË = AYP np, (11.2) 


which is the general form of Navier-Stokes’s law for viscous Newtonian fluids. 
The total number of components of the tensor of intrinsic mechanical properties 
A¥8 is 81. However, taking into consideration the symmetry of the stress and the 
rate of deformation tensors, p/ = př, £5 = Ega, the number of AU"? can be reduced 
to 36. In case of isotropic fluid, only nine non-zero components remain: 
AH —A?22 = A993. AU —A 1133 = A93. (Alia =A, AU £ 0, i £z j). Then 
introducing parameters 4, 4,, such that A = 2u, + 4,, Al” = Au, expression 
(11.2) can be written in the principal axes of the stress and the rate of deformation 
tensors, é and t(,jj, as 


3 
ti = Žu X èj + 2,8), i= 1,2,3. (11.3) 
j=1 


Multiplying (11.3) by (dx y and after simple rearrangements (Sedoy 1976) for 
Navier-Stokes’ law in the general curvilinear coordinate system, we have 


toj = Auli (E) Bi + 2Hutij, (11.4) 
and for contravariant components 


u 


Substituting (11.5) into (11.1) for an isotropic viscous fluid in the general 
curvilinear coordinate system, we get 


p! = —pg! + Aug divo, + 2,2", (11.6) 


where J; (é) : = divo, and 0,(uUy,Vu,W) is the velocity vector. 
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The fluid is considered to be incompressible if an arbitrary volume of it remains 
constant throughout the space and time. Then, from the law of conservation of 
mass, we obtain 


dp, Op, 
dí ôt 


TU Vp, — 


+V.: (Puu) = 0, (11.7) 


where dp„dt is the material derivative and the meaning of other parameters is as 
described above. In the case of homogenous (p, = const at any point) incom- 
pressible fluid, p,, becomes a single-valued parameter and (11.7) takes the simplest 
form 


V.D, — 0, (11.8) 


which is the statement of conservation of volume. 

To derive the system of Navier-Stokes equations that describe the dynamics of 
viscous incompressible Newtonian fluid, we proceed with the law of motion. In the 
general curvilinear coordinate system, it is given by 


pā = p.f. + Vip! (11.9) 
where à is acceleration 
. 0 Ov! f) 
qi = n D= +o (C40). (11.10) 


Here P is the density of the mass force, V; is the sign of the covariant 
derivative, and T d are the Christoffel symbols of the second kind. Covariant 
differentiation of (11.6) with respect to pË, yields 


Vip! = -Vip F (Au F 1,)V'(V i Du) + ju, Av’, 


where V', A are the contravariant derivative and the Laplace operator, respec- 
tively. In the above, we have assumed that 1, and u, are constants. Additionally, if 
the fluid is incompressible then 


Vip! = -Vip + u, Av. (11.11) 
Substituting (11.11) in (11.9), we obtain 


p.a! = p.f! — Vp + p, Av. (11.12) 


Equations (11.12), (11.8) constitute the governing system of equations of motion 
of viscous incompressible Newtonian fluid in the general curvilinear coordinate 
system with constant dynamic viscosity u„. The formulation of pertinent initial and 
boundary conditions closes the problem. 

In cylindrical coordinates {r, 9, z]. (11.12) and (11.8) take the form (Sedov 
1976) 
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205; Op u, 20w 

pur =F,- P+ m (Au - RARE). 

MEET 1 ôp Uo 2 Ou; 

Pile =F roe ty( Ay ~ 28-4 55"), (11.13) 
_ = Op 

Puz =f, ag z” H, Au; 

and 

Op, , OPullr?  OPuYo , ÕPuWz _ 

ro E= + 0o +r Dz = 0. (11.14) 


Here, the Laplacian is 


9 o 101 oo 
Oo? rr r0g9 O2 


and the components of the acceleration vector are 


|. Ou, Ou, | vo Ou, Ou, ve 
ep S ui r’ 
Ovo Ovo | Vg Ov Ovo Urvo 
E E! 4 , 
te= a ar r 0p i Oz r 
_ ôw, Ow, vo Ow; Ow; _ 
az = CTS Ur E F dp T Wz ae Du(Ur, Vo, Wz) - 


Urine is approximately 95 % water with the remainder 5 % of dissolved non- 
organic and organic solutes. These include Na*,K^, Ca?*, Mg^*, Cl, 
ammonium, phosphate, and sulfate ions, along with urea, creatinine, and uric acid. 
Their concentration depends on metabolic state, diet, and administered medica- 
tions. Normal urine density is p, ~ 1.003—1.035 (g/cm?), while that of water is, 
Pw ~ 0.998 (g/cm?). Any deviations are associated with urinary disorders, 
e.g., Py < 1.003 g/cm? is observed in overhydration, antidiuretic hormone defi- 
ciency, and the loss of ability by the kidneys to concentrate urine; 
Pu > 1.035 g/cm? is commonly seen with albuminuria, glucosuria, pyuria, etc. 
These variations are linearly related to the solute weight fraction, [...], and are 
given by (Putnam 1971) 


p, = 0.04775]. . .] + 0.99325. (11.15) 


The estimated dynamic viscosity, u,, of urine at normal body temperature 
T = 37 °C, is 0.6531 < u„ < 0.7978 (mPa-s) and is within the range of that of 
water, 0.548 < uy, < 1.002 (mPa-s). Experiments have also revealed that yu, is 
sensitive to the variation of the fraction of solutes 
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0.9 exp(1.5 14h) for [...]<0.5 


1.42 exp(1.25 EL) for [...] > 0.5. 


Hy = (11.16) 


Although the dynamic viscosity of fluids changes with temperature, the human 
body operates within 35—40.5 °C. The decrease or increase in T above the phys- 
iological limits is seen in extreme pathological conditions—severe hypothermia, 
malignant hyperthermia—and is incompatible with life. Additionally, the kine- 
matic viscosity, v,, (v, = u/p,), varies between 0.6518 < v,< 0.795 and 
V, ~ Vy. From this analysis, it is obvious that urine can be treated as a viscous 
incompressible Newtonian fluid. 

Consider the urinary bladder to be containing urine, i.e., a soft bioshell filled 
with viscous incompressible Newtonian fluid. The system of (10.13, 10.14, 10.19, 
10.21, 11.13 and 11.14), complemented by auxiliary relationships—constitutive 
relations for the tissue, ion currents, the activation and deactivation parameters of 
ion gates—along with initial and boundary conditions, comprises the mathematical 
formulation of the problem of urine-bladder interaction. Only a combined solution 
to the system of equations will offer a missing insight into and help to unravel the 
intricate processes of signal transduction in neuronal assemblies, generation and 
propagation of myoelectrical waves in the detrusor, mechanical responses with 
changes in stress-strain states, and configurations and the flow of urine in the 
bladder. At the moment when the book was written, no solutions to the problem of 
combined analysis of urinary flow in the biologically active bladder were 
available. 


11.2 Brain-Bladder Axis 


In research we like to assume a linear, causal relationship between two or more 
factors. This is a valid way of testing and exploring spatiotemporal organization 
and pattern formation in the simplest, but not real, biological setting. The medical 
practice in the West has, for more than a century, been dominated by a dualistic 
view of human nature, i.e., biological and psychological, and reflects the reduc- 
tionist methodology to study complex phenomena. Today, science is moving into a 
position of integration. The systems biology approach is developing rapidly and is 
the subject of debate at scientific forums and meetings. The main argument of this 
discussion is the interplay between psychological stresses and somatic symptoms 
in the development of diseases. It has become convincingly clear that one cannot 
experience an emotion or think a thought without biological correlates. Unresolved 
mental conflicts lead to activation of the central nervous system, the autonomic 
system, and, finally, the organism as a whole. 

In this respect, some progress has been made in gastroenterology. The brain-gut 
axis relationship is being entertained as a theoretical explanation to the patho- 
genesis of irritable bowel syndrome (Mayer and Tillisch 2011; Wilhelmsen 2000; 


168 11 Challenges in Human Urinary Bladder Mechanics 


Camilleri and Di Lorenzo 2012; Prins 2011; Aerssens et al. 2007). It stipulates that 
stress signals from the gut pass to the designated areas in the brain, where they are 
integrated with other interoceptive signals from the body and contextual infor- 
mation from the environment. The integrated response is sent back to various 
target cells within the gut to assure its homeostasis during external/internal per- 
turbations and to adapt its function to the overall state of the body. 

The brain-urinary bladder interactions in overactive urinary bladder, nocturnal 
incontinence, painful bladder, and detrusor hyperreflexia, among others, have 
never been addressed in the scientific literature, although the associations appear to 
be obvious. Questions related to the pathogenesis of visceral (bladder) hypersen- 
sitivity, either chemical or mechanical in nature, its neuropathways and molecular 
mechanisms have never yet been studied in detail. There is experimental evidence 
that a dense nexus of suburothelial nerves located in the neck of the bladder are 
critical in the sensory function of the urothelium (Birder et al. 2001; Wiseman 
et al. 2002). Traditionally viewed merely as a passive barrier, it has now become 
clear that the urothelium vividly responds to exogenous/endogenous stimuli and 
transmits them within the bladder wall (de Groat 2004). 

Nociceptors are sensory end organs that selectively respond to noxious or 
potentially tissue-damaging stimuli. Surprisingly, little is known about their ana- 
tomical distribution within the wall of the bladder, biochemical heterogeneity, 
mechanical, electrophysiological, and pharmacological properties, although, these 
features underlie different aspects of pain. Populations of nociceptors are often 
defined on the basis of receptortype expression. Cholinergic, adrenergic, brady- 
kinin, purinergic, and vanilloid receptors, along with their respective chemical 
mediators, e.g. ACh, NA, ATP, and NO, have been confirmed experimentally in 
the suburothelium. However, facts about their co-localization and co-transmission 
remain an elusive, if not unobtainable, objective. 

In functional bladder disorders, the major subjective complaints include dis- 
comfort and pain in the suprapubic area. The following chain of events can be 
offered as an explanation: chemical, mechanical, or thermal activation of the 
intramural nociceptors results in low-amplitude signals which are transmitted 
through the pelvic nerves to the spinal cord where they are integrated to produce 
the concurrent excitation is passed to the periaqueductal gray in the brain via the 
interneurons in the lumbosacral spinal cord. The presumed connection between 
central and lateral PAG enables them to control excitation of the PMC. Thus, 
during normal bladder filling, PMC neurons and the parasympathetic pathway are 
turned off; but at a critical level of sensitization' and visceral hypersensitivity,” the 
afferent signaling switches them to maximal activity. 

In the brain, many neuron populations are distributed through different nuclei 
and centers: the medullary raphe nuclei, the locus coeruleus, the brain stem, 


! Sensitization is defined as increased reactivity to stimuli in pain pathways. 


? Visceral hypersensitivity is the exaggerated experience of pain in response to mildly painful or 
even normal visceral stimuli. 
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Barrington’s nucleus, the caudal and preoptic hypothalamus, and the cerebral 
cortex are all involved in regulation of lower urinary tract function. Psychological 
factors influence the synaptic mechanisms and feed forward loops from the brain. 

It has been suggested that the descending central pathways modulate the 
transmission of nociceptive information at the spinal cord level and lead to 
somatization.^ Sensitization occurs in the dorsal horn of the spinal cord, and it may 
also originate at the level of the limbic structures. At the moment, the exact 
mechanisms and pathways engaged are unclear. 

Existing investigative techniques are also rather limited. The use of imaging 
studies such as positron emission tomography and magnetic resonance imaging is 
still in its infancy in urology and little quantitative information is available 
regarding the central nervous system activation in functional bladder disorders. 
Unraveling the mysteries of the brain-bladder interaction will provide novel 
treatment targets and therapies to patients suffering from these debilitating 
conditions. 


11.3 Bladder Substitutes 


Bladder replacement and reconstruction remain the mainstream modalities in 
treatment of bladder disorders arising from bladder cancer, exstrophy, neurogenic 
dysfunction, chronic inflammatory diseases of the bladder such as interstitial 
cystitis, tuberculosis and schistosomiasis, and ketamine abuse, among others. 
Surgical options include orthotopic conduit, continent diversion, and bladder 
replacement. All existing procedures—the creation of intestinal conduits or the 
neo-bladder—are aimed to enhance the reservoir capacity of the organ formed. 
However, none of them offer the required functionality—controlled contractility— 
and thereby do not fulfill the need of regulated inflow and outflow of urine. 
Questions related to surgical procedures, short- and long-term complications are 
beyond the scope of this book. Here, we shall focus our attention only on problems 
that are pertinent to biomechanics of the organ. 

Currently, the most reasonable alternative substitutes for the bladder are seg- 
ments of the small and large intestines and tissues grown in vitro. The search 
continues for successful biological and artificial grafts, which replicate myoelec- 
trical properties of the detrusor and, thus, respond effectively to the neuromodu- 
latory signals in the body. Dedicated comparative studies of smooth muscle layers 
of the ileal and recto-sigmoid regions of the gastrointestinal tract and the detrusor 
have demonstrated similarities in electrophysiological and pharmacological reac- 
tions to electrical, mechanical, and chemical stimulations. Functional correlations 
that are advantageous for finding the biomechanically fair construct include the 


? Somatization is a descriptive term meaning somatic complaints not fully explained by any 
known medical condition. 
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development of phasic contractions in response to treatment with ACh, and 
relaxation to NA, the lack of spontaneous and/or baseline electromechanical 
activity in the absence of interstitial cells of Cajal, and the ability of the circular 
smooth muscle layer to generate tonic-type contractions. It should be borne in 
mind that the performance of the intestinal conduit depends mainly on the function 
of the intramural ganglionic network of cells rather than the extrinsic innervation. 
While the enteric nervous system is well developed and can operate independently 
from the central nervous system, little is known about the capabilities and prop- 
erties of the autonomic nervous system of the bladder. Investigations into its 
function are very fragmentary and descriptive. The unresolved question of pace- 
maker activity and the paucity of morphophysiological data on intraganglionic 
arrangements and signal transduction processes within the bladder wall hamper the 
development of a successful artificial pacemaker and the coordination of filling and 
micturition of the reconstructed bladder. 

Active (contractile) and passive properties of the substitute tissue present 
another hurdle in reconstruction of the bladder. We have only recently started 
gaining insight into uniaxial and biaxial biomechanical characteristics of the 
gastrointestinal tract and the bladder. Most experiments are being conducted on 
animal tissues and only a few studies are related to human organs (Alexander 
1971, 1976; Andersson et al. 1989; Bider et al. 2001; Coolsaet et al. 1975a,b, 1976; 
Egorov et Al 2002; Finkbeiner 1999; Gloeckner 2003; Gloeckner et al. 2002; 
Kondo and Susset 1973; Kondo et al. 1972; Korossis et al. 2009; Miftahof and 
Nam 2010; Nagatomi et al. 2008; Parekh et al. 2010; Sacks 2000; van Mastright 
and van Duyl 1978; van Mastrigt and Nagtegaal 1981; Venegas et al. 1991; Wagg 
and Fry 1999; Wognum 2010). Results demonstrate nonlinear force-ratio of 
deformation relationships, however, with a great degree of variability in response 
to dynamic loading. Thus, the intestinal tissue is softer, i.e., it stretches by 
20-30 % of its original length at T.; c 0, compared to the bladder wall that 
extends by 5-10 % at T., ~ 0; the intensity of active forces generated by the 
detrusor in response to electrical and chemical stimulations exceed those produced 
by the intestinal smooth muscle; the intestinal wall withstands a maximum stress 
of 0.6-1.5 kN/cm?, and the bladder 1.0-4.1 kN/cm?. This explains the current 
success in actual surgical recreation of the neo-bladder (pouch) which serves as a 
reservoir for urine but with impaired or little functionality. Thus, patients do not 
appreciate its state of fullness and have to squeeze the "bladder" to drain urine. 

Biomimetics is a term frequently used to describe the use of concepts and 
principles from nature and their application to creating new materials, devices, and 
systems. Investigative methods of biomimetics have uncovered many of the fun- 
damental biological processes of self-assembly, cellular mechanisms, and inte- 
grative biology. Key advances have also been made in the ability to engineer and 
grow tissues and organs, including the urinary bladder, and the design of the 
"smart bladder pacemaker". The progress in this field in the last decade has been 
remarkable and is promising, although more research is needed before these 
innovative discoveries can be used in clinical practice. 
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The subject of the book is amazingly rich and broad, and I offer my apologies to 
those whose favorite topics have been omitted. However, I hope that readers will 
have been sufficiently inspired to dedicate themselves further to this fascinating 
area for research. 
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